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Abstract  of  Dissertation  Presented  to  the  Graduate  School 
of  the  University  of  Florida  in  Partial  Fulfillment  of 
the  Requirements  for  the  Degree  of  Doctor  of  Philosophy 

PHYSICAL  MODELING  OF  MO S- CONTROLLED  HIGH-VOLTAGE 
DEVICES  FOR  INTEGRATED  CIRCUIT  COMPUTER-AIDED  DESIGN 

By 

YEONG-SEUK  KIM 
August  1990 

Chairman:  Dr.  J.  G.  Fossum 

Major  Department:  Electrical  Engineering 

This  dissertation  presents  methodology  for  physical  charge- 
based  modeling  of  MOS-controlled  high-voltage  (HV)  devices 
for  integrated  circuit  (HVIC)  computer-aided  design  (CAD) . 
New  models  for  two  different  MOS-controlled  HV  devices, 
the  insulated-gate-bipolar  transistor  ( IGBT)  and  the  double- 
diffused  MOS  transistor  (DMOST) , are  developed.  The  effects 
of  the  buffer  layer  and  static/dynamic  latch-up  in  the 
IGBT  are  characterized,  and  quasi-saturation,  space-charge- 
limited  current  flow,  and  effects  of  the  inherent  BJT 
in  both  vertical  and  lateral  DMOSTs  are  modeled.  These 
effects,  which  are  not  properly  represented  in  conventional 
equivalent- (sub) circuit  models,  are  physically  and  sometimes 
semi -numerically  accounted  for  in  our  models.  Two-dimensional 
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numerical  device  simulations  were  used  extensively  to  study 
the  effects  and  to  aid  the  model  development. 

The  developed  models  are  implemented  in  the  circuit 
simulator  SPICE  via  FORTRAN  subroutines  (UDCSs) . With  only 
known  structural  (device)  parameters  and  crudely  extracted 
model  parameters,  device/circuit  simulations  favorably 
predicted  measured  characteristics  of  test  devices . The  models 
in  SPICE  provide  a capability  of  mixed-mode  device/circuit 
simulation,  which  is  not  afforded  by  other  equivalent- 
(sub) circuit  models,  and  hence  can  facilitate  computer-aided 
optimal  device/circuit  design  of  HVICs . 
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CHAPTER  1 
INTRODUCTION 

Recent  advances  in  high-voltage  (HV)  device  technology 
have  enabled  the  integration  of  HV  devices  with  low-voltage 
control  circuitry  on  the  same  chip.  Such  high-voltage 
integrated  circuits  (HVICs)  now  commonly  use  MOS-controlled 
HV  devices,  exploiting  the  high  input  impedance  of  the  MOS 
gate  relative  to  that  of  power  bipolar- junction  transistors 
(BJTs)  . The  MOS  channel  is  usually  formed  by  a double- 
diffusion process;  sequential  diffusion  of  p-  and  n-type 
impurities  in  the  epi  (or  substrate)  region  yields  the  MOS 
channel  as  shown  in  Fig.  1.1.  The  double-diffusion  process  is 
widely  used  for  HV  devices  because  of  easier  process  steps 
to  achieve  short  channels  while  realizing  high  breakdown 
voltages.  MOS-controlled  HV  devices  built  with  the  double- 
diffusion process  include  vertical  and  lateral  insulated- 
gate-bipolar  transistors  (IGBTs)  and  double-diffused  MOS 
transistors  (DMOSTs) . The  difference  between  these  two 
structures  is  due  to  the  anode  (drain)  design;  the  IGBT 
is  a modification  of  the  DMO ST  structure  made  by  forming  a 
pn  junction  at  the  drain  contact  as  shown  in  Fig.  1.1  (for 
n-channel  devices) . 
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Gate 


Drain/ 


Fig.  1.1  Cross-sectional  view  of  the  IGBT  (p+— diffusion  in 
the  anode)  and  the  DMOST  (n+— diffusion  in  the 
drain) . Typical  gate  oxide  thickness  and  channel 
length  are  0.1  //m  and  2 //m,  respectively. 
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Along  with  the  matured  process  technology  of  HV  devices  and 
integrated  circuits,  CAD  simulation  tools  should  be  developed 
to  support  optimal  device  and  circuit  design.  The  equivalent - 
(sub) circuit  SPICE  models  for  HV  devices  that  are  used  today 
are  too  empirical  and  do  not  properly  account  for  the  unique 
behavior  of  the  HV  devices,  such  as  latch-up  in  the  IGBT  and 
quasi-saturation  and  space-charge-limited  current  flow  in  the 
DMOST . Therefore  the  development  of  physically  based  models 
for  MOS-controlled  HV  devices,  and  implementation  into  the 
simulation  tool,  e.g.,  SPICE,  are  needed. 

The  purpose  of  this  dissertation  is  to  develop  and 
implement  physical  models  of  MOS-controlled  HV  devices,  i.e., 
IGBTs  and  DMOSTs,  for  HVIC  CAD.  The  effects  of  the  buffer  layer 
and  static/dynamic  latch-up  in  the  IGBT  are  characterized, 
and  quasi-saturation,  space-charge-limited  current  flow,  and 
effects  of  the  inherent  BJT  in  both  vertical  and  lateral 
DMOSTs  are  modeled.  The  developed  models  are  implemented  in 
the  circuit  simulator  SPICE  via  FORTRAN  subroutines  (UDCSs) . 
The  major  contributions  made  in  this  dissertation  are 

(1)  the  formulation  and  solution  of  the  steady— state 
and  transient  carrier  transport  problems  in  the  IGBT 
including  a buffer  layer; 

(2)  the  study  of  static  and  dynamic  latch-up  in  the  IGBT 
based  on  two-dimensional  numerical  device  simulations 
and  the  development  of  models  for  circuit  simulation; 
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(3)  the  development  and  experimental  verification  of 
physical  charge-based  models  for  both  vertical  and 
lateral  DMOSTs; 

(4)  the  implementation  of  the  IGBT  and  DMO ST  models  into 
SPICE  using  user-defined  controlled  sources  (UDCSs) . 

In  Chapter  2,  a theoretical-experimental  study  of  the 
effects  of  the  buffer  layer  on  steady-state  and  transient 
vertical  IGBT  (VIGBT)  performance  is  presented.  Results 
indicate  that  previous  studies  are  incomplete  and  hence 
misleading,  and  that  in  fact  optimal  design  criteria  for 
the  buffer  layer  must  be  defined  based  on  more  comprehensive 
analyses.  The  study  gives  physical  insight,  needed  in  such 
definition,  and  describes  proper  model  parameter-extraction 
techniques  not  previously  considered. 

In  Chapter  3,  an  insightful  study  of  static  and  dynamic 
latch-up  in  the  lateral  IGBT  (LIGBT) , based  on  extensive  two- 
dimensional  numerical  device  simulations,  is  described.  The 
insight  is  then  used  as  a basis  for  developing  a physical  SPICE 
LIGBT  model,  which  is  useful  for  device  simulation  and  design 
(e.g.,  for  latch-up  immunity)  as  well  as  for  HVIC  simulation. 
The  basic  mechanisms  underlying  latch-up  in  the  LIGBT  are 
identified  for  various  excitations.  The  significance  of  the 
effective  pin  diode  that  materializes  via  the  conductivity 
modulation  in  the  regenerative  process  is  stressed  in  the 
simulations,  as  is  the  importance  of  non-quasi-static  bipolar 
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transistor  behavior,  heretofore  unrecognized.  The  SPICE  model 
is  supported  by  measurements  of  test  devices  and  by  the 
numerical  device  simulations. 

In  Chapter  4,  a physical,  semi-numerical  charge-based  model 
for  vertical  DMOSTs  (VDMOSTs)  is  developed  and  implemented 
in  SPICE.  The  model  is  derived  from  regional  quasi-static 
analyses  of  carrier  transport  which  implicitly  characterize 
the  device  currents  and  charges  and  which  require  numerical 
solution.  Newton-Raphson  iterative  device  solutions  are 
derived  within  the  circuit  nodal  analysis  framework  of  SPICE. 
PISCES  simulations  and  measurements  of  test  devices  support 
the  model,  which  is  demonstrated  in  dc  and  transient  SPICE 
simulations  of  VDMOSTs  and  HVICs . 

In  Chapter  5,  two  different  structures  of  the  lateral  DMO ST 
(LDMOST) , i.e.,  the  LDD  and  RESURF  LDMOSTs , are  studied 
extensively  using  the  two-dimensional  device  simulator 
PISCES.  The  PISCES  simulations  provide  new  physical  insights 
on  the  normal-  and  reverse-mode  operations  of  the  LDMOSTs, 
which  are  used  for  developing  the  composite  LDD  LDMOST  model. 
In  the  modeling  methodology,  the  LDD  LDMOST  is  regionally 
partitioned  into  three  main  components  (the  channel,  the  drift 
region,  and  the  npn  BJT) , and  carrier-transport  problems  in 
each  region  are  solved.  The  composite  model  is  implemented 
in  SPICE  for  HVIC  CAD  and  is  supported  by  measurements.  The 
modeling  methodology  is  also  applicable  to  the  RESURF  LDMOST. 
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In  Chapter  6,  the  main  accomplishments  of  this  dissertation 
are  summarized,  and  suggestions  for  further  research  are 
discussed. 

Appendix  A,  related  to  Chapter  3,  describes  the  algebraic 
manipulations  of  a base  charge  partitioning  to  represent  non- 
quasi-static  behavior  in  the  bipolar  transistor  of  the  IGBT . 

Appendix  B lists  the  UDCS  source  code  of  the  LDMOST  SPICE 
model  developed  in  Chapter  5 . 


CHAPTER  2 

EFFECTS  OF  THE  BUFFER  LAYER  ON  VIGBT  PERFORMANCE 

2 . 1 Introduction 

This  chapter  presents  a buffer-layer  model  for  the  vertical 
IGBT  (VIGBT)  [Bali84] . The  model,  based  on  the  one-dimensional 
representation  of  the  device,  accounts  for  buffer-layer 
effects  on  both  steady-state  and  transient  characteristics 
of  the  VIGBT.  The  carrier-transport  (continuity/ambipolar) 
equations  in  the  quasi-neutral  regions  of  the  device  are 
solved  to  yield  the  dc  current-voltage  dependence  I.\(V\)  and 
the  turn-off  transient  characteristics  I\(t)  of  the  VIGBT.  Some 
of  the  modeling  is  applicable  as  well  to  the  lateral  IGBT 
with  a buffer  layer,  which  is  discussed  in  the  next  chapter. 

The  VIGBT,  also  called  the  COMFET  [Rus83],  has  recently 
emerged  as  a useful  high-voltage  power  switch.  It  is  a merged 
MOS/bipolar  structure  that  exploits  the  advantages  of  MOS  and 
bipolar  devices  to  yield  low  forward  (ON)  voltage  drop,  fast 
switching  time,  and  high  input  impedance. 

The  basic  structure  of  the  vertical  n-channel  IGBT,  shown 
in  Fig.  2.1,  has  a DMO ST  merged  with  a wide-base  pnp  B JT . The 
n+— buffer  layer  [Goo83]  between  the  n~— base  and  the  p+-emitter 
is  commonly  included  to  prevent  punch-through  in  the  forward- 
blocking mode  and  indeed  to  enhance  the  device  performance 
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Fig.  2.1  Unit-cell  structure  of  vertical  n-channel  IGBT  with 
n+— buffer  layer.  The  external  circuitry  and  gate 
bias  indicated  correspond  to  a transient  turn-off 
from  the  steady-ON-state . 
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[Goo83,  Nak85] . Such  enhancement  has  been  experimentally 
demonstrated  as  an  increase  in  the  (static)  current  required 
to  latch  the  pnpn  structure  [Goo83]  and  as  an  improvement  in 
the  trade-off  between  the  steady-ON-state  voltage  drop  and 
the  turn-off  time  [Nak85] . A recent  theoretical  study  [Kuo86] 
implied  design  criteria  for  the  epitaxial  (n-  and  n+)  layers 
of  the  device  to  effect  this  trade-off.  The  study  was  not 
sufficiently  rigorous  however  to  guarantee  optimal  designs, 
as  demonstrated  in  this  chapter. 

A comprehensive  theoretical  study  [Hef86]  of  buffer- 
layer  effects  on  VIGBT  performance,  predicting  that  the 
buffer  layer  decreases  the  switching  energy  loss,  was 
recently  presented.  This  study  involved  much  mathematical 
detail  but  little  experimental  support,  and  hence  its 
results  cannot  be  put  into  proper  perspective  because  the 
relative  significance  of  various  physical  mechanisms  cannot 
be  ascertained.  In  this  chapter  we  use  a similar,  but 
simpler  and  more  insightful  model  to  fully  explain  performance 
data,  including  new  experimental  results  given  herein  that 
seem  to  contradict  [Nak85,  Kuo86]  and  that  suggest  optimal 
design  criteria  for  the  buffer  layer.  The  simple  analytic 
model  accounts  for  buffer-layer  effects  on  both  steady-state 
and  transient  characteristics  of  the  VIGBT.  We  support  the 
model  expermentally , introducing  useful  parameter-extraction 
techniques,  and  we  discuss  its  utility  alluded  to  above, 
explaining  the  seemingly  contradictory  data. 
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2.2  VIGBT  Model 

Our  study  of  the  buffer  layer  is  done  effectively 
using  a simple  model  to  describe  both  steady-state  and 
transient  characteristics  of  the  VIGBT.  The  model  is  based 
on  a regional  analysis  [Hef86,  McD85,  Fos86]  of  the  one- 
dimensional representation  of  the  vertical  p+n+n“p  BJT  in  the 
VIGBT  illustrated  in  Fig.  2.2. 

2.2.1  Steady-ON-State 

In  the  steady-ON-state,  a voltage  is  applied  as  shown 
between  the  emitter  (anode)  and  collector  (cathode)  terminals, 
and  the  base  is  driven  by  the  DMOST,  the  gate  of  which 
is  biased  above  the  threshold  voltage  to  enable  ON-state 
conduction. Consistent  with  typical  operation,  we  assume  high- 
level  injection  (p  ~ n)  in  the  quasi-neutral  n~— base  region, 
and  low-level  injection  in  the  n+-buffer  and  p+-emitter 
regions.  We  further  neglect  recombination  in  the  relatively 
narrow  buffer  layer,  an  assumption  which  we  have  found  to  be 
reasonable  for  typical  devices.  We  stress  that  the  primary 
effect  of  the  buffer  layer  on  the  device  performance  is 
manifested  through  the  carrier  transport  across  the  layer, 
irrespective  of  any  recombination  in  the  layer.  Additional 
simplification  of  the  steady-state  analysis  is  effected  by 
noting  that  the  BJT  is  constrained  to  operate  in  (near)  the 
forward-active  mode,  with  negligible  base-width  modulation, 
since  the  voltage  drop  across  the  DMOST  is  small  [Fos86] . 


11 


x=o  wf+w  f- 


p- 

ri n 

i i 

i i 

i i 

n+  j 

i i 

i i 

i i 

i i 

n- 

T 

1 

1 

1 

P 

( 

'A 

(Emitter)  | 

i 

i 

i 

j 

! (Butter)  | 

i i 

i i 

i i 

i i 

(Base) 

i i 

i i 

i i 

i i 

[(Collector) 

i 

i 

i 

J 

U — wt— ►! 

◄ WB ► 

2.2  One-dimensional  representation  of  vertical  p+n+n~p 
BJT  in  the  VIGBT.  The  junction  space-charge  regions 
are  indicated  by  the  dashed  lines . 


Fig. 


12 


Using  a regional  modeling  approach  (region  boundaries  are 
defined  in  Fig. 2.2),  we  first  solve  the  steady-state  ambipolar 
transport  equation  in  the  quasi-neutral  n“-region,  obtaining 
[Fos88 ] 


p(.r) 


■ ; / Wf  +H'”— * 

smn  I — L-j-[ 


sinl‘  (rr) 


(2.1) 


where  L \ is  the  ambipolar  diffusion  length  and  TF  = TT'/j  - is 
the  quasi-neutral  base  width.  Note  in  (2.1)  that  the  carrier 
density  drops  to  zero  at  the  edge  of  the  base-collector 
depletion  region  (x  = W~  + IT")  . The  analogous  solution  in  the 
quasi-neutral  n+-buffer  layer,  for  a uniform  doping  density 
and  negligible  recombination,  is 

p(»?)  -p(«) 

i,|x|  = ifc *+P«>)  • (2.2) 


In  (2.2)  we  have  used  the  fact  that  the  width  of  the 
n*n  space-charge  region  is  much  narrower  than  the  buffer- 
layer  width  ir ,,  but  have  properly  accounted  the  significant 
difference  between  p (ir+)  and  p(\Vj).  In  fact,  for  quasi- 
equilibrium  in  the  space-charge  region, 

p{Wf)NDf=  [P(WF)Y  (2.3) 
where  NDf  is  the  buffer- layer  doping  density. 

To  derive  the  dc  current -voltage  dependence  IA(VA)  of  the 
VIGBT  and  relate  it  to  the  buffer-layer  properties,  we  must 
connect  the  boundary  hole  densities,  p(  0),  p (V+) , p (uj)  in 
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(2.1)  and  (2.2),  to  the  terminal  conditions.  To  make  this 
connection,  we  first  express  the  hole  current  density  in  the 
buffer  layer  as 


Jp  (o  < .r  < ir+ 


) - ~<iDvr 


p(h 7)  -pm 


w 


(2.4) 


where  D,,/  is  the  buffer-layer  hole  diffusivity.  The  electron 
current  density  at  x = 0 is,  for  negligible  recombination  in 
the  p+n+  space-charge  region,  defined  by  the  recombination 
in  the  quasi-neutral  p+-emitter  [McD85] : 

Jn (°)  — ~i)T-  pm  (2.5) 


where  ./,v o is  the  electron  saturation  current  density  for  the 
p+-emitter.  Combining  (2.4)  and  (2.5)  then  yields  the  total 
device  current  density  in  terms  of  p( 0)  and  p ^H'+ \ ; 


.7 


l = 


•7,V0  qDpf\ 

”i/Xuf  Wf  ) 


(2.6) 


A second  expression  for  J A in  terms  of  the  boundary  hole 
densities  is  derived  as  follows.  Applying  JA  = JP  + JXf  with 

p-n'  at  the  ed9e  (;r  = Wj)  of  the  quasi-neutral  base,  and 
noting  that  JN  ~ ./jV  (0),  we  obtain  [McD85] 

JN  (0)  ~ —JA  + (jDa~  (,•  = TT  - ) (2.7) 

where  b is  the  electron-hole  mobility  ratio  and  D A is  the 
ambipolar  diffusivity  in  the  base  region.  Now  (2.1),  (2.5), 

and  (2.7)  give 
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J A = 


l + b JX0 

b "i/Nof 


p(  0)  + 


1 + b coth  ( W/  LA ) 


■qD 


- A 


(”’/)  • 


(2.8) 


The  three  equations  (2.3),  (2.6),  and  (2.8)  define  the 

boundary  hole  densities  in  terms  of  JA  but  also  involve 
TT',  another  unknown  variable.  To  define  I,\(VA)  then,  two 
additional  expressions  relating  VA  to  the  unknown  variable 
must  be  derived.  Referring  to  Fig.  2.1,  we  write  V\  as  the 
sum  of  the  internal  voltage  drops : 

^A  = y'j  ” + I . " " + Vf„,  + IaHs  + IbRmos  (2.9) 


where  IA  = AJAl  Iq  — AJb,  and  A is  the  device  area.  In  (2.9), 


+ ..+ 


v; " = VTln 


P(  0)NDf 


(2.10) 


is  the  drop  across  the  emitter-buffer  junction, 

(">') 


= -VTlv 


A 


Ue 


(2.11) 


is  the  drop  across  the  buffer-base  junction,  XDe  is  the  epi 
doping  density,  and 


JAf1r 

+ b)nvp{x) 


- It 7 — -In 

6 + 1 


p(Trw) 


(2.12) 


is  the  drop  across  the  (conductivity-modulated)  n_-epi  region, 
which  is  approximated  as  the  integral  of  the  electric  field 
[McD85  ] from  IT  ^ to  IT  the  point  at  which  />(;»•)  in  (2.1)  equals 
Dr  • The  last  two  terms  in  (2.9)  are  ohmic  drops  across  the 
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extrinsic  series  resistance  /?<?  and  across  the  DMO ST  through 
which  the  BJT  base  current,  /#,  flows  in  the  steady  state. 
Note  that  the  DMO  ST  drop  equals  (for  normal  operation  below 
latch-up)  the  reverse  bias  on  the  base-collector  junction, 
which  we  have  assumed  is  small.  The  final  relation  needed 
to  characterize  Ia(Va)  is 


which  equates  the  base  current  to  the  integrated  recombination 
in  the  quasi-neutral  base  and  emitter  regions.  In  (2.13),  r// 
is  the  high-injection  carrier  lifetime  in  the  base. 

The  resulting  system  of  equations  underlying  (2 . 9)  can  be 
solved  numerically  to  give  IA  ( V'a),  and  incidently  the  BJT 
current  gain,  /?  = (IA  - Ib)  /Ib/  which  indicates  (for  the  one- 
dimensional structure)  how  the  VIGBT  (BJT  emitter)  current 
is  divided  between  the  DMO  ST  and  the  BJT  collector.  The 
results  depend  on  the  structural  device  parameters  and  doping 
densities,  including  the  buffer-layer  width  Wf  and  (uniform) 
doping  density  NUf,  as  well  as  r the  high-injection 
carrier  lifetime  in  the  n“-base  region,  and  JY 0,  the  electron 
saturation  current  density  in  the  p+^emitter  region.  To  make 
the  one-dimensional  analysis  representative  of  actual  VIGBTs, 
we  use  a semi-empirical  constant  K A «1)  that  accounts  for 
multidimensional  carrier  flow.  Crudely  I{A  is  the  ratio  of 
the  collector  and  emitter 


HT+ir 


(2.13) 


li 


areas  (Fig.  2.1),  and  therefore 
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/* -♦  A.4/J/[l + (1  - AVj)/J].  we  stress  that  models  based  strictly  on 
one-dimensional  analyses  [Nak85,  Kuo86,  Hef86]  can  be  quite 
misleading. 


2.2.2  Transient  Turn-Off 

The  unique  transient  turn-off  characteristic  of  the  VIGBT, 

comprising  an  initial  rapid  drop  in  forward  current  followed 

by  a slower  decay,  was  previously  modeled  [Fos86]  without 

accounting  for  the  buffer  layer.  An  extended  version  of 

this  model,  derived  straightforwardly  following  Fossum  and 

McDonald  [Fos86],  is  used  here  to  account  for  the  buffer- 

layer  effects.  The  initial  rapid  drop,  A/,  is  modeled  as  in 

Fossum  and  McDonald  [Fos86],  but  is  based  on  the  steady-state 

characterization  including  the  buffer  layer  described  above. 

The  slower  transient,  /,(,),  is  modeled  as  Qp(t) /r,p  [Fos86]  , 

where  rlp  is  the  (average)  base  transit  time,  but  now  the 

hole  charge  in  the  quasi-neutral  base,  Qpr  is  governed  by 

a modified  [Hef86]  differential  equation  (in  time  t)  which 

accounts  for  the  buffer  layer  in  the  relation  between  Qp 

and  Qnt  the  minority  electron  charge  in  the  quasi-neutral 
emitter : 


Ir 


Qp(*)  , 

Tll 


'IQp  ( * ) , Qn(t)  <I(J„  (t) 

'It  + r»  +— — ="• 


(2.14) 


Where  r„  is  the  effective  recombination  time  for  the  electrons 
in  the  emitter  and  rfg./rf,  is  negligibly  small.  Using 
the  steady-state  analysis  in  the  previous  subsection,  in 
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Par^4cular  (2.3),  (2.6),  and  (2.8),  and  assuming  quasi-static 
conditions  for  the  transient,  we  can  show  that 

Qn  — C i Q~  + C-tQp 


(2.15) 


where  C\  and  C'>  are  constants  given  by 

4rH./,Y0  1 


Ci  = 


Aq1  ( IT'/1  - .r nj  1 - TT/./yVoiV/J f/lxiDj, f U J 


(2.16) 


and 


C>  = 


4H  / t„  of 


(2.17) 


<i  ( H'b  - )-  nr  1 - WfJ d // ItqDpfnj 

In  (2.16)  and  (2.17),  rdm  is  the  width  (r d in  Fig.  2.1)  of 
the  base-collector  junction  depletion  region  corresponding 
to  IA  = /,  [Fos86]  . 

The  differential  equation,  which  is  quasi-static,  is 
nonlinear,  but  can  be  solved  analytically  to  give 


I.\  (*)  - 


1 + i1  - ™i>{-t/TH((ff))\ 


(2.18) 


where  /,  =/„-A/  is  the  current  at  the  beginning  of  the  slower 
transient  phase  of  the  turn-off  [/n  is  characterized  by  the 
steady-ON-state  analysis  underlying  (2.9)],  and  Tll{(ff)  is  an 
effective  carrier  lifetime  which  differs  from  ru  because  of 
the  presence  of  the  buffer  layer: 

T f-^  Df 

(2.19) 


J: ~ + 4./,V"II  fXpf 


Tmeff)  -H  q ( ivl3  - Xdn, 


2 ■) 
nr 


The  remaining  terminology  in  (2.18)  and  (2.19)  is  defined 
in  Fossum  and  McDonald  [Fos86]  . We  note  that  if  H>  approaches 
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zero  (viz.,  no  buffer  layer),  TIf(eff)  approaches  rn , and  (2.18) 
simplifies  to  the  result  derived  in  Fossum  and  McDonald 
[Fos86]  . Physically,  the  second  term  on  the  right  side 
of  (2.19)  represents  the  effect  of  recombination  in  the 
p —emitter  on  IA(f)r  influenced  by  the  buffer  layer  through 
which  the  recombining  electrons  must  flow.  Hence  even  with  a 
buffer  layer  in  place,  if  ./Vo  becomes  negligibly  small,  then 
aiso  TI1{(ff)  ~ t,{  and  I t (/)  ~ I\txp{ -t/rH  ) . 

2 • 3 Theoretical-Experimental  Interpretations 

VIGBTs  fabricated  at  RCA/DSRC  were  measured.  Typical  IA(V\) 
data  taken  from  an  n-channel  device  (with  the  gate  biased 
at  10  V,  about  8 V above  the  threshold  voltage)  are  plotted 
in  Fig.  2.3.  The  simulated  characteristic,  also  shown,  was 
derived  from  the  model  outlined  in  Section  2.2,  using  known 
structural  device  parameters  and  doping  densities.  The  key 
model  parameters,  i.e.,  rH.  Jv„,  /?.v,  and  B]IOSr  were  extracted 

from  additional  steady-state  and  transient  measurements  as 
follows . 

The  transient  turn-off  characteristic,  modeled  by  (2.18), 
was  measured  for  an  ON-state  current  /„  low  enough  to  ensure 
that  IA(f)  is  exponential  as  described  by  the  numerator  of 
(2.18).  This  measurement  gives  r//(f//)  directly,  subject  of 
course  to  the  validity  of  the  model  assumptions.  (We  have 
found  that  if  the  OFF-state  voltage  is  high  enough  (>  200  V)  to 
cause  a significant  transient  variation  of  the  quasi-neutral 
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Fig.  2.3 


Measured  and  calculated  steady-ON-state  current 

V^GB?1'  <t ^"-°12t5a°geo  Cht5^ct-eristics  for  n-channel 
wst  (T  - 25  C)  . Model  parameter  values  used 

are  rf{  12  // s,  JXo  = l.i  x 10  13  A/cm2,  i?s  = o 08 

, and  Ituos  = 0.04  (l,  also,  X,„  - 5 x 10“  cm"3, 

"/a  "o'b^'  = 1 * 10,4  cm~3'  - « m. 
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base  width  (IT  n-  r,i),  then  rH{t ff)  can  be  over-estimated.)  Using 
(2.19)  then,  th  and  JN()  were  evaluated  by  fitting  the  model 
prediction  for  the  turn-off  time  T0  and  its  dependence  on  /„. 
T0  is  the  time  required  for  IA(t)  to  drop  to  0.1/„;  it  was 
measured  by  pulsing  the  gate  ON  (10  V)  and  OFF  (0  V),  while 
the  anode  was  biased  through  a resistive  load.  The  model- 
predicted  T0  vs.  was  derived  by  first  modeling  the  steady 
state,  and  using  the  results  to  characterize  the  transient 
turn-off  as  described  in  Section  2.2. 

Plotted  in  Fig.  2.4  are  the  measured  T0  vs.  data 

along  with  the  model  characteristic  (which  corresponds  to 
th  ~ 12  [is  and  ./jVo  ~ 1.1  x 10-13  A/ cm2  at  T = 25°C)  . The 
theoretical-experimental  fits  illustrated  in  Figs.  2.3  and 
2.4  also  involve  iterative  evaluation  of  Rs  (~  o.08  ft) 

and  Ruos  (~  0.04  ft)  , the  dependence  of  which  on  gate 

voltage  enables  separation  from  the  series  resistance.  Since 
the  T0(l o)  characteristic  is  not  strongly  dependent  on  the 
small  resistances  Rs  and  RMOSr  their  extraction  is  made 
mainly  from  the  IA(VA)  characteristic.  in  all  cases  the 
model  predictions  matched  the  measurements  well,  without  much 
empiricism  involved.  The  theoretical -experimental  agreement 
gives  credence  to  our  model  and  demonstrates  the  utility  of 
the  parameter-extraction  techniques. 


We  now  show  how  the  model  can  be  used  to  give  physical 
insight  and  thereby  aid  the  optimal  design  of  the  buffer 
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2.4 


Measured  and  calculated  turn-off  time  vs.  ON-state 

~ 25°  C)  ’ ThS  model  Parameter  values 
used  are  the  same  as  those  given  in  Fig.  2.3 


Fig. 
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Fig.  2.5  Measured  forward-voltage  drop  vs.  turn-off  time 
taken  from  p-channel  VIGBTs  having  different 
llfetimes  in  the  epi  base  region.  The 
ON-state  current  was  held  constant  as  indicated 
°f.  a pe}.le^  araa  of  120  X 120  mil2;  the  OFF-state 
oltage  ] B - 400  V.  The  data  points  are  connected 

dependences*!  t0  emPhasi^  th*  Parametric 
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layer.  Plotted  in  Fig.  2.5  are  additional  data  taken  from 
p-channel  VIGBTs  fabricated  at  RCA/DSRC.  The  measured  0N- 
state  voltage,  VAr  versus  T0  for  three  buffer-layer  doping 
densities,  j\Af,  was  obtained  from  devices  in  which  the 
carrier  lifetime  was  varied  by  heavy-metal  diffusion.  The 
ON-state  current  I0  was  held  constant.  The  data  show  that 
predominantly  for  a given  Ta , ]A  increases  as  NAf  is  increased, 
thus  suggesting  that  NAf  be  low  for  optimal  design.  This 
suggestion  is  contrary  to  that  emanating  from  corresponding 
(n-channel  VIGBT)  data  [Nak85],  which  show  that  for  a given 
T°'  } [ decreases  as  XDf  is  increased,  and  from  theoretical 
predictions  [Kuo86] . This  apparent  contradiction  can  be 
explained  with  reference  to  our  analysis.  We  note  that 
for  a given  th,  an  increase  in  NUf  will  tend  to  decrease 
the  carrier  injection  level  in  the  n'-base  region,  which 
concomitantly  tends  to  decrease  T0,  but  to  increase  VA  because 
of  the  higher  Vfpi  in  (2.9).  if  the  decreased-^  tendency  is 
predominant,  and  the  increase  in  VA  ±8  less  significant 
because,  for  example,  Vept  is  negligibly  small,  then  the 
implication  regarding  optimal  XDf  [Nak85,  Kuo86]  is  obtained. 
Otherwise,  the  opposite  implication  in  Fig.  2.5  is  obtained. 
Based  on  our  studies,  we  deduce  that  this  opposite  effect  is 
evident  only  when  the  injection  level  in  the  n~-base  region 
is  predominantly  low,  thus  yielding  a very  high  i;,,,  as  in 
the  ordinary  power  MOSFET.  In  fact  the  VIGBT  with  a very 
high  NDf  (and/or  wide  TT'y)  behaves 


as  a power  MOSFET  with 
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no  benefit  from  conductivity  modulation.  Thus  optimizing  the 
buffer-layer  design  is  more  involved  than  previously  stated 
[Nak85,  Kuo86,  Hef86]  and  must  involve  a more  comprehensive 
study  of  the  effects  of  all  regions  of  the  VIGBT . 

Such  a (theoretical)  study  is  illustrated  in  Figs.  2.6,  2.7, 
and  2 . 8 where  we  plot  calculated  VA  - Tn  curves  for  different 
values  of  key  device  parameters,  Nuf.  Wj,  and  ./jVo . The  curves 
were  derived  by  varying  rH  between  0.1  and  100  /is  while 
holding  J0  and  Vb  constant.  Note  in  each  figure  that  points 
corresponding  to  a specific  th  define  a locus  on  the  VA  - T0 
coordinate  system  having  a negative  slope;  viz.,  shortening 
T0  via  the  parametric  changes  indicated  necessarily  increases 
1 4 . In  Fig.  2.6,  the  insight  discussed  above  regarding  the 
influence  of  Nof  is  illustrated.  The  convergence  of  the 
dif f erent  — curves  varies  as  the  key  device  parameters, 
especially  JNOf  are  varied.  Additional  calculations  show  that 
increasing  J ^ 0 moves  the  point  of  convergence  to  the  left, 
thereby  enhancing  the  tendency  shown  in  Nakagawa  and  Ohashi 
[Nak85]  and  predicted  in  Kuo  and  Hu  [Kuo86] . Decreasing 
•7;vo  moves  the  point  to  the  right,  enhancing  the  tendency 
reflected  in  Fig.  2.5.  The  convergence  of  the  curves  in  Fig. 
2.6,  although  predicted  by  our  analysis  based  on  high-level 
injection  in  the  n —base  region,  nevertheless  implies  the 
onset  of  low  injection  which,  as  mentioned  above,  intensifies 
the  opposite  \A  — T0  vs.  Nd/  tendency  shown  in  Fig.  2.5. 


25 


Fig.  2.6 


Simulated  forward-voltage  drop  vs.  turn-of- 
for  two  values  of  Nnfr  with  ?Vf  m °* 

X 10  A/cm  » an°i  U'B  = 45  //m. 


time 
= 1 
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Fxg.  2.7  Simulated  forward-voltage  drop  vs 
for  two  values  of  JXo,  with  XDf 
Uf  ~ 10  /'"W  and  WB  = 45  /,m. 


turn-off  time 
2 x 1017  cm-3, 


27 


Fig. 


2.8  Simulated  forward-voltage  drop  vs. 
for  two  values  of  \Yf,  with  XDf  = 

and  for  no  buffer  layer,  all  with 
A/cm  , and  II  q = 45  /im . 


turn-off  time 
2 x 1017  cm’3, 
■/.Vo  = 5 x 10-13 
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The  significance  of  J ^ u is  further  emphasized  in  Fig.  2.7. 
For  a specified  buffer  layer,  the  calculated  data  plotted 
show  that  increasing  ./,Vn  decreases  T„  but  also  increases  V'i, 
for  the  same  /() . This  is  similar  to  the  effect  of  the  buffer 
layer:  the  BJT  f3  is  lowered,  thereby  reducing  the  quasi-Fermi 
level  separation  (forward  bias)  at  the  emitter-base  junction 
and  the  injection  level  in  the  n“— base  region. 

Finally  in  Fig.  2.8,  calculated  dependences  on  Wf  are 
shown.  Increasing  IT/,  analogous  to  increasing  ./v„,  decreases 
T0  but  increases  VA.  For  comparison,  the  case  \Vf  = 0 (no 
buffer  layer)  is  also  shown,  indicating  that  adding  the  buffer 
produces  the  same  tendencies  as  increasing  ./iVo . However,  as 
can  be  inferred  from  our  discussion,  the  effects  of  the  buffer 
layer  on  the  VIGBT  performance  cannot  be  ascertained  without 
considering  the  significance  of  all  key  device  parameters . 
We  note  for  example  that  the  benefit  of  the  buffer  layer  is 
enhanced  by  a reduction  of  TT'#  because  V,p,  in  (2.9)  tends 
to  become  less  significant,  and  hence  the  T„  - VA  trade-off 
discussed  is  less  critical. 

2 . 4 Summary 

A simplified  physical  model  for  the  VIGBT  has  been  used 
to  interpret  both  steady-ON-state  and  transient  turn-off 
measurements  reflecting  significant  effects  of  the  buffer 
layer.  The  theoretical-experimental  study  has  (i)  provided 
good  physical  insight,  (ii)  implied  optimal  buffer-layer 
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design  criteria,  (iii)  described  model  parameter-extraction 
techniques,  and  (iv)  indicated  shortcomings  in  corresponding 
results  of  previous  work  due  to  a lack  of  comprehensiveness 
in  the  analyses . 


CHAPTER  3 

STATIC  AND  DYNAMIC  LATCH-UP  IN  THE  LIGBT 


3 . 1 Introduction 

This  chapter  presents  static  and  dynamic  latch-up  models 
for  the  lateral  IGBT  (LIGBT)  [Darw84,  Pat86] . in  Chapter 
2,  an  analytic  buffer-layer  model  for  the  VIGBT,  which  is 
similar  in  structure  to  the  LIGBT  except  the  anode  layout, 
was  developed  and  experimentally  verified.  The  physical 
insight  gained  from  the  modeling  of  the  VIGBT,  and  from 
PISCES  simulations  discussed  here  is  used  to  develop  the 
latch-up  models  for  the  LIGBT. 

The  LIGBT  is  a high-voltage  power- switching  device  that 
can  be  used  in  high-voltage  integrated  circuits  (HVICs) . 
However  the  LIGBT,  and  VIGBT  as  well  can  be  latched  into  an 
undesirable  state  in  which  gate  control  is  lost.  This  state, 
which  can  be  triggered  either  statically  (by  sufficiently 
high  quasi-static  current)  or  dynamically  (by  sufficiently 
high  displacement,  or  transient  charging  current),  results 
from  the  latch-up  [Rob86]  of  the  parasitic  pnpn  structure 
in  the  LIGBT.  The  LIGBT  must  therefore  be  designed  to  avoid 
latch-up  in  the  operation  of  the  HVIC.  Computer-aided  design 
that  includes  latch-up  modeling  for  both  device  and  circuit 
simulation  is  required. 


30 


31 


Very  little  previous  work  on  the  simulation  of  latch- 
up  (especially  dynamic  latch-up,  which  is  the  predominant 
problem  in  HVICs)  in  LIGBTs  has  been  done  [Pat86,  Darw87, 
Chow87] . In  this  chapter,  we  describe  an  insightful  study 
of  the  static  and  dynamic  latch-up  phenomena  in  the  LIGBT, 
which  culminates  in  the  development  of  a physical  model  for 
SPICE  simulation  of  HVICs.  The  study  is  based  on  extensive 
simulations  of  latch-up  in  the  LIGBT  using  the  two-dimensional 
device  simulator  PISCES  [Pin84] . The  simulations  reveal  the 
internal  mechanisms  underlying  static  and  dynamic  latch— up 
£ot  a variety  of  excitations,  and  they  show  that  the  lateral 
pnpn  structure  becomes  an  effective  pin  diode  after  latch-up 
because  of  the  excessive  conductivity  modulation  near  the 
surface.  In  fact,  the  conductivity  modulation  is  observed 
to  be  the  governing  mechanism  in  the  latching  and  in  the 
latched  state  (in  accord  with  analyses  of  latch-up  in  CMOS 
ICs  [Sei87,  Cou88]),  and  hence  must  be  properly  accounted 

for  in  the  simulations.  Also  noteworthy  is  the  discovered 
significance  of  non-quasi-static  bipolar  behavior  [Foss86] 
in  the  dynamic  latch-up,  as  well  as  the  possibility  of 
different  controlling  charging  currents,  depending  on  the 
transient  terminal  conditions . Dynamic  latch-up  can  occur 

for  very  fast  gate-voltage  switching  and/or  for  very  high 
blocking  voltages. 

The  SPICE  model  for  the  LIGBT,  including  latch-up,  is 
developed  by  extending  the  charge-based  VIGBT  model  in  Chapter 
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2,  which  comprises  one-dimensional  physical  representations 
of  the  ambipolar  carrier  transport  in  different  regions  of 
the  device.  The  model  is  implemented  in  SPICE  via  user-defined 
controlled  sources  (FORTRAN  subroutines)  and  is  verified 
by  measurements  of  test  devices  fabricated  at  AT&T  Bell 
Laboratories , as  well  as  by  the  PISCES  simulations. 

3.2  Numerical  Simulations 

To  gain  the  physical  insight  needed  to  develop  the  latch- 
up  model,  PISCES  [Pin84]  was  used  to  study  both  static  and 
dynamic  latch-up  in  the  representative  LIGBT  structure  shown 
in  Fig.  3.1.  PISCES  is  a two-dimensional  device  simulator,  for 
dc  or  transient  terminal  conditions,  that  utilizes  finite- 
element  approximations  and  a Gummel/Newton  numerical  method 
to  solve  the  discretized  semiconductor  equations  defining 
the  carrier  transport  [Cou88],  including  contemporary  models 
for  the  semiconductor  parameters.  Doping-density  profiles 
and  structural  dimensions  corresponding  to  our  test  devices 
from  AT&T  Bell  Laboratories  (see  Section  3.3)  were  assumed. 
Although  these  LIGBTs  are  dielectrically  isolated  [Luc88],  the 
insight  afforded  by  the  simulations  and  modeling  we  describe 
is  helpful  in  analyzing  the  junction-isolated  resurf  device 
structure,  which  includes  an  additional  vertical  constituent 
bipolar  transistor  and  which  could  include  a shorted  anode 
[Darw87,  Chow87,  Fos88]  . As  indicated  in  Fig.  3.1,  in 
the  steady  on-state,  the  constituent 


DMO ST  gate  is  biased 
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V8>Vt 


Fig.  3.1  LIGBT  test  structure  with  external  circuitry  and 
gate  drive  for  a turn-off  transient. 
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above  the  threshold  voltage  and  the  anode  is  biased  above 
the  source/cathode  (ground)  through  a load  resistor.  In  the 
transient  turn-off,  the  gate  voltage  drops,  in  a finite  fall 
time,  to  zero. 

3.2.1  Static  Latch-Up 

The  static  latch-up  is  demonstrated  in  Fig.  3.2  where 
a PISCES-simulated  dc  current-voltage  characteristic  of  the 
LIGBT  is  shown.  To  enable  simulations  of  the  entire  (multi- 
valued current)  characteristic,  in  reasonable  computation 
time  [Cou88] , the  applied  anode  voltage  was  increased 
incrementally  to  derive  the  sub-latch  current,  but  after 
latch-up  (onset  at  point  where  dVA/dIA  = 0)  , the  anode  was 
driven  with  an  increasing  current  and  the  voltage  was  derived 
from  the  simulation.  Above  the  onset  of  static  latch-up,  a 
negative-resistance  region  reflects  the  regenerative  latching 
mechanism  involving  the  lateral  pnp  transistor  (the  anode  is 
the  emitter)  and  the  parasitic  npn  transistor  at  the  cathode 
(the  DMO ST  source  is  the  emitter);  see  Fig.  3.1. 

The  regeneration  occurs  because  holes  injected  from  the 
forward-biased  anode  are  collected  and  flow  laterally  though 
the  npn  base  pinch  resistor  RpBr  producing  a forward  bias 
(predominantly  under  the  DMO ST  gate)  on  the  npn  emitter- 
base  junction.  Electrons  thereby  injected  (predominantly 
laterally)  from  the  source  drive  the  npn  base  harder, 
producing  additional  forward  bias  on  the  anode  junction  and 
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LA.  — VA  (PISCES) 

Fig.  3.2  PISCES-simulated  LIGBT  current -voltage  character- 
istic  showing  static  latch-up.  The  gate  bias  is 
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additional  hole  injection.  The  process  thus  enhances  the 
conductivity  modulation  of  the  pnp  base  region,  including 
the  adjacent  space-charge  region  at  the  npn  base-collector 
junction.  This  junction  is  ultimately  neutralized  as  I\ 
increases  by  the  hole-electron  plasma,  and  the  steady  latched 
state  is  characterized  by  the  effective  pin  diode  that 
materializes  [Sei87,  Cou88] . 

For  a point  in  Fig.  3.2  below  latch-up,  Fig.  3.3  shows 
the  simulated  current  flowlines  within  the  device.  Normal 
LIGBT  operation  is  evident;  the  channel  current  of  the  DMO ST 
supplies  electrons  to  the  base  of  the  lateral  pnp  transistor 
which  support,  through  recombination,  the  hole  transport 
current.  The  LIGBT  (anode)  current  equals  the  channel  current 
plus  the  pnp  collector  current . Note  that  most  of  the  hole 
current  flows  laterally  through  RpH  in  the  npn  base  region, 
but  the  induced  ohmic  drop  is  not  sufficient  to  significantly 
forward-bias  the  source  junction;  the  electron  current  is 
confined  predominantly  to  the  DMO ST  channel . 

The  simulated  current  flowlines  in  Fig.  3.4  reflect 
operation  above  the  onset  of  latch-up  in  Fig.  3.2.  In 
contrast  to  the  condition  pictured  in  Fig.  3.3,  we  see  in 
Fig.  3.4  significant  injection  of  electrons  from  the  n+-source 
(emitter)  into  the  npn  base  and  collector  (pnp  base)  . The  hole 
current  flowing  through  Rrn  produces  a voltage  drop  which  in 
this  case  is  high  enough  to  activate  the  npn  transistor. 
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Fig.  3.3 


PISCES -simulated  current  flowlines  in  the  LIGBT  in 
. *91  q * 1 ^ased  below  latch-up.  The  anode  current 
mu  X,  .10  A//,m  and  the  gate  voltage  is  10  V. 

arS  plotted  such  that  equal  currents 
, r between  any  two  adjacent  flowlines.  For 
the  whole  device  (a),  A/  = 0.2  x 10~5  A//,m;  in  the 
vicinity  of  the  DMO ST  (b)  , A/  = 0.1  x 10"5  A//,m 
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Fig.  3.4  PISCES -simulated  current  flowlines  in  the  LIGBT  in 
Fig.  3.1  biased  above  latch-up.  The  anode  current 
is  3.4  x 10  A///m  and  the  gate  voltage  is  10  V.  For 

the  whole  device  (a).  A/  = 0.2  x 10“5  A//,m;  In  the 
vicinity  of  the  DMO ST  (b)  , A/  = 0 . 1 x 10~5  A//,m. 
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and  initiate  the  regenerative  process  in  the  pnpn  structure. 
The  predominant  electron  injection  occurs  laterally  from  the 
source  under  the  DMOST  gate  where  the  forward  bias  is  highest. 
For  sufficiently  high  current,  the  base-collector  junction  of 
the  npn  and  pnp  transistors  is  highly  conductivity-modulated 
(as  in  a quasi-saturated  transistor  [Sze81])  by  holes  and 
electrons  just  below  the  surface,  and  an  effective  pin  diode 
exists  between  the  anode  and  the  cathode/source  of  the  LIGBT . 
The  current-voltage  characteristic  of  this  pin  diode  defines 
the  negative  resistance  in  Fig.  3.2,  as  well  as  the  conductance 
of  the  latched  device  [Sei87,  Cou88] . In  this  condition, 
removal  of  the  gate  voltage  does  not  necessarily  result  in 
the  turn-off  of  the  LIGBT. 

As  a lead-in  to  the  next  subsection  on  dynamic  latch- 
up,  we  note  here  that  the  simulated  steady-state  Ii{V.\) 
characteristic  in  Fig.  3.2  is  not  the  characteristic  which 
would  be  observed  in  common  current -voltage  measurements.  In 
a curve-tracer  measurement,  the  voltage  is  swept  up  and  down, 
and  hence  a multi-valued  current  characteristic  like  that  in 
Fig.  3.2  would  show  an  apparent  hysteresis.  This  hysteresis 
could  be  eliminated  by  sweeping  the  current  instead  of  the 
voltage,  but  the  dynamic  nature  of  these  measurements  also  can 
cause  a hysteresis.  In  fact,  the  hysteresis  observed  in  curve- 
tracer  measurements  of  latch-up  can  be  attributed,  in  part, 
to  the  dynamic  nature  of  the  current.  PISCES  and  SPICE  (using 
the  LIGBT  model  presented  in  Section  3.3)  simulations  of  a 
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sweeping-current  measurement  show  the  hysteresis,  with  the 
latch-up  occurring  at  a higher  current  when  I\  is  increasing 
with  time  than  when  it  is  decreasing  with  time.  The  actual 
steady-state  characteristic  lies  within  the  hysteresis  loop. 
These  results  are  consistent  with  HP-4145  current -driven 
measurements  we  have  made  and  can  be  understood  based  on 
the  insight  gained  in  the  next  subsection.  Another  reason 
for  the  discrepancy  between  a (two-dimensionally)  simulated 
steady-state  latching  characteristic  and  a measurement  is  the 
(three-dimensional)  inhomogeneities  in  the  actual  device  that 
can  cause  nonunif orm  latching  and  hence  also  a hysteresis . 

3.2.2  Dynamic  Latch-Up 

With  regard  to  actual  HVIC  operation,  dynamic  latch-up  of 
the  LIGBT  is  the  real  problem.  When  the  gate  voltage  is 
dropped  to  turn  off  the  LIGBT  from  an  unlatched  on-state, 
the  device  can  latch-up  and  continue  to  conduct.  There  are 
two  distinct  conditions  which  induce  dynamic  latch-up.  One  is 
when  the  gate  voltage  drops  very  fast  and  induces  excessive 
displacement  current  through  the  gate  oxide  that  drives  7?,,^. 
The  other  one  is  when  the  (off-state)  blocking  voltage  is 
quite  high  and  induces  excessive  charging  currents  within 
the  LIGBT  during  the  switching  transient  that  drive  /?,,#. 
Of  course  the  probability  of  dynamic  latch-up  increases  as 
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the  steady  on- state  condition  approaches  the  static  latch-up 
onset  (depicted  in  Fig.  3.2). 

For  a resistive  load  as  shown  in  Fig.  3.1,  Fig.  3.5 
shows  PISCES-simulated  transient  turn-off  characteristics, 
I.\(t),  for  a short  gate-voltage  fall  time  (t f = 0.5  ns)  and 
for  two  different  initial  steady  on-state  currents  (Iu/Il  = 
0.25,  0.86,  where  Ii  is  the  static  latch-up  onset  current). 
The  off -state  voltage  used  is  low  (Vjj  = 10  V)  . For  the 
lower  initial  current  It),  the  turn-off  is  normal,  but  for  the 
higher  J0,  a delay  is  seen  Fig.  3.5.  The  device  current 
flowlines  for  the  latter  case  at  three  different  times, 
shown  in  Fig.  3.6,  reveal  the  nature  of  the  transient  delay. 
The  induced  gate-oxide  displacement  current  (predominantly 
electrons  flowing  off  the  gate-drain  overlap  capacitance  into 
the  pnp  base)  tends  to  drive  the  regenerative  process  in  the 
pnpn  structure.  Although  initially  {t  = 0.44  ns)  the  gate 
depletion  (dis)  charging  current  diverts  holes  away  from  /?/)0, 
later  (/  = 2.55  ns),  long  after  the  gate  voltage  has  been 
removed,  the  transient  voltage  drop  across  Rl>n  is  sufficient 
to  activate  the  npn  transistor,  and  significant  conductivity 
modulation  prevails  in  the  device.  The  reason  for  this 
delayed  activation  is  that  the  discharging  of  the  depletion 
capacitance  occurs  relatively  fast,  and  precedes  a transient 
(electron  current)  drive  on  the  pnp  base  from  the  discharging 
overlap  capacitance,  which  increases  the  hole  current  through 
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IA/IO  - t (PISCES) 


Fig.  3.5 


PISCES -simulated  transient  turn-off  characteris- 
tics  for  a short  gate-voltage  fall  time  of  0.5  ns 
and  for  two  different  initial  on-state  currents, 
h/Ir  where  //.  = 2.2  x 10~5  A//,m  is  the  static 
latch-up  onset  current . 
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Fig.  3.6  Current  flowlines  (A/  = 0.1  x 10-5  A///m)  in  the 
vicinity  of  the  DMO ST  at  three  different  times  in 
the  transient  turn-off  simulation  in  Fig.  3.5  for 

/it///,  = 0.86. 
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RpB  • However,  in  this  simulation,  this  transient  (quasi- 
dynamic)  latch-up  is  not  sustained  because  of  the  low  VB 
and  the  simple  circuit  configuration  which  are  insufficient 
to  force  a steady  latched  condition.  The  LIGBT  ultimately 
turns  off  after  the  delay  caused  by  temporary  activation  of 
the  effective  pin  diode.  In  real  HVICs,  this  quasi-dynamic 
latch-up  can  be  converted  to  an  actual  dynamic  latch-up  by 
stray  capacitance  and/or  inductance  and  by  higher  VB . 

PISCES-simulated  transient  turn-off  characteristics  for 
higher  (50  V,  400  V)  and  longer  /y  (50  ns)  are  shown  in 
Fig.  3.7.  In  the  case  of  the  higher  VB,  dynamic  latch-up 
is  evident,  even  though  the  gate  displacement  current  is 
insignificant.  As  indicated  by  the  current  flowlines  in  Fig. 
3.8,  the  latch-up  in  this  case  results  because  of  substantial 
charging  currents  induced  in  the  LIGBT  which  flow  through  RpB. 
Note  in  Fig.  3.8  that  the  lateral  current  flowing  through  Rj)B 
increases  with  time.  This  is  explained  as  follows.  During  the 
first  phase  of  the  turn-off  transient,  IA(t)  drops  rapidly  due 
to  the  removal  of  the  DMO ST  channel  current  [Fos86] . Then,  in 
the  resistive-load  circuit,  the  anode  voltage  VA  increases 
accordingly,  which  causes  three  distinct  internal  charging 
currents  to  flow,  all  of  which  are  manifested  by  holes  flowing 
through  RpB  which  tend  to  induce  the  latch-up.  One  is  the 
displacement  current  that  charges  the  base-collector  junction 
as  its  reverse  bias  increases.  This  current  is  defined  by 
the  transient  depletion  capacitance  of  the  junction  [Fos86] . 
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Current  flowlines  (AI  = 0.1  x 10~5  A/,/m)  in 
vicinity  of  the  DMO ST  at  two  different  times 
the  transient  turn-off  simulation  in  Ficr  3 7 
VB  = 400  V. 
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Fig.  3.8 
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Another  current  is  the  transient  base  transport  (collector) 
current  of  the  pnp  transistor,  which  increases  because  of  the 
shrinking  quasi-neutral  base  width  [Fos86] . The  third  current, 
which  has  previously  not  been  recognized  as  being  significant, 
is  the  component  of  the  quasi -neutral  base  discharging  current 
that  flows  out  the  collector.  Conventional  quasi-static 
transistor  models  [Sze81]  do  not  account  for  this  charge 
partitioning  [Foss86] , which  actually  reflects  a non-quasi- 
static behavior,  and  simply  associate  all  the  base  charging 
current  with  the  emitter.  For  fast  transients,  like  latch-up 
in  the  LIGBT,  the  base  charge  partitioning  must  be  used  to 
account  for  transient  charge  redistribution  [Foss86] . Our 
SPICE  simulations  (discussed  in  the  next  section)  reveal 
that  the  non-quasi-static  collector  component  of  the  base 
discharging  current  is  critical  in  simulating  dynamic  latch- 
up  in  the  LIGBT.  In  fact,  this  component  is  typically  (for 
our  devices)  the  largest  discharging  current  flowing  during 
the  turn-off  transient;  it  controls  the  dynamic  latch-up. 

3.3  SPICE  Model 

Using  the  physical  insight  gained  from  the  PISCES 
simulations  described  in  Section  3.2,  we  extend  the  physical 
SPICE  model  for  the  LIGBT  [Fos88]  to  properly  account  for 
latch-up.  The  network  representation  of  the  new  SPICE  model 
is  shown  in  Fig.  3.9,  superimposed  on  the  LIGBT  structure.  All 
the  elements  in  the  model,  the  underlying  physics  of  which  is 
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Fig.  3.9  Network  representation  of  the  SPICE  LIGBT  model 
superimposed  on  the  device  structure . 
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described  below  or  in  Fossum  et  al . [Fos88] , are  implemented 
in  SPICE  via  user-defined  controlled  sources  (UDCSs)  [Vee86] , 
which  are  FORTRAN  subroutines  in  SLICE  [Har84],  an  enhanced 
version  of  SPICE2.  The  new  model  in  SPICE  provides  a mixed- 
mode device/circuit  simulation  capability,  which  in  fact  is 
essential  for  realistic  simulations  of  latch-up  in  HVICs 
and  which  is  prohibitive  in  PISCES  because  of  limitations 
on  the  allowed  terminal  drive  conditions  and  because  of 
the  prodigious  computation  times  required  for  the  device 
simulation.  (Computation  times  are  compared  in  Section  3.4.) 

The  new  model  in  Fig.  3.9  includes  basic  LIGBT  elements 
defined  in  Fossum  et  al . [Fos88]  . They  were  derived  by 
regionally  partitioning  the  LIGBT  into  one-dimensional  wide- 
base  bipolar  transistors  coupled  to  the  DMOST . The  ambipolar 
transport  problem  for  high  injection  in  the  transistor  was 
solved  analytically,  but  not  in  closed  form,  to  define  a 
charge -based  module.  The  system  of  module  equations  is 
solved  by  the  circuit  simulator  during  the  nodal  analysis  to 
characterize  all  elemental  currents,  including  the  charging 
currents  dQ/dt,  based  on  the  quasi-static  approximation.  The 
model  extension  presented  here  includes  (a)  a module  for  the 
parasitic  npn  transistor,  including  the  emitter  as  well  as  the 
base  resistances,  (b)  non-quasi-static  charge  partitioning  in 
one  pnp  transistor  module  (the  second  pnp  module  [Fos88]  is 
not  crucial  here) , and  (c)  representation  of  the  effective 


50 


pin  diode  that  forms  during  latch-up.  These  extensions,  needed 
for  the  proper  simulation  of  latch-up,  are  now  described. 

The  current  sources  Ic„.  and  Ig„  model  the  parasitic  npn 
transistor  at  the  cathode  which,  in  latch-up,  is  activated  by 
the  voltage  drop  across  the  base  pinch  resistor  Rvg.  Following 
conventional  bipolar  modeling  [Sze81],  we  describe,  to  first 
order,  the  collector  current  as 

Icn  = IsCXp[<lVBEn/kT]  (3.1) 

and  the  base  current  as 

I Bn  = ICn/Jn  (3.2) 

where  /,  (which  is  proportional  to  Rpb)  and  ii„  are  constants. 
In  (3.1),  the  base-emitter  voltage  VgEn  is  given  by  the  drop 
across  Rpb  minus  the  lateral  voltage  drop  across  R„e  in  the 
emitter  (n+— source) , caused  by  the  emitter  current  (Icn  plus 
I Bn)  plus  the  DMO ST  channel  current  (see  Fig.  3.9) . Note 
that  increasing  R„e,  which  could  be  done  by  lengthening  the 
source  and/or  by  lowering  the  doping  density,  tends  to  improve 
latching  immunity  but  also  results  in  higher  on-resistance 
for  the  LIGBT  (see  Section  3.5). 

As  indicated  in  Fig.  3.9,  the  hole  charge  stored  in 
the  quasi-neutral  base  of  the  pnp  transistor  is  partitioned 
into  two  components,  Qbe  and  Q BC  r associated  with  the 
emitter  and  collector.  The  partitioning,  which  is  physically 
representative  of  bipolar  transistors  [Foss86] , is  defined  by 
integrating  the  carrier  continuity  equations  in  the  base  (see 
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the  Appendix  A)  to  give  the  best  approximation  for  the  non- 
quasi-static  behavior  associated  with  the  distributed  charge 
dynamics.  As  described  in  the  Appendix  A for  high  injection 
and  significant  base  recombination, 

Qbe  — (2/3 )Qb  (3.3) 

and 

Qbc  - (l/3)Qfl  (3.4) 

where  Qb  is  the  total  quasi-static  hole  charge  stored  in  the 
quasi-neutral  base  region  [Fos88] . Interestingly  the  charge 
partitioning  in  (3.3)  and  (3.4)  is  approximately  the  same  as 
that  derived  for  the  case  of  low  injection  in  a uniformly 
doped  base  with  negligible  recombination  [Foss86] . We  stress 
that  the  conventional  quasi-static  model  [Sze81] , in  which 
Qbe  = Qb  an<*  Qbc  ~ 0,  will  give  significantly  erroneous 
simulations  of  latch-up  in  the  LIGBT. 

In  general  (see  the  Appendix  A and  Fossum  et  al . [Fos88]), 
Qbe  and  Qbc  depend  on  both  the  emitter-base  and  base- 
collector  junction  voltages.  This  means,  for  example  in  the 
predominant  forward-active  mode  of  the  pnp  transistor,  that 
the  charging  current  dQ  bc  l dt  cannot  be  represented  by  a simple 
capacitor  [Fos88];  in  fact,  the  component  (dQ Bc / OVeb) dVEBl dt 
reflects  a transcapacitance,  and  the  network  in  Fig.  3.9 
inherently  accounts,  to  first  order,  for  any  nonreciprocity 
and  non-quasi-static  behavior  in  the  device. 


52 


In  a transient  turn-off  simulation  of  the  LIGBT,  dQgc/dt 
is  important  because  it  reflects  the  finite  carrier  transit 
time  in  the  quasi-neutral  base  [Foss8  6]  . Thus,  as  we  will 
illustrate  in  Section  3.4,  it  tends  to  reduce  the  fast  initial 
drop  in  LIGBT  current  by  the  removal  of  the  DMO ST  channel, 
and  it  places  an  additional  transient  drive  on  which 

renders  the  LIGBT  more  vulnerable  to  dynamic  latch-up. 

When  the  LIGBT  latches,  the  effective  pin  diode  between  the 
anode  and  the  cathode  materializes  because  of  the  conductivity 
modulation  of  the  base-collector  junction  of  the  pnp  and  npn 
transistors.  This  pin  diode,  which  actually  defines  the 
latching  characteristics  [Sei87],  is  simulated  in  our  model 
by  a semi-empirical  nonlinear  resistor  Rpin  as  indicated  in 
Fig.  3.9.  The  resistor,  which  effectively  shunts  the  base- 
collector  junction  when  the  npn  current  is  sufficiently  high, 
is  characterized  in  accord  with  quasi-saturation  modeling  of 
the  npn  transistor  [Sze81] . if  we  assume  that  the  electrons 
drift  through  the  junction  space-charge  region  with  saturated 
drift  velocity  va,  then  the  electron  (and  hole)  density  in 
the  region  is  approximately  ICn/qAvs  where  A represents  an 
area  of  the  junction  (predominantly  under  the  DMO ST  gate) 
that  is  modulated.  Consequently  we  assume  that  the  modulated 
region  is  a shunting  resistance  (for  hole  flow)  given  by 

Rpin  — scrV9/ flplcn 


(3.5) 
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where  ir,rr  represents  the  width  of  the  high-field  space-charge 
region,  which  is  assumed  to  be  a constant  extrapolated  crudely 
from  the  equilibrium  depletion-region  width.  The  important 
dependence  in  (3.5)  is  that  on  Ic„ . Below  latch-up,  the  npn 
transistor  is  not  activated  significantly,  and  the  low  ICll 
defines  a large  Rv,n  that  is  inconsequential  in  the  model 
network.  As  If,,  increases,  decreases,  simulating  the 

effective  removal  of  the  junction  potential  barrier  by  the 
hole-electron  plasma.  This  semi-empirical  model  for  the 
effective  pin  diode  is  simple,  but  as  we  show  in  the  next 
section  it  portrays  LIGBT  latch-up  adequately. 

3 . 4 Model  Verification 

LIGBT  test  structures  with  two  different  pnp  base  widths 
(N  B = 78  /'m,  128  //m)  were  fabricated  using  the  dielectric- 

isolation  bipolar -CMOS -DMOS  (BCDMOS)  technology  [Luc88]  at 
AT&T  Bell  Laboratories.  Each  device  is  contained  in  a separate 
silicon  tub  dielectrically  isolated  from  the  substrate  as  well 

as  other  components,  and  all  terminals  are  accessible  from 
the  top  surface. 

Measured  and  SPICE-simulated  (quasi-) static  latch-up 
current -voltage  characteristics  are  shown  in  Fig.  3.10 
The  measured  data  in  Fig.  3.10(a),  taken  with  an  HP-4145 
parameter  analyzer  using  a current-source  driver,  show  an 
abrupt  transition  from  normal  LIGBT  operation  to  the  latched 
pin-diode  mode,  in  which  the  gate  control 


is  lost  (the 


54 


IA 

CmA) 


Fig. 


3.10 


LIGBT 

static 

widths 


iU^rtnt~VOltage  characteristics,  showing 
latch-up,  for  two  different  lateral  base 
(a)  Measured;  (b)  SPICE-simulated. 
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device  cannot  be  turned  off  by  removing  the  gate  voltage) . 
The  simulation  results  in  Fig.  3.10(b)  show  the  latch-up, 
not  so  abruptly  but  representatively.  No  rigorous  model 
parameter  extraction  was  done;  representative  values  [Fos88] 
were  assumed.  Note  that  the  latch-up  onset  current  decreases 
with  decreasing  TT'y,  in  both  Figs.  3.10(a)  and  3.10(b),  and 
that  the  latched-mode  conductance  is  higher  for  the  narrower 
irfl.  These  results  are  due  to  the  dependence  of  the  pnp 
current  gain  /3p  on  TT  Br  which  is  implicitly  modeled  by  the 
system  of  pnp  module  equations  [Fos88]  underlying  IE,  ICr 
and  ]fpi  in  Fig.  3.9.  Decreasing  1VB  increases  lr  (from 
about  0.6  to  1.0  in  this  case),  which  means,  for  a specific 
anode  current  IAr  a larger  fraction  of  it  flows  through 
RvBt  thereby  promoting  the  latch-up  at  lower  I A . The  good 
correlation  between  the  measurements  and  simulations  in  Fig. 
3.10  provides  sufficient  support  to  conclude  that  the  SPICE 
model,  with  the  semi -empirical  in  (3.5),  does  indeed  give 
reasonable  simulations  of  static  latch-up. 

A dynamic  latch-up  simulation,  for  the  case  of  very  high 
blocking  voltage  (\ B = 400  V)  as  discussed  in  Section  3.2, 
is  shown  in  Fig.  3.11.  For  contrast,  the  simulation  of  the 
turn-off  transient  of  the  same  device  but  with  a lower  VB 
(=  50  V)  is  also  shown;  this  turn-off  does  not  result  in 
latching.  In  the  case  of  the  higher  VBr  the  LIGBT  latches 
because  of  higher  charging  currents  flowing  through  RpB . 
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IA/IO  — t (SLICE) 

Fig.  3.11  SPICE -simulated  transient  turn-off  characteris- 
tics for  a gate-voltage  fall  time  of  50  ns  and 
for  two  different  blocking  voltages  ]'B . 
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The  significance  of  the  base  charge  partitioning  (non-quasi- 
static discharging)  in  the  pnp  transistor  is  demonstrated 
in  Fig.  3.12  where  we  plot  the  excess  (above  the  steady 
on-state  bias)  transient  base-emitter  voltage  of  the  npn 
transistor  taken  from  the  SPICE  simulation  of  the  ]'B  = 400 
V turn-off.  For  comparison,  the  component  of  this  excess 
voltage  due  to  the  current  < IQ  DC M flowing  through  n,,B  is 
also  plotted  in  the  figure.  This  component  in  fact  exceeds 
transient  voltages  due  to  the  other  two  discharging  currents . 
Other  simulations  reveal  that  the  maximum  controllable  (on- 
state)  current  for  the  LIGBT  is  overestimated  when  Q Bc  is 
neglected  as  in  conventional  quasi-static  transistor  models. 
Comparison  of  the  transient  simulations  in  Fig.  3.11  with 
the  corresponding  PISCES  simulations  in  Fig.  3.7,  which  show 
qualitative  agreement,  gives  support  for  the  SPICE  model  in 
dynamic  latch-up  simulation.  (No  attempt  was  made  to  tune  the 
SPICE  model  parameters  based  on  the  PISCES  simulations.)  We 
note  further  that  the  SPICE  simulations  required  only  about 
1/600  of  the  computation  time  of  the  PISCES  simulations. 

3 . 5 Summary 

An  insightful  study  of  static  and  dynamic  latch-up  in  the 
LIGBT,  based  on  extensive  PISCES  simulations,  has  revealed  the 
important  underlying  physical  mechanisms  and  has  facilitated 
the  development  of  a physical  SPICE  LIGBT  model  for  HVIC 
simulation.  The  SPICE  model 


was  supported  by  measurements 
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AVBEN  — t (SLICE) 


Fig.  3.12 


The  excess  transient  base-emitter  voltage  (solid 
pb  °f  the  npn  transistor  in  the  VB  = 400  V 

^“Ulati°"  in  Fi9'  3-11'  “d  the  "component 
Of  this  excess  voltage  (dotted  curve)  due  to  the 
partitioned  discharging  current  <IQnc/< '//. 
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of  test  devices  as  well  as  by  the  PISCES  simulations. 
The  study  showed  that  the  latch-up  is  controlled  by  an 
effective  pin  diode  that  materializes  between  the  anode  and 
the  cathode/source  due  to  excessive  conductivity  modulation 
by  the  hole-electron  plasma  near  the  surface  of  the  device.  In 
the  dynamic  case,  the  regenerative  process  in  the  lateral  pnpn 
structure  can  be  triggered  by  either  displacement  current  in 
the  DMO ST  gate-drain  overlap  capacitance  (for  very  fast  gate- 
voltage  switching)  and/or  by  internal  charging  currents  in  the 
constituent  pnp  bipolar  transistor  (for  very  high  blocking 
voltages).  Both  PISCES  and  SPICE  simulations  of  dynamic 
latch-up  in  the  LIGBT  demonstrated,  for  the  first  time,  the 
significance  of  non-quasi-static  bipolar  behavior  (pnp  base 

charge  partitioning)  in  defining  the  maximum  controllable 
on-state  current. 

The  charge-based  SPICE  model  is  useful  in  LIGBT  as  well  as 
in  HVIC  design,  since  the  model  parameters  are  predominantly 
physical,  structural  designs  to  maximize  latch-up  immunity  can 
be  derived  and  checked.  The  mechanisms  by  which  a fast  falling 
gate  voltage  and/or  a high  blocking  voltage  result  in  dynamic 
latching,  subject  to  the  proximity  to  static  latch-up  in  the 
steady  on-state,  can  be  simulated  in  such  designs.  Indeed 
dynamic  latching  under  all  possible  circuit  transients  can  be 
simulated  representatively  and  with  computational  efficiency. 
Such  simulations  can  lead  to  optimal  designs,  for  a given 
technology,  which  could  involve  control  of  not  only  which 
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previous  studies  [Rob86]  emphasize,  but  also  of  R„F.  y;).  j„. 
and  even  the  DM0 ST  properties  [Cho88]  . For  example,  we  have 
done  simulations,  exemplified  in  Fig.  3.13,  which  show  that 
latch-up  can  be  avoided  for  normal  operating  conditions  by 
increasing  the  DMO ST  source  resistance  R,lF  to  values  which, 
albeit  increase  the  on-state  resistance  of  the  LIGBT,  are  not 
infeasible.  In  such  a design,  the  voltage  drop  across  RllF  is 
exploited  to  counteract  that  across  Rl>Dr  and  thereby  reduce 
the  bias  on  the  parasitic  npn  transistor  (see  Fig.  3.9). 


IA[A] 
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Fig. 


IA  - VA  (SLICE) 

3'13  sp ICE -simulated  LIGBT  current -voltage  ch»r»rt„ 
rstrcs  showing  the  effect  of  varying  the  DMO-si 
source  resistance  R„E.  y g tne  DM0ST 


CHAPTER  4 

PHYSICAL  VDMOST  MODELING  FOR  HIGH-VOLTAGE  IC  CAD 

4 . 1 Introduction 

This  chapter  presents  a charge-based  model  for  the  vertical 
DMO ST  (VDMOST)  [Huc84] . As  in  Chapter  3,  the  two-dimensional 
device  simulator  PISCES  is  used  to  gain  the  physical  insight 
the  device  operation  and  to  aid  the  model  development . The 
device  is  regionally  partitioned  into  several  one-dimensional 
components  and  the  carrier-transport  equations  are  solved  in 
each  region.  The  model  is  implemented  in  SPICE  using  FORTRAN 
subroutines . 

DMOSTs  can  now  be  merged  with  low-power  control  circuitry 
to  produce  high-voltage  integrated  circuits  (HVICs) . Effective 
CAD  of  such  HVICs,  and  of  those  containing  IGBTs  described  in 
Chapters  2 and  3,  requires  a physical  DMO ST  model  implemented 
in  a circuit  simulator  such  as  SPICE.  However  the  equivalent- 
circuit  MOST  models  in  SPICE  today  are  too  empirical  and  do 
not  properly  account  for  important  unique  features  of  the 
high-voltage  DMO ST  such  as  quasi-saturation,  space-charge- 
limited  (SCL)  current,  channel  doping-density  gradient,  and 
local  heating.  In  this  chapter,  we  present  a physical  modeling 
methodology  for  DMOSTs  intended  for  HVIC  CAD.  A semi-numerical 
model,  implemented  in  SPICE,  is  developed  for  vertical  DMOSTs. 


62 


63 


Some  of  the  modeling  is  applicable  as  well  to  lateral  DMOSTs 
(LDMOSTs) , which  are  discussed  in  the  next  chapter. 

In  the  methodology,  the  DMO ST  is  regionally  partitioned 
into  two  main  components:  the  double-diffused  MOS  transistor 
(viz.,  the  channel)  and  the  lightly  doped  drain  (drift 
region) . The  constituent  components  are  represented  by  charge- 
based  models  in  which  the  current  and  charges  are  implicitly 
related  to  internal  node  voltages . The  characteristics  of  the 
composite  DMO ST  model  are  derived  semi-numerically  within  the 
nodal  analysis  framework  of  SPICE.  The  model  is  supported  by 
two-dimensional  numerical  simulations  done  with  PISCES  and  by 
measurements  of  test  VDMOSTs.  The  PISCES  simulations  also  aid 
the  model  development  by  giving  physical  insight  regarding 
the  unique  field  distribution  in  the  channel,  the  excess 
carrier  distribution  in  the  drift  region,  and  depletion- 
region  modulation  due  to  free  carriers . 

The  composite  charge-based  VDMOST  model  is  not  a simple 
equivalent  circuit,  but  is  a physical,  semi -numerical 
representation  (in  SPICE)  of  the  underlying  carrier-transport 
equations.  The  critical  elements  are  defined  by  an  implicit 
set  of  nonlinear  equations  that  require  numerical  solution, 
which  is  done  using  a Newton-Raphson  iterative  method  within 
the  model  routine.  In  addition  to  the  constituent  channel  and 
drift-region  models,  the  (n-channel)  VDMOST  model  accounts 
for  the  parasitic  npn  BJT,  which  is  especially  important 
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with  regard  to  the  transient  charge  dynamics  of  the  VDMOST. 
Further,  local-heating  effects  are  included  in  the  electron 
mobility  model,  which  enables  simulation  of  the  negative 
resistance  in  the  DMOST  saturation-region  characteristics . The 
intrinsic  charge  dynamics  are  accounted  for  directly  via  total 
time  derivatives  of  stored  charges,  which  are  characterized 
semi -numerically  by  the  regional  analyses.  SPICE  simulations 
of  VDMOSTs  and  HVICs  are  presented  to  demonstrate  the  model. 

4.2  Model 

The  VDMOST  (shown  in  Fig.  4.1)  model  is  derived  from 
regional  analyses  of  carrier  transport  in  the  channel  and 
regions.  The  parasitic  BJT  is  accounted  for  in  the 
composite  model  by  physically  based  current  elements,  which 
can  be  activated  by  the  lateral  voltage  drop  across  the 
p— body  resistance. 

4.2.1  Channel  Region 

The  double-diffused  n-MOS  transistor  has  a nonuniformly 
doped  [-V-i  (.r)]  channel  region  defined  by  the  compensation 
between  the  p-body  and  n-source  lateral  diffusions.  -Y,(.,)  has 
a maximum  (iV4„)  near  the  source  and  decreases  sharply  toward 
the  dram.  Consequently,  as  quantified  by  PISCES  simulations 
for  Vqs  > 0,  the  inversion  charge  is  minimum  at  the  source, 
and  hence  the  longitudinal  electric  field  is  maximum  there 
(Irh  oc  ExQn)  . Fig.  4.2  shows  PISCES-simulated  equipotential 
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Fig.  4.1  VDMOST  structure.  The  partitioning  of  the  drift 
region  is  shown  and  geometrical  parameters  are 
defined.  The  common  deep  p+— body  region,  which 
does  not  influence  our  regional  analyses,  is  not 
shown . 


66 


PISCES-simulated  equipotential  contours  (AT  = 0 1 

(V  in_tco  7rlcdinity  of  the  ^OST  [LUC88]  channel 
O/rs;  - 50  V,  ] r;s  = 10  V)  . On  the  and  ,/-axes 
1 division  = 0.1  ,,m.  ’ es' 


Fig.  4.2 
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contours  in  the  channel  of  a typical  DMO ST  [Luc88] . The 
high  density  of  lines  near  the  source  indicates  the  high 
electric  field  there.  Our  channel  model  is  representative 
of  these  properties,  and  involves  an  appropriate  gradual- 
channel  approximation  that  yields  the  channel  current-voltage 
characteristic,  with  the  threshold  voltage  defined  by  iV  \o  and 
the  saturation  controlled  by  drift-velocity  saturation  at  the 
source  (rather  than  at  the  drain  as  in  the  conventional  MOST) 

We  approximate  the  nonuniform  doping  density  along  the 
channel  as  an  exponential  function  [XA  (.v ) = -i].r/L )] , 

and  we  use  a piecewise-continuous  electron  velocity  model 
[Sod84]  : 


ntf  f 

1 + ( Une  f fEr  / 2 V9ai ) 

= v.wi  otherwise  . 


for  v < v$al 


(4.1) 


In  (4.1)  Er  — (J\  /dx  is  the  longitudinal  electric  field 
(magnitude),  (~  10'  cm/s)  is  the  saturated  drift  velocity, 
and  I'veff  = /<»n/[l  + 0Oc;5  - VV )J  where  is  the  low-field 

mobility,  h is  the  fitting  parameter  which  defines  the 
transverse— field  dependence,  and 


1 T = Vfb  + 2 \&b  \ + 


IobII172 

or 


(4.2) 


is  the  threshold  voltage.  For  strong  inversion  and  negligible 
short-channel  effects,  the  quasi-static  channel  current  is 
described  by  solving  for  E,  in  Ich  = -WxVQ„  (TH  is  the  device 
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width  and  Q„  is  the  channel  charge  density)  and  integrating 
it  from  source  (.r  = 0)  to  drain  (r  = Z): 


lrh  - ^ :lln<  f f 

V idneff/-t'satL)  Irft]  ( C or 

where  1 r/,  is  the  channel  voltage 
capacitance, 


[Ql\L)-  Qj;(0)]  (4.3) 

<-^roP  • The  body  depletion 


Cd 


<JQd^ 

(I  V - 


1/2  0 _ '2>i\?b\\ 

aWd\\  V Vch  ) 


(4.4) 


is  assumed  to  be  invariant  along  the  channel  to  obtain  a 
first-order  characterization  of  the  depletion  charge  Qd  in 
terms  of  the  voltage  ( (~  Yr,t.v/L)  along  the  channel.  In  (4.3), 

Qn(L)  = -C'ox  [I as  — Vt-  Vch  ( 1 + C\i/ Cor )]  (4.5) 

and 

Qn  (0)  = — C‘or  (I  GS  — I'r)  • (4  6) 


Combining  (4.3) -(4. 6),  we  get  the  channel  current-voltage 
characteristic  in  the  triode  region  of  operation: 


Irh  W:l'vrff 

2L  l1  + {t'nfff/2r,(llL)  rf/lJ  (h 

[2C„r  ( \ GS  -Vt)-  (C0T  + 0/(, ) 1 + 2/, 0,(1  \oB\\ 


where  0/o  = \/^X. W?R • The  effect  of  the  channel  doping- 
density  gradient  is  reflected  in  (4.7)  by  when  /,  = 0, 

(4.7)  simplifies  to  the  conventional  MOST  expression.  Note 
in  (4.5)  that  0/  is  typically  negative,  and  much  larger  in 
magnitude  than  C„, . This  unusual  "depletion  capacitance" 
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?/  > 0,  simulates  the  added  channel  conductance  due  to 
the  smaller  depletion  charge  away  from  the  source. 

As  VDS  increases,  ]'rh  increases,  and  ultimately  the  electron 
velocity  saturates  near  the  source  (at  c=0)  as  characterized 
by  (4.1) . Then  the  saturation  current  is  expressed  as 

Ich{.iu1)  ^ r I'sat^Ju  (•*’  — 0) 

= N s l'sa  I Cox  ( I as  -Vt)  . (4.8) 

Unlike  in  the  conventional  MOST,  the  saturation  current  in 
(4.8)  is  independent  of  C\i  and  cannot  be  modeled  by  a channel 
pinch-off  approximation.  Combining  (4.8)  with  (4.7)  yields 
^ rh{.mi)  for  a given  VGS,  above  which  I,.h  = Ic , . 

4.2.2  Drift  Region 

Depending  on  the  VDMOST  structure,  the  (quasi-)  saturation 
characteristics  can  be  controlled  by  the  drift  region  instead 
of  the  channel  [Dar86]  . Whereas  the  channel  model  is  applicable 
to  both  vertical  and  lateral  DMOSTs,  as  well  as  to  the  IGBT, 

different  drift-region  models  are  needed  for  the  VDMOST  and 
the  LDMOST . 

For  the  VDMOST  in  the  normal  mode,  the  PISCES  simulations 
(Fig.  4.3)  reveal  that  the  drift  region  can  be  partitioned 
into  three  sections  as  shown  in  Figs  4.1  and  4.3  to 
simplify  the  analysis.  The  voltage  drops  in  the  sections, 
which  are  characterized  semi -numerically  based  on  reasonable 
approximations  for  the  carrier  transport,  are  added  to  define 
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^SC?S7simulated  current  flowlines  (AI  = 9 7xl0“6 
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the  composite  quasi-two-dimensional  model.  In  section  A, 
excess  electrons  supplied  through  the  channel  (for  VGS  > \j) 
accumulate  to  support  the  negative  gradient  of  the  vertical 
electric  field  for  ] ns  > 0 [Dar86]  . These  excess  electrons, 
which  accumulate  even  before  the  drift  velocity  saturates, 
produce  SCL  current.  The  underlying  Poisson  equation  is 
complex,  but  can  be  solved  analytically  to  give  an  implicit 
characterization  of  the  voltage  drop  across  the  region.  In 
section  B,  the  electron  current  spreads  laterally  with  an 
angle  n (~  45°)  which  can  be  assumed,  in  accord  with  the  PISCES 
simulations,  to  be  independent  of  the  bias  condition  and 
the  structural  geometry  [Huc84].  The  boundary  locations  are 
defined  by  the  depletion  region  of  the  p-body  junction,  which 
varies  with  the  voltage  drop  across  section  A.  The  region  is 
assumed  to  be  quasi-neutral  (//  ~ ND ) . Section  C is  also  neutral 
but  has  a constant  cross-section.  The  (moving)  boundary 
between  sections  B and  C is  also  controlled  by  the  depletion- 
region  width  of  the  p-body  junction.  Conventional  MOST  models, 
which  do  not  account  for  SCL  current,  underestimate  the  drain 
current.  Our  simulations  show  that  this  underestimation  can 
be  quite  significant  at  high  YJ)S . 

In  each  section,  the  drain  current  is  generally  expressed 
using  a continuous  version  of  (4.1) : 


n ii  it 


1 + EJEC 


In  = Ayqn 


for  Ey  > 0 


(4.9) 
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where  Ec  = raai/nnd  and  //,„<  is  the  low-field  electron  mobility. 
Combining  (4.9)  with  Poisson's  equation  in  one  (vertical) 
dimension  yields 


In  (4.9)  and  (4.10),  Av  is  the  cross-sectional  area  at  depth 
U and  E,,  = dV/dy  is  the  magnitude  of  the  (vertical)  electric 
field. 

In  section  A,  Ay(=WzL,i)  is  assumed  to  be  constant.  This 
assumption  is  based  on  the  PISCES  simulations,  which  reveal 
a slight  JFET  effect  in  section  A at  low  Vus  but  show  it 
compensated  at  higher  ] due  to  predominant  excess  electrons 
[Dar86] . Equation  (4.10)  is  solved  analytically  by  extending 
the  analysis  in  [Lam70]  (which  does  not  account  for  velocity 
saturation)  to  yield  expressions  for  the  position  y and  the 
voltage  1 at  y in  terms  of  Eu . The  nonlinear  y(Ev)  and  V(E„) 
expressions  then  implicitly  give  the  section  voltage  drop. 


which  requires  numerical  solution;  U'7  is  defined  in  Fig. 
4.1  and  TT(/  is  the  depletion  region  width  of  the  body-drift 
junction.  Note  that  the  depletion-region  width  depends  on 


Eu  is  everywhere  increasing  with  .?/,  so  the  velocity  of  the 
electrons  increases  along  the  vertical  direction  from  the 
surface.  However  the  excess  electron  density  decreases  with 


(4.10) 


I (y  — IT j + IT ,( ) = . 


(4.11) 


0’i  + l>*),  the  effective  reverse  bias  on  the  p-body  junction. 
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increasing  y because  the  gradient  of  £„  becomes  smaller. 
In  our  analysis,  we  neglected  the  diffusion  component  of 
current,  which  is  confirmed  to  be  negligible  by  both  PISCES 
simulations  and  analytic  solutions. 

In  sections  B and  C,  n ~ i\  d . (The  error  resulting  from 
the  quasi-neutrality  assumption  is  insignificant;  results 
of  accurate  analytic  solutions  show  little  difference  with 
those  obtained  using  the  approximation.)  With  A „ = 
TT-  [L,i  + (y  - \Yj  - ]V,i)rota]  in  section  B and  constant  in  section 
C,  (4.9)  is  solved  for  Ev,  and  the  result  is  integrated  over 
the  thicknesses  of  sections  B and  C to  obtain 


Vb  + Vc 


[d 

l J (1  N Dl'nd  cot  O 


+ 


Ip 

M - ( Ld  + Lp)qNi)li„,i 


(L,i  + Lt, ) — Ip/Er 
W.qNDfivdLd- Id/Ec 

Ip/Er  _ Lptano).  (4.12) 


The  geometric  parameters  in  (4.12)  are  defined  in  Fig.  4.1. 
Note  that  current  continuity  is  effected  in  section  B by  the 
electron  velocity  decreasing  with  increasing  cross-sectional 

area.  This  occurs  within  the  realm  of  the  quasi-neutrality 
condition . 

Now  the  total  drift-region  voltage  drop  of  the  VDMOST  is 


^ drift  — l .1  + 1 B + Ec  ■ (4.13) 

For  a given  current  Iu,  which  equals  /,,,  from  the  channel- 
region  analysis,  (4.13)  must  be  solved  numerically  because 
of  the  implicit  nature  of  in  (4 . 9) - (4 . 11) . The  Newton- 

Raphson  iterative  method  is  used  (in  the  model  routine)  to 
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solve  these  equations.  Note  that  the  voltage  drop  in  the 
accumulation  layer  under  the  extended  gate,  which  varies 
inversely  with  \\;s,  has  been  neglected  in  (4.13). 

4.2.3  Parasitic  BJT 

The  parasitic  npn  BJT  in  the  n-channel  VDMOST,  (up-mode) 
activated  in  the  PISCES  simulations  of  Fig.  4.4,  plays  an 
important  role  in  both  normal-  and  reverse-mode  operations. 
In  the  normal  mode,  where  a positive  voltage  is  applied  to 
the  drain,  the  npn  BJT  near  the  channel  (down-mode:  see  Fig. 
4.4)  can  be  activated  either  statically  (by  impact-ionization 
current)  or  dynamically  (by  displacement  or  charging  current) 
due  to  the  p-body  voltage  drop.  (in  typical  devices,  this 
voltage  drop  is  minimized  by  using  a deep  p+-body  diffusion  to 
reduce  the  body  resistance  Rp0.)  m the  reverse  mode,  where  a 
negative  voltage  is  applied  to  the  drain  (which  could  happen 
during  transient  switching  of  the  VDMOST  in  an  HVIC) , the 
BJT  near  the  body  contact  (up-mode:  see  Fig.  4.4)  becomes  an 
essential  part  of  the  VDMOST  operation . The  reverse  recovery 
is  very  slow  due  to  large  stored  charge  in  the  drift  region 
(where  high  injection  prevails) , and  can  sometimes  drive  the 
BJT  to  saturation,  which  would  lengthen  the  recovery  and 
could  possibly  lead  to  destructive  breakdown  [Bal87]  . 

The  npn  BJT  is  modeled  following  conventional  bipolar 
transistor  analysis  for  the  normal  mode,  and  accounting  for 
high-level  injection  in  the  drift  region  for  the  reverse 
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Up-mode  BJT  Down-mode  BJT 


Fig.  4.4 


PISCES -simulated  current  flowlines  (A/  = 2.4xl0-3 

thu-  3fevefse-mode  operation  of  the  VDMOST 
[Luc 8 8 ] , which  show  the  activation  of  the  up- 

mode  BJT  {\DS  = -2.0  V,  Vas  = 0 V).  The  regions 
functioning  as  the  emitter  and  collector  of^oth 
the  up-mode  and  down-mode  BJTs  are  indicated.  On 
the  .i—  and  //—axes,  1 division  = 1.0  //m. 
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mode.  The  base  transport  current  is  characterized  in  terms 
of  the  internal  junction  voltages : 


Ict  - IcCi)C*'P  ( 1 BEn/]  T ) — lEC'0?rp(2VBCn/VT)  (4.14) 


where  two  different  saturation  currents  are  used  because  the 
BJT  action  takes  place  at  different  locations  depending  on 
the  mode  of  operation.  In  the  reverse  mode  of  the  VDMOST, 
the  body-drift  region  junction  is  forward-biased  (Vue,  > 0) 
predominantly  near  the  body  contact;  and  I ecu  in  (4.14)  is 
relevant.  In  the  normal  mode,  the  source-body  junction  can  be 
forward-biased  {VBt:n  > 0),  predominantly  near  the  channel,  due 
to  the  lateral  voltage  drop  across  the  body  resistance;  and 
/cco  in  (4.14)  is  relevant.  We  neglect  the  BJT  base  current 
for  the  latter  case.  For  the  reverse  mode,  it  is  defined 
predominantly  by  the  hole  transport  in  the  drift  region. 

The  transport  mechanism  in  the  drift  region,  where  high- 

level  injection  occurs  in  the  reverse  mode,  is  characterized 

by  the  ambipolar  transport  equation.  For  p ~ „,  the  solution 
[Fos88]  is 


P(y)  = [or,  - !I)/La\  + p(W,)*h,h  \(y  - W,)/L  ,1 


(4.15) 


sink  [(IF,  - T Vj)/LA] 

where  LA  = ^Dat„  is  the  ambipolar  diffusion  length,  DA  is 
the  ambipolar  diffuaivity,  and  r„  is  the  high-injection 
carrier  lifetime.  The  boundary  carrier  densities  are 

;><"•,)=  (nl/ND)c.rPiVtt€n/VT)  and  p(ir,)  = XD(r,,(VJn.„t/Vr)  where 
iS  the  voltage  drop  across  the  low-high  junction  at  the 
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drain,  which  is  solved  from  the  following  boundary  conditions. 
The  hole  and  electron  current  expressions  are  combined  to  give 
two  boundary  conditions,  neglecting  recombination  in  the  p- 
body  and  n+-drain  regions:  /„({/  = TT>)  ~ -ICT  and  Ip(y  = TT',)  ~ 0 

[McD85] . For  given  \ BEn  and  I oCnr  these  boundary  conditions 
are  solved  to  give  p(Wj)  and  ;>(!!',)  in  (4.15).  Then,  the 
reverse-mode  base  current  is 


where  QP  is  the  hole  charge  stored  in  the  drift  region: 

u , 

Qp  = <flVzLp J p(y)dy  . (4.17) 

The  voltage  drop  across  the  drift  region  in  the  reverse 
mode  can  be  obtained  by  integrating  the  electric  field,  using 
P(y)  in  (4.15),  and  by  adding  Vjn-„+ : 


where  E(y)  is  related  to  p{y)  through  a combination  of  the  hole 
and  electron  current  expressions  [Fos88] . The  integration 
in  (4.18),  done  analytically,  yields  Vdrift  as  a function  of 

p{lVj)  and  />(Tr,). 

The  VDMOST  charge  dynamics  are  defined  predominantly  by 
Qp  in  (4.17)  in  the  reverse  mode.  In  the  normal  mode,  the 
body-drift  region  junction  depletion  charge. 


Ibr  = Qp/ry 


(4.16) 


(4.18) 


Q.J  - qWzLvXD\Yd  . 


(4.19) 
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is  significant,  as  well  as  charge  associated  with  the  gate 
discussed  in  the  next  section. 

4.3  SPICE  Model 

Based  on  the  analysis  of  the  VDMOST  in  Section  4.2,  we 
propose  a physical  semi -numerical  SPICE  model,  the  basis  of 
which  is 

^ DS  ^ ch  ^ drift  • (4.20) 

where  V is  given  either  by  (4.13)  or  (4.18)  depending  on 
the  mode  of  operation.  The  complete  network  representation 
of  the  new  SPICE  model  is  shown  in  Fig.  4.5.  The  model 
is  not  a simple  equivalent  circuit.  The  elements  shown 
represent  physical  components  of  current  and  voltage  in  the 
VDMOST,  described  implicitly  by  the  analysis,  which  requires 
numerical  solution.  Note  in  Fig.  4.5  that  VBEn  (the  voltage 
drop  across  the  p-body  resistance  R,)R)  and  VBCn  (=  -Ych ),  which 
control  ICT  and  IBB,  are  defined  in  accordance  with  the 
analysis  of  the  parasitic  BJT . 

The  transient  behavior  of  the  device  is  simulated  using 
both  the  quasi-static  currents/voltages  and  charging  currents 
(IQ /dt.  The  charges  within  the  device  are  expressed  in  terms  of 
node  voltages  as  described  in  Section  4.2,  and  the  transient 
charging  currents  are  represented  by  the  time-derivative  of 
the  charges,  using  the  quasi-static  approximation: 

<]Qi_  v 

' bVj  (It 


(It 


(4.21) 
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Body/Source 


Cfls  T ■oh  (P 
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-VBEn+ 

— vw- 
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vBCn©IBR  ©dQp/dt  ©dQj/dt 


dQoD 

dt 


“Voh+ 


Olci 


0+vdrift 
Drain  1 


Fig.  4.5  Network  representation  of  the  charge-based  VDMOST 
model . 
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where  i = P,  ./,  and  GD  and  j = GS,  rh , BEn.  Included 

in  (4.21)  is  the  DMOST  gate-drain  charging  current  dQ^o/dt. 
This  current  arises  from  the  charge  near  the  surface  of 
the  drift  region  under  the  gate  electrode  and  is  important 
in  the  switching  performance  of  the  VDMOST . The  associated 
gate  charge  Qgq  is  characterized  by  basic  MOS  theory  for  two 
different  modes  of  operation.  Whenlcs-lrfc  — VfB  > 0/  the  surface 
is  accumulated  with  electrons,  and  Qcw  ~ L(,WzC'ox  ( Vas  ~ Kh  ~ VFW) 
with  the  approximation  that  the  surface  potential  i.'s  ~ 0.  When 
t as  ~ ^ ch  ~ ^FB  < 0,  only  depletion  charge  can  exist  under  the 
gate  electrode  due  to  the  usual  high-frequency  operation  of 
the  device.  In  this  case,  \qs  — ^ ch  ~ ^ FD  — 'i'S  + Qcw/C'or  and 
Qg'D  = ~ L,/]]  s yj2(ieNj)  ( —u's)  define  the  charge. 

The  effective  gate-drain  capacitance  associated  with 
•IQcn/M  is  amplified  by  the  Miller  effect  into  an  equivalent 
input  gate  capacitance,  thus  slowing  down  the  gate-controlled 
switching.  Furthermore,  (IQgd/^G  with  gate-circuit  resistance 
can  cause  dV/dt- induced  turn-on  [Bal87]  . When  VDS  is  abruptly 
increased  while  1 as  < 1 Tr  the  charging  current  through  the 
gate  resistance  can  increase  the  gate  potential  above  VT,  thus 
turning  on  the  VDMOST.  Also  included  in  the  model  is  Ca s' , the 
gate-source  capacitance,  the  most  significant  component  of 
which  is  that  due  to  gate  electrode  overlapping  the  n+— source. 

All  of  the  elements  in  the  model  are  implemented  in 
SPICE2  via  user-defined  controlled  sources  (UDCSs) , which  are 
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FORTRAN  subroutines  in  SLICE  [Har84] , an  enhanced  version  of 
SPICE2.  Each  UDCS  accesses  the  VDMOST  model  routine,  which  is 
flowcharted  in  Fig.  4.6.  The  model  routine  numerically  solves, 
usincf  the  Newt on—Raphson  iteration  method,  the  implicit  system 
of  quasi-static  equations  developed  in  Section  4.2.  The  UDCSs, 
for  a given  Newton-Raphson  iterate  on  VGS.  Vch , and  VBE„  in 
the  SPICE2  nodal  analysis,  provide  the  element  currents,  or 
voltage  for  and  the  needed  partial  derivatives  with 

respect  to  the  node  voltages . The  derivatives  are  calculated 
via  difference  equations,  based  on  multiple  calls  of  the 
model  routine.  In  the  charging-current  UDCSs,  backward 
Euler  method  for  numerical  integration  is  used  to  get 
best  numerical  stability.  The  semi-numerical  VDMOST  model, 
implemented  in  SPICE,  provides  then  a two-level  Newton  method 
for  device/circuit  simulation.  Note  that  the  model,  once 
finalized,  can  be  written  directly  into  the  SPICE2  source 
code  by  modifying  the  algorithm  in  Fig.  4.6  to  accept  the 

device  terminal  voltages  as  input  from  the  circuit  nodal 
analysis . 

The  negative  resistance  often  observed  in  I,)S  - Vns 
measurements  of  the  VDMOST  is  mainly  due  to  the  local 
heating  of  the  device.  When  } 'Ds  is  increased  for  a given 
\GS  in  the  saturation  region,  the  local  power  dissipation 
elevates  the  device  temperature,  and  concomitantly  reduces 
the  carrier  mobility  as  well  as  the  saturation  velocity; 
hence  IDS  decreases.  If  the  duty  cycle  for  the  pulsed-mode 


82 


Given  VGS  ,Vch,and  VBEn 
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Semi -numerical  VDMOST  model-routine  flowchart.  Both 
norma1 -and  reverse-mode  paths  are  shown  within  the 
outer  iteration  loop  for  local  heating 


Fig.  4.6 
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measurements  is  reduced,  less  heating  occurs  and  consequently 
the  negative  slope  in  the  ID$  — ]'DS  characteristics  tends  to  be 
eliminated.  To  account  for  this  effect  for  general  DMOST/HVIC 
simulation,  we  implement  the  concept  of  thermal  resistance 
in  the  SPICE  model.  The  device  junction  temperature  (T;)  is 
different  from  the  ambient  temperature  ( Ta ) as  follows: 

Tj  = Ta  + n, I,  \ dsIds  (4.22) 

where  Rlh  is  the  thermal  resistance  of  the  device.  Empirical 
relationships  between  mobility  and  saturation  velocity  and 
temperature  [Mul77]  are  used  to  implement  (4.22)  in  parameters 
I1""'  l1”*'  and  v™t  • This  implementation  (which  does  not  include 
the  BJT  temperature  dependence)  requires  another  iterative 
solution  around  the  basis  model  routine  as  shown  in  Fig.  4.6. 
It  can  be  superseded  if  local  heating  is  negligible. 

4 • 4 Model  Verification 

Test  VDMOSTs  are  fabricated  at  AT&T  Bell  Laboratories  using 
the  dielectric  isolation  (DI)  BCDMOS  technology  [Luc88] . Our 
initial  measurements  were  done  using  a very  low  duty  cycle 
(1%)  to  avoid  the  negative-resistance  effect  in  the  Ius  - yl)s 
characteristics.  The  duty  cycle  is  a critical  parameter 
in  the  measurements  as  we  note  below  in  the  discussion  of 
local-heating  simulations . 

in  Fig.  4.7,  measured  dc  IDS(Vos.  rcs)  data  for  two 
VDMOSTs  having  drift-region  widths  L,,  are  compared  and 
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Fig.  4 


.7  Simulated  (curves)  and  measured  (points)  VDMOST 
current-voitage  characteristics  for  two  different 
drift-region  widths  L,i . 
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show  good  agreement  with  model-simulated  current-voltage 
characteristics . The  simulated  characteristics  were  derived 
from  SPICE  with  our  new  DMO ST  model,  using  known  structural 
device  parameters  and  crudely  extracted  model  parameters  that 
fit  the  measured  data.  The  channel  length  L is  2 /im.  The 
device  width  TJ\  was  taken  as  the  sum  of  sub-cell  widths  and 
thus  effectively  accounts  for  the  device  layout.  The  key 
model  parameters  which  determine  the  (channel)  saturation 
characteristics  are  jV.i„,  and  vsnt  in  the  channel  region. 
Parameters  A and  r„(7  in  the  drift  region  are  the  key 
parameters  for  the  quasi-saturation  characteristics.  Quasi- 
saturation, reflected  by  no  dependence  of  the  drain  current 
on  the  gate  bias  and  by  high  output  conductance  without  the 
current  saturation,  is  predominant  for  the  Ld  = 4 /in  device. 

In  Fig.  4.8  we  show  SPICE-simulated  transient  character- 
istics of  the  VDMOST  with  an  inductive  load.  As  the  gate 
voltage  is  dropped  to  turn  off  the  device,  the  drain  voltage 
Od.v)  is  abruptly  increased  above  the  off-state  blocking 
voltage  due  to  the  negative  voltage  drop  across  the  inductor 

(Li) ' Due  to  this  ^ick  Crease  of  VDS,  the  charging  current 
of  the  depletion  region  of  the  body-drift  region  junction 
(dQ. l /(It)  flowing  through  RllB  produces  a voltage  drop  VBh„, 
which  can  turn  on  the  down-mode  npn  BJT.  In  addition  to  Lf , 
the  load  resistance  RL  also  affects  the  turn-off  time.  As 
UL  increases,  dQj/df  slows  down  due  to  small  drain  current; 
thus  the  turn-off  time  is  increased. 


CAD  G)A 
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Simulated  inductive  turn-off  transients  for  the 
L,i  - 10  /i m VDMOST . The  inductor  (Li  = 5 //H)  in 
series  with  RL  _ (=  250  12)  is  connected  to  the  drain 
from  the  blocking  voltage  source  {Y^  = 20  V),  and 
\c;s  is  pulsed  as  indicated.  The  simulated  1 'Ds  ( / ) 

and_  \nEn(t),  the  voltage  drop  across  the  p-bodv 
resistance,  are  shown.  J 


Fig.  4.8 
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Simulated  reverse-recovery  characteristics  of  the  VDMOST 
is  shown  in  Fig.  4.9.  Initially  the  device  is  in  the  reverse 
mode,  in  which  the  up-mode  BJT  is  activated  and  a high 
concentration  of  carriers  (holes  and  electrons)  is  stored  in 
the  drift  region.  Therefore  the  reverse  recovery  is  slow, 
controlled  by  the  carrier  lifetime  ( 77/  = 0.2  //.s)  in  the  drift 
region.  The  device  is  then  switched  from  the  reverse  mode 
to  the  normal  blocking  mode.  For  RpB  = 0,  the  simulation 
shows  about  a 0.1  Bs  delay.  But  for  R)tB  = 2.5  ohm,  the 

down-mode  BJT  is  turned  on  by  the  voltage  drop  across  Rj)B 
{\<lQf>/dt\npH>  0.7V)  1 the  delay  is  longer  (0.2  /is).  For  larger 
nvn  = 10  ohm,  the  BJT  becomes  heavily  saturated,  and  the 
VDMOST  does  not  turn-off  for  times  up  to  1 /is.  It  will 

eventually  turn  off  unless  destructive  breakdown  occurs. 

Figure  4.10  shows  Ids  — Vds  characteristics  of  the  VDMOST 
measured  with  a higher  duty  cycle  (3%)  , which  show  local- 
heating effects.  The  simulated  characteristics,  also  shown 
^-n  the  figure,  were  fitted  to  the  measurements  by  varying 
the  parameter  Rlh  in  (4.22).  As  VGS  and  IDS  increase,  the 
local  heating  becomes  more  severe  and  the  negative  resistance 
becomes  more  predominant  as  indicated  by  the  measurements  and 
predicted  by  the  simulations.  Note  that  R,h  tends  to  be  higher 
for  DI  devices  than  for  junction-isolated  devices. 
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Fig.  4.9  Simulated  reverse-recovery  current  transients  for 
the  L,i  = 10  / im  VDMOST  with  different  p-body 

resistances.  The  drain-source  voltage,  applied 
through  a 10-ft  series  resistor,  changes  from  —5  V 
to  +4  V with  a 0 . l-/,s  rise  time  starting  at  0.1  // s. 
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vds  cv: 


Fig.  4.10  Simulated  (curves)  and  measured  (points)  VDMOST 
current-voltage  characteristics  which  show  local- 
heating effects.  The  duty  cycle  used  in  the 
measurements  is  3%  and  the  extracted  thermal 
resistance  /?,/,  is  0.5. 
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4 . 5 Summary 

The  VDMOST  model  developed  herein  is  useful  for 
HVIC  CAD,  which  requires  reliable  physical  device/circuit 
simulation.  Some  of  the  modeling  is  applicable  to  LDMOSTs 
as  well.  The  charge-based  models  are  quasi-static;  non- 
quasi-static  (NQS)  effects,  which  can  be  significant  in 
HVICs,  could  be  accounted  for  via  physics-based  extensions 
[Fos88] . Comprehensive  transient-model  verification  based  on 
device/circuit  measurements  would  be  required  to  ascertain 
the  significance  of  NQS  behavior.  The  models  are  semi- 
numerical,  an  unavoidable  trait  of  good  physical  models,  but 
are  implemented  in  SPICE.  The  UDCS  implementation  described 
is  CPU-intensive  and  is  intended  only  for  preliminary  SPICE 
model  verification.  Indeed  once  the  models  are  finalized, 
they  can  be  written  directly  into  SPICE  source  code,  creating 
a semi-numerical  device/circuit  simulator  for  MOS  HVICs. 


CHAPTER  5 

NEW  PHYSICAL  INSIGHTS  AND  MODELS  FOR  LDMOST  IC  CAD 

5 . 1 Introduction 

This  chapter  presents  a physical  composite  model  for 
lateral  DMOSTs  (LDMOSTs)  [Cla86] . The  modeling  methodology 
developed  in  previous  chapters  culminates  in  this  chapter: 
the  physical  insights  gained  from  PISCES  simulations  are 
used  as  a basis  for  developing  the  composite  LDMOST  model, 
in  which  the  device  is  regionally  partitioned  into  several 
one-dimensional  components,  and  carrier-transport  equations 
in  each  region  are  solved.  The  model  is  implemented  in  SPICE 
using  FORTRAN  subroutines  and  is  supported  by  measurements. 
Device  and  circuit  simulations  show  the  superiority  of  the 
model  over  equivalent-circuit  models. 

High-voltage  integrated  circuits  (HVICs)  including  LDMOSTs 
are  now  viable  alternatives  to  hybrid  circuits  comprising 
discrete  high-voltage  devices.  Optimal  design  of  these  HVICs 
requires  physical  LDMOST  models  implemented  in  a circuit 
simulator  like  SPICE.  Equivalent- (sub) circuit  SPICE  models 
for  the  high-voltage  LDMOST  [Wil87]  have  some  utility,  but  are 
fundamentally  deficient;  they  cannot  simulate  unique  power 
device  characteristics  such  as  quasi-saturation  and  reverse- 
recovery  transients. 
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Several  physical  investigations  of  the  on-resistance  of 
LDMOSTs  [C0I8I,  Men86,  Poc76]  have  been  done,  but  a complete 
physical  model  has  not  been  reported.  LDMOSTs  in  HVICs  are 
subjected  to  unusual  operating  conditions,  and  hence  robust 
models  are  needed.  For  example,  an  inductive  load  in  power 
circuits  will  force  the  LDMOST  into  high  current /voltage 
regimes  where  reverse  current  flow  through  the  device  may 
occur.  Claessen  and  van  der  Zee  [Cla8  6]  presented  a model 
for  the  LDMOST  covering  the  high  current /voltage  region.  But 
their  model  is  limited  to  dc  conditions,  and  it  does  not 
account  for  quasi-saturation  nor  the  inherent  BJT,  which  is 
commonly  activated  in  HVICs. 

In  this  chapter  we  develop  a new  physical  charge-based 
LDMOST  model  which  accounts  for  the  BJT.  We  study  both 
the  lightly  doped-drain  (LDD)  and  the  reduced-surface- 
field  (RESURF)  n-channel  LDMOST  structures,  using  the  two- 
dimensional  device  simulator  PISCES  [Pin84]  for  guidance. 
The  PISCES  simulations  provide  new  insights  on  normal- 
and  reverse-mode  dc  operation  of  the  LDMOSTs  and  on  the 
significance  of  the  inherent  npn  BJT  in  transient  reverse- 
recovery  characteristics.  We  emphasize  the  LDD  LDMOST 
structure  (see  Fig.  5.1(a)),  which  we  regionally  partition 
into  three  components:  the  channel,  the  drift  region,  and 
the  npn  BJT.  The  composite  model  is  implemented  in  SPICE 
for  HVIC  CAD.  The  modeling  methodology  is  also  applicable  to 
the  RESURF  LDMOST  structure  (see  Fig.  5.1(b)),  and  with  our 
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Fig.  5.1  Unit-cell  structure  of  (a)  the 
RESURF  LDMOST . 


LDD  and  (b)  the 
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physical  VDMOST  modeling  in  Chapter  4 defines  LDMOST  models 
for  either  technology. 

5 . 2 New  Insights  By  Numerical  Simulations 
PISCES  [Pin84 ] simulations  have  been  done  to  gain  physical 
insight  regarding  both  the  LDD  and  the  RESURF  LDMOST 
structures.  This  new  insight  facilitates  the  device  modeling 
for  circuit  simulation  described  in  the  next  section. 
PISCES  is  a two-dimensional  numerical  semiconductor  device 
simulator,  for  either  dc  steady-state  or  transient  bias 
conditions.  It  solves  the  basic  semiconductor  device 
transport  equations,  using  plausible  physical  models.  For  the 
LDD  LDMOST,  we  used  doping-density  profiles  and  structural 
dimensions  obtained  from  Siliconix,  Inc.  [Wil87] . The  LDMOSTs 
were  biased  into  both  normal  mode  0 DS  ^ 0)  and  reverse  mode 
07xs  < 0)  in  the  dc  simulations;  the  LDMOSTs  were  switched 
from  the  reverse  mode  into  the  normal  mode  in  the  transient 
reverse-recovery  simulation. 

5.2.1  LDD  LDMOST 

The  self -isolated  LDD  LDMOST  is  shown  in  Fig.  5.1(a). 
The  channel  region,  which  functions  as  a fast  switch  for  the 
electron  current  flow,  consists  of  two  different  regions:  the 
double-diffused  channel  and  the  substrate  channel.  The  double- 
diffused  channel  is  nonuniformly  doped  by  the  difference  of 
lateral  diffusions  between  the  p-body  and  the  n+-source.  The 
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substrate  channel  has  constant  doping  density,  defined  by  the 
substrate  doping.  PISCES  simulations  show  that  the  electron 
velocity  saturates  near  the  source  in  the  double-diffused 
channel,  not  near  the  drain  as  in  conventional  MOSFETs . We 
discussed  this  unusual  behavior  in  Chapter  4 . 

The  LDD  (n_— drift)  region  is  used  to  block  high  voltages 
in  the  OFF-state  without  breaking  down,  but  it  also  affects 
the  current-voltage  characteristics  in  the  ON-state.  The 
simulated  ON-state  {\cs  >^Tr  ^ DS  > 0)  current  flowlines  shown 
in  Fig.  5.2  indicate  that  the  current  is  confined  near  the 
surface  of  the  LDD  region  due  to  the  LDD  — substrate  junction 
electric  field.  As  1 ps  is  increased,  the  junction  depletion 
region  widens,  and  hence  reduces  the  effective  cross-section 
of  current  flow  ( JFET  action)  . At  a critical  VDSf  for  a 
specific  VGs,  the  carriers  (electrons)  attain  their  maximum 
(saturated)  drift  velocity  near  the  drain,  and  begin  to 
accumulate  to  support  incremental  current . This  electron 
accumulation  counteracts  the  JFET  action  near  the  drain. 

The  npn  BJT  inherent  in  the  LDD  LDMOST  (see  Fig.  5.1(a))  is 
activated  when  a (transient)  reverse  current  flows  through  the 
LDMOST,  as  is  forced  by  some  power-switching  circuits,  e.g., 
switching  power  supplies  [Bal87] . The  BJT  causes  a reverse- 
recovery  delay  and  a troublesome  abrupt  current  change,  which 
we  discuss  later.  There  are  two  different  modes  of  operation 
of  the  BJT,  depending  on  the  bias  conditions  as  discussed 
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Fig. 


5.2  PISCES -simulated  current  flowlines  (Al=4  OxlO'6 
A//,m)  in  the  vicinity  of  the  drift  region  of  the 
LDD  LDMOST  (I'ps  - 16  V,  1D,  = 20  V,  . The  "unction 
depletion  region  is  indicated  by  the  dashed  curves . 
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in  Chapter  4.  To  study  the  two  modes,  we  have  done  PISCES 
simulations  of  the  reverse-recovery  characteristics.  For  the 
transient  simulations,  ]'dS  increases  from  — 5V  to  5V  through 
a load  resistor  in  a finite  rise  time,  with  T >/s-  = OV  (<  } j ) . 

PISCES -simulated  reverse-recovery  characteristics  is  shown 
in  Fig.  5.3.  To  indicate  the  activation  of  the  two  BJT  modes, 
we  have  separated  out  three  current  components:  Is(t), 

and  IB(t)  are  the  currents  flowing  into  the  n+-drain,  the 
n+— source,  and  the  p-body,  respectively.  For  0 < / < 45  ns, 
Is(t)  > 0,  implying  inverse-mode  operation  with  the  drain- 
substrate  junction  forward-biased.  For  45  ns  < / < 05  ns, 
Is(t)  < 0,  implying  forward-mode  operation  with  the  source- 
body  junction  forward-biased  (near  the  channel) . with  a 
negative  voltage  applied  to  the  drain  (source  and  body 
contacts  are  tied  to  ground)  at  / = 0 sec,  all  the  junction 
(P_P  ' P -n  , and  n~-n+)  potential  barriers  are  lowered  such 
that  conductivity  modulation  (high-injection)  is  prevalent 
in  the  p —substrate  and  n_-drift  regions.  A BJT  action, 
with  zero  collector-base  voltage,  supports  the  current: 
~ -Us  + Id)  = —{fll  + 1 )Ib-  Figure  5.4(a)  shows  PISCES- 
simulated  current  flowlines  at  t = 0 sec . It  reveals  the 
inverse-mode  BJT  action.  About  75%  of  the  total  current  (In) 
is  collector  current  (Is)  , electrons  flowing  in  the  n+-source. 
The  rest  of  the  total  current  is  base  current  (IB)  , holes 
flowing  in  from  the  body  contact  to  support  the  high-injection 
recombination  in  the  substrate/drift  region.  High  gain  ( j,)  of 


I(t)  [A///m] 
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Fig. 


5.3  PISCES -simulated  reverse-recovery  current  tran- 
^n?tSH°^the  LDM0ST-  ^e  drain-source  voltage, 
? L M v °^9  ^ <=  54  KS7— //m)  , changes  from  -5 

r arp  ,Whlth  a 0-1-/'«  rise  t^e.  ID(t),  7S(/),  and 
.J1  f e the  currents  flowing  into  the  n+-drain, 
the  n -source,  the  p-body,  respectively. 
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Fig.  5 


. 4 


PISCES-simulated  current  flowlines  (a)  at  / = 0 s 
in  the  reverse-recovery  transients  in  Fia  5 3 

of  the  inverse-mode  Bji 
(Al-7.4xl0  A///m)  and  (b)  at  /=  50  ns,  showino 

A/^mtCtlVati°n  °f  thS  forward-mode  BJT  (Al=8.1xl0~^ 
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this  inverse-mode  transistor  decreases  the  reverse-recovery 
delay  because  there  is  less  charge  stored  corresponding  to  a 
specific  reverse  current.  Optimal  design  of  the  device  would 
hence  require  a high  gain  for  the  inverse-mode  BJT,  but  a 
low  gain  for,  or  suppression  of  the  forward-mode  BJT  (using 
the  deep  p+-body  diffusion),  which  we  now  discuss. 

The  direction  of  /y(/)  changes  during  the  reverse  recovery 
(at  t ~ 7 ns  in  Fig.  5.3)  due  to  the  increased  current  ( <IQp/dt ) 
discharging  the  substrate  (and  body).  This  (dis) charging 
current  flowing  through  the  lateral  body  resistance  B,,b 
produces  a voltage  drop,  hence  forward-biasing  the  source- 
body  junction  near  the  channel  and  inducing  electron  injection 
into  the  p-body  and  p~-substrate . The  forward-mode  BJT 
is  thereby  activated.  Figure  5.4(b)  shows  PISCES-simulated 
current  flowlines  at  / = 50  ns  in  the  reverse  recovery  showing 
this  BJT  action.  About  35%  of  the  total  current  (ID)  is  due 
to  the  forward-mode  BJT  (/,)  at  this  time.  This  BJT  mode 
is  undesirable  for  several  reasons.  First,  it  increases  the 
reverse-recovery  delay  because  of  added  charge  stored  in  the 
p-— substrate  region.  Second,  it  induces  an  abrupt  change  in 
current  ID(t)  evident  in  Fig.  5.3  because  of  the  fast  switching 
of  /,(/).  This  phenomenon,  called  "snap  recovery"  [Bal87],  can 
generate  a very  high  voltage  drop  across  stray  inductance, 
which  can  induce  dV/dt  turn-on.  Finally  the  turn-on  of  the 
forward-mode  BJT  can  result  in  saturation  of  the  npn  BJT  and 
possibly  destructive  breakdown  of  the 


LDMOST  [Bal87] . 
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Since  HVICs  require  high  voltage/current  operation  of  the 
LDMOST,  impact  ionization  in  the  LDD  region  can  be  important. 
We  briefly  discuss  this  phenomenon.  When  the  LDMOST  operates 
at  high  drain  bias,  carriers  (electrons)  moving  through  the 
high-field  LDD  region  acquire  high  enough  kinetic  energy 
to  generate  electron-hole  pairs  via  impact  ionization.  The 
generated  electrons  add  to  the  drain  current  and  generated 
holes  are  swept  into  the  p-body  region,  increasing  the  body 
potential  and  thus  activating  the  forward-mode  npn  B JT . Our 
PISCES  simulations  reveal  that  the  impact  ionization  occurs 
in  the  LDD  region,  and  not  in  the  channel  region.  The 
impact  ionization  occurs  first  near  the  drain  where  the 
electric  field  is  highest  and  spreads  toward  the  channel  as 
Vds  is  increased.  Then  both  the  LDD  and  substrate  regions 
are  conductivity-modulated  with  electron-hole  plasma,  which 
results  in  lower  voltage  drop  across  the  device.  This  is  the 
so  called  "latch-back"  effect  [Huc82] . 

5.2.2  RESURF  LDMOST 

The  RESURF  LDMOST,  shown  in  Fig.  5.1(b),  is  similar  in 
structure  to  that  of  the  LDD  LDMOST  with  the  exception  of 
the  drift-region  design.  The  whole  n'-epi  layer  functions 
as  a drift  region,  and  its  thickness  and  doping  density  are 
controlled  to  produce  low  surface  electric  field  to  provide 
high  breakdown  voltage  [App79] . We  discuss  the  n'-drift  region 
in  the  normal  mode  of  operation  and  the  npn  B JT  in  the  reverse 
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mode  of  operation.  The  double-diffused  channel  is  as  same 
as  that  of  the  LDD  LDMOST . 

Under  normal  operation,  the  reverse-biased  epi-substrate 
junction  has  a depletion  region  which  effectively  reduces 
the  cross-sectional  area  of  electron  current  flow,  as  in 
the  JFET.  At  high  currents,  excess  electrons  accumulate  in 
the  drift  region  due  to  drift-velocity  saturation  (which  was 
ignored  [Cla86]>.  The  excess  electrons  produce  space-charge- 
limited  (SCL)  current  flow  and  modulate  the  epi-substrate 
junction  barrier,  thus  compensating  for  the  widening  of  the 
depletion  region.  This  quasi-saturation  is  similar  to  that 
in  the  LDD  LDMOST,  but  a main  difference  between  the  two 
versions  of  the  LDMOST  is  the  magnitude  of  on-resistance; 
the  RESURF  LDMOST  shows  smaller  on-resistance  due  to  the 
thicker  drift  region.  Figure  5.5(a)  shows  PISCES-simulated 
current  flowlines  in  the  RESURF  device  for  normal  ON-state 
operation.  The  gate  voltage  is  biased  very  high  to  avoid 
channel-current  saturation.  The  current  flow  in  the  drift 
region  near  the  channel  is  controlled  by  the  reverse-biased 
epi-substrate  junction  (JFET  action) . Near  the  n*^irain,  the 
cross-sectional  area  of  electron  current  is  smaller  due  to 
larger  reverse  bias  across  the  epi-substrate  junction,  and 

" >'V°  dUe  t0  the  ^rift-velocity  saturation.  This  conductivity 
modulation  compensates  the  JFET  action. 

The  inherent  npn  BJT  in  the  RESURF  LDMOST  is  slightly 
different  from  that  in  the  LDD  LDMOST  (see  Fig.  5.1).  when 
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Fig.  5 


IT 


(a) 


(b) 


PISCES -simulated  current  flowlines  (a)  in  the 
normal-mode  operation  of  the  RESURF  LDMOST 
6 ' ],)S=40Vf  and  ^I=7.8xl0-6  A/,,m)  and  (b) 

n-  T7reyerse"mode  operation  of  the  RESURF  LDMOST 
(Ic.’S-OV,  \ ds=~  0.8V,  and  Al=1.8xlQ-6  A//>m)  . 
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a negative  voltage  is  applied  to  the  drain  (source,  body, 
and  substrate  contacts  are  shorted  to  ground) , inverse-  and 
forward-mode  BJT  action  can  be  triggered,  as  in  the  LDD 
LDMOST.  But  the  RESURF  LDMOST  has  a vertical  pin  diode,  which 
is  like  an  extrinsic  base.  Furthermore,  high-injection  occurs 
in  the  p~— substrate  where  significant  recombination  occurs. 
Figure  5.5(b)  shows  PISCES-simulated  current  flowlines  in 
the  reverse-mode  operation  (VGS  = 0 V,  VDS  = -0.8  V).  There 
are  three  important  current  components  in  this  mode:  1)  the 
collector  current  flowing  in  the  n+— source  (/s-  ~ 0.40  Id);  2) 
the  recombination  current  flowing  in  the  p-body  (If),  ~ 0.15 
Ip) / 3)  the  recombination  current  flowing  in  the  p~— substrate 
{I lie  — 0.45  Id)  • The  ratio  of  two  base  recombination  currents 
(Ilii/Ilie)  varies  depending  on  the  device  structure.  For 
example,  a thick  substrate  reduces  /gf. 

The  reverse-recovery  characteristic  of  the  RESURF  LDMOST 
is  similar  to  that  of  the  LDD  LDMOST,  but  it  exhibits  a 
longer  delay  for  the  same  carrier  lifetime.  This  is  due  to 
the  additional  charge  stored  in  the  vertical  pin  diode,  i.e., 
in  both  the  epi  and  substrate  regions . 

5 . 3 LDMOST  Model 

Based  on  physical  insights  gained  from  the  simulations 
discussed  in  Section  5.2,  we  model  the  LDMOST  regionally 
partitioning  it  into  three  main  components:  the  channel,  the 
drift  region,  and  the  inherent  BJT.  We  describe 


a physical 
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composite  model  for  the  LDD  n-channel  LDMOST,  which  is 
implemented  in  SPICE  and  hence  has  utility  for  HVIC  CAD. 
The  modeling  methodology  is  also  applicable  to  the  RESURF 
LDMOST . 


5.3.1  LDMOST  Channel 

The  LDD  LDMOST  channel  is  modeled  based  on  the  VDMOST 
channel  model  in  Chapter  4 , but  expanded  to  account  for  the 
extended  channel  over  the  substrate  as  illustrated  in  Fig. 
5.1(a).  We  thus  approximate  the  channel  doping  density  as 

i\A  % NAi)exp( -ifx/Li ) 0 < r < Lx 

= Nas  Li  < r < L\  + L>  (5.1) 

where  L\  and  Z2  are  the  double-diffused  and  substrate  channel 
lengths.  We  use  a piecewise-continuous  electron  velocity 
model : 


PneffEr 

1 + [PneffEr/2vs„i) 
~ l' sat 


for  v < 
othtnvisf 


(5.2) 


where  Ex  - |c/T  /dx\  is  the  longitudinal  electric  field, 
(as  lO'cm/s)  is  the  saturation  velocity,  /inrff  = 
Pao/[l  + 0(Vas-VT)]  where  is  the  low-field  mobility,  and 

0 is  the  fitting  parameter  which  defines  the  transverse-field 
dependence . 
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Following  the  methodology  used  in  Chapter  4,  we  describe 
the  strong-inversion  channel  current  by  integrating  Irh  = 
-N'V’Qn  from  r = 0 to  Ij  + L>: 

^ zUntf  f 


2(1]  +L2)  [1  + (untff/ 2rsa(  (Li  + L,))  T'r/,] 
Qt(L\)  — Qj  ( 0 ) ( Qi[L{+L,)-Qi(Lx) 

Cox  4*  Crfi  C or  "I-  C|/,i 

where  Vck  is  the  total  channel  voltage  drop,  and 
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(5.4) 


(5.5) 


are  the  body  depletion  "capacitances"  for  the  double-dif fused 
and  substrate  channels,  respectively.  In  (5.4),  Tj  is  the 
double-diffused  channel  voltage  drop,  which  is  defined  by 
equating  the  following  two  equations  for  channel  current 
obtained  from  integrations  along  the  respective  channels: 

W zlhie  f f 
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Qi(Li)-QUO) 

L\  [l  + | ) ( j]  2 ( Cor  + C ) 
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Li  [l  4-  (//„f / / /2vs,uL  >)  0 ch  — ( | )]  2 (C'0.F  + C',/2 

In  this  formalism, 


(5.6) 

(5.7) 


Qn  (0)  — ~C0- r ( 1 QS  — 1 7 


Qn  ( ■/- 1 ) — —Cor 


1 f/s  - 1 r - 1 1 ( 1 4-  77 


C 


C or 


(5.8) 

(5.9) 


Qn  (-£1  4-  L >)  = —C0 


rGS-rr-r,(i  + £kL)_(,-4-r|,/i  + ^ 


^ o.r 


Co, 


(5.10) 


107 


and  the  threshold  voltage  VT  is  defined  by  the  doping  density 
at  the  source  as  in  Chapter  4 : 


\/4</f  <Y.l0  \<Pd\ 


n r 


(5.11) 


The  channel -current  saturation  is  caused  by  electron 
velocity  saturation  near  x=0,  which  implies  a simple 
expression  for  the  saturated  current  as  in  Chapter  4: 

^ch(sat)  W n ( •('  = 0 ) 

= i vsnl  Cof  ( \ };s  — ] ) . (5.12) 

Equating  (5.12)  and  (5.3)  gives  Vch{sut)  for  a given  Vns,  above 
which  the  channel  current  saturates. 

5.3.2  LDP  Region 

The  LDD  drift  region,  in  the  general  case,  is  further 
divided  into  two  sections:  section  A,  near  the  substrate 

channel,  is  the  f ield-dependent-velocity  region  where  the 
JFET  action  takes  place;  section  B,  near  the  n+-drain,  is 
the  (high-field)  saturated-velocity  region.  At  low  currents, 
the  LDD  region  exhibits  only  JFET  action,  viz.,  only  section 
A.  At  high  currents,  section  B materializes.  Let  us  consider 
two  regions  of  operation  defined  by  ILDD  relative  to  ILDD[sat)t 

the  drain  current  at  which  the  electron  velocity  saturates 
at  the  drain  in  the  LDD  region  (,•  = £,,  in  Fig.  5>6)> 
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P 


Fig.  5 


.6  Magnified  view  of  the  LDD  region  showing  the 
Gaussian  volume  (thick  solid  lines) . 
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I LDD  < lu)P(Mt) 

In  this  case,  Er  is  not  high  to  cause  velocity  saturation. 
Hence  the  carrier  density  does  not  exceed  the  doping  density 
Nfj,  and  the  LDD  region  (only  section  A exists)  is  described 
by  conventional  JFET  theory  [Sze81] . The  drain  current 
is  expressed,  using  a continuous  electron  velocity-field 
function,  as 


1 1 ,DD  = 11%  [llj 


yd  ( •<’ )]  qXo 


I1  ml 'Er 

1 + Er/Er 


(5.13) 


where  yj  is  the  metallurgical  LDD  junction  depth,  y,i(.r)  is  the 
extent  of  the  junction  depletion  region  into  the  LDD  at  .»•, 
/'»(/  is  the  low-field  electron  mobility,  and  Er  = r.„, ////„, / . 

Integrating  (5.13),  we  get  the  relationship  between  the 
drain  current  and  the  LDD  region  voltage  drop  Vldd : 


I LDD  = 


(i + fe) 


2e  jS 


AS 


<1^  u ( Nas  + , 


1 LDD  + 1 bi ) ,/2  - 1 ’b/' 


(5.14) 


which,  for  a given  IujDr  is  solved  analytically  for  Yldd 
Note  from  (5.13)  that 


hDD(sai)  = IV,  [f)j  - !J,I  (L,i)]qXDcs,lt  (5.15) 

where  )/,/(£,/)  is  given  by  the  depletion  approximation  and 
^ LDU{$at)r  which  follows  from  equating  (5.15)  to  (5.14). 
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I LDP  > hpp(  sal) 

In  this  current  range,  the  LDD  region  is  subdivided  into  two 
sections,  A and  B,  since  the  electron  velocity  saturates  near 
the  drain.  In  section  A,  electrons  move  by  the  field-dependent 
velocity,  but  in  section  B by  the  saturation  velocity.  At 
the  boundary  r = L.\  between  these  two  sections,  the  drain 
current  is,  from  (5.13),  ILDD  = IT’,  [yj  - yj  (L  .4)]  qN/jr,,,,  where 
ij,i(La)  depends  on  the  voltage  VA  at  v = L \ through  the  depletion 
approximation.  This  current  expression  yields 

h.pp  V (1^p  ( Has  + ftp ) 


'a  = [yj- 


- r. 


bt 


(5.16) 


H pl'sat ) 2ei\ 

where  i-3  the  sxabstrate  doping  density  used  in  (5.1). 

With  (5.16),  the  JFET  theory  applied  to  section  A defines 
the  section  A length,  L A . 

In  section  B,  we  apply  Gauss's  theorem  to  a volume  enclosed 
by  planes  of  unit  width  (see  thick  solid  lines  of  Fig.  5.6)  : 

fh  Hj  r !lj 

6 I Ejd-V  + f I ( —Ea  ) <hj  + f I Efdy  = q J J (Xu  — n ) dyd.c 


(5.17) 


0 0 [,  A 0 

where  the  integral  component  along  the  top  surface  of  the 
LDD  is  neglected  due  to  the  small  transverse  electric 
field  there.  To  further  simplify  the  analysis,  we 

assume  that  the  LDD-substrate  junction  electric  field, 
Ej  — 1/2  [E,  (.<•  = L A)  + Ej  ( ,c  = L,/)],  is  an  average  value  independent 
of  r . The  depletion  approximation  is  used  to  evaluate  this 
average.  Also  E,  1 and  Er  are  assumed  to  be  independent  of 
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position  y.  All  of  these  assumptions  concerning  the  electric 
field  are  consistent  with  PISCES  simulations. 

Differentiating  both  sides  of  (5.17)  with  respect  to  x 
gives 


T” 

fEj  - njj  = (ilXuyj  + Qn 


= <1X'DIJ)  - 


hup 

H : I'.itit 


(5.18) 


This  equation  is  integrated  twice  to  give  the  section  B 
voltage  drop 


' 7j  = 


( L,i  - L \ )~ 
'2eyj]Vzrslll 


(Ildd 


IT  f«yi\  pijjVsai  + f E j 1 1 ) -f  ( L,i  — ) E,  . (5.19) 


Now  the  complete  LDD-region  voltage  drop  is 

XLDD  = XA  + XB  < (5.20) 

where,  for  a given  IldDi  1 .1  and  1#  are  characterized  by 
(5.16)  and  (5.19).  The  equations  constitute  the  LDD-region 
model.  Because  Ej  in  (5.19)  depends  on  YldDi  equations 

require  iterative  numerical  solution,  which  is  done  within 
the  model  routine  implemented  in  SPICE  (see  Section  5.4). 

5.3.3  npn  BJT 

The  inherent  npn  BJT  of  the  LDD  LDMOST  is  modeled  based  on 
physical  insights  from  Section  5.2  and  on  our  VDMOST  BJT  model 
in  Chapter  4 . Because  most  of  the  BJT  current  flows  through 
the  p —substrate  as  shown  in  Fig.  5.4(a),  the  one-dimensional 


112 


n+— p-p~— n+  structure,  shown  in  Fig.  5.7,  is  representative 
of  the  basic  device  operation.  As  discussed  in  Section  5.2, 
the  BJT  can  possibly  operate  in  the  forward  mode,  where 
the  source-body  (n+— p)  junction  is  forward-biased,  or  in  the 
inverse  mode,  where  the  drain- substrate  (n+— p~)  junction  is 
forward-biased.  In  the  latter  case,  carriers  (electrons)  are 
hardly  injected  through  the  LDD-substrate  (n~— p_)  junction 
due  to  the  voltage  drop  (debiasing)  across  the  LDD  region. 

The  integral  charge-control  relation  (ICCR)  is  applied 
only  to  the  p-body  base  region  to  obtain  the  base  transport 
current.  Actually  there  is  a composite  base  region  for  the  npn 
BJT,  comprising  the  p-body  base  and  the  p-— substrate  base.  But 
we  focus  on  the  p-body  base  where  Jp  ~ 0;  the  predominant  hole 
current  flows  vertically  through  the  body-substrate  junction 
under  the  body  contact,  removed  from  the  intrinsic  p-body 
base,  as  shown  in  Fig.  5.4(a).  The  ICCR  in  the  quasi-neutral 
p-body  base  region,  for  low  injection,  gives 


gD"  [»  (rFp)  p (rEp)  - n (■)•/?,,)  p (.vr,,)] 
tb>. 

f Na  ( .r ) </.r 

J'Ep 


where  for  quasi-equilibrium  in  the  junctions. 


(5.21) 


”(•'■£>)/>  (•'•£>)  = "hrp(VuEn/VT)  (5.22) 

and 

n (*BP)p (*bp)  - p1  ( n I,,-) 

— ^ as (XP  (2Vjpp-/Vt) 


(5.23) 
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Fig.  5.7  The  one-dimensional  n+-p-p--n+  transistor  of  the 
LDD  LDMOST. 
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The  junction  voltages  and  boundaries  in  (5 . 21) - (5 . 23)  are 
defined  in  Fig.  5.7.  High  injection  was  assumed  in  the 
p~— substrate  (p~n>NAS)  to  derive  (5.23).  Combining  these 
equations  yields  the  net  p-base  transport  current, 


r A'<iD„n r 

IcT  ~ 7STZ BcJ' r)  ~ ' <-rp(2Vj„-/ IV) 


r r T 13  p 

f A t ( r ) d.r  I jyl  ( ) flx 

r*p  rLp  ' 

= IcCOf-rpCl'BEv/^T)  ~ I 'rP  (21  JpV-  /l  V)  • 


(5.24) 


where  superscripts  F and  I designate  different  parameters 
corresponding  to  the  forward  and  inverse  BJT  modes . 

The  BJT  base  (recombination)  current  in  the  forward  mode 
is  typically  negligible  relative  to  the  hole  charging  current 
UlQ,,/(lt)  that  triggers  the  forward-mode  BJT  in  transient 
operation  of  the  LDMOST . The  current  is  thus  not  modeled. 
In  the  inverse  mode  however,  the  recombination  in  the 
p -substrate  defines  a significant  component  of  base  current. 
To  characterize  this  recombination,  we  solve  the  ambipolar 
transport  equation  in  the  substrate,  where  we  assume  high 
injection.  The  solution  [Fos88]  is 

„ ( r)  _ " (*Bp-)shth  (~f,~r)  + » {*CP-)*inh 


where  the  boundary  carrier  densities  are  defined  by  the 
junction  voltages : 


" ( vBp~)  - Pi-van-)  = y,\sejrp{Vfpp-/Vr)  . 

n (rO>~)  = ("i/-x  \s)  f.rp(VJp-n  + /Yr)  . 


(5.27) 
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The  boundary  current  conditions  are 


h U‘Bi>~)  — ~h'T 


(5.28) 


and 


IP  (scp-)  - 0 


(5.29) 


where  we  assume  negligible  recombination  in  the  n+— drain 
region.  Equations  (5 . 26)  - (5 . 29)  define  n (. vBp~)  and  » (.rC/,-)  for 
given  \ bEh  and  Yjpp-  [McD85]  . Now  the  recombination  current 
for  inverse-mode  operation  is  given  by 


where  77/  is  the  (high-injection)  carrier  lifetime  in  the 
p~— substrate . 

The  substrate-region  voltage  drop  is  significant  for 
inverse-mode  BJT  operation.  It  is  characterized  by  integrating 
the  electric  field  in  the  substrate,  derived  from  the 
ambipolar- transport  analysis  [McD85] , and  by  adding  r,  + ; 


Note  here  that  the  two  bases  (body  and  substrate)  are 
linked  by  (5.24)  and  (5.28);  Icr  is  controlled  by  Yjpp-  and 
Ibi  is  a function  of  ICT.  Thus  for  specified  junction  voltages 


(5.30) 


(5.31) 


^BEh  and  YJpp-r  the  above  formalism  analytically  describes  the 
pertinent  BJT  characteristics,  i.e.,  ICT,  IBjf  and  Y,lth. 


116 


5 . 4 SPICE  Simulations  And  Discussion 
The  physical,  semi-numerical  composite  model  for  the  LDD 
LDMOST  described  in  Section  5.3  is  implemented  in  SPICE. 
Figure  5.8  shows  the  network  representation  of  the  model, 
based  on  integration  of  the  regional  analysis  with  I,h  = Ildi )• 
This  is  not  a simple  equivalent-circuit,  but  represents  the 
physical  charge-based  model.  All  the  elements  in  Fig.  5.8 
are  implemented  in  SPICE  via  user-defined  controlled  sources 
(UDCSs)  which  access  the  model  routine  that  numerically  solves 
the  implicit  system  of  equations  developed  in  Section  5.3. 
UDCSs  are  FORTRAN  subroutines  in  SLICE  [Har84],  an  enhanced 
version  of  SPICE2 . 


The  sources  in  the  model  network  in  Fig.  5.8  are  voltage- 
controlled  current/voltage  sources  defined  in  terms  of  node 
voltage-differences  1 qs  > ^ rh  > and  1 z?ir»  • The  voltage  source  Y,infi 
represents  ] ldd  in  normal-mode  operation  of  the  LDMOST  and 
Yaub  in  reverse-mode  operation.  The  transient  characteristics 
of  the  device  are  simulated  by  the  charging— current  UDCSs, 


']Qi  = v 0Q,'lYj 

( It  DYj  (It 


(5.32) 


where  i — P,  ./,  and  j — GS,  eh,  DEn.  This  charge  dynamics 
is  based  on  the  quasi-static  approximation  with  the  charges 
derived  from  the  analysis  in  Section  5.3: 


Q, 


w 


n ( .r  ) <l,v  . 


nr 


(5.33) 
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Gate 


Fig.  5.8  Network  representation  of  the  LDD  LDMOST  model. 
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with  n (~  p)  given  by  (5.25),  is  the  hole  charge  stored  in  the 
p~— substrate  for  inverse-mode  operation  of  the  inherent  BJT; 


with  \ tqs  = \\h + 1 LDDr  is  the  LDD-substrate  junction  depletion 
charge  for  normal-mode  operation  of  the  LDMOST . (The  excess 
electron  charge  in  the  LDD  due  to  velocity  saturation  is 
negligible  as  discussed  in  Chapter  4 . ) Also  included  in 
the  model  is  C'qsi  the  gate-source  capacitance,  the  most 
significant  component  of  which  is  that  due  to  the  gate 
electrode  overlapping  the  n+— source.  The  charging  currents  are 
calculated  numerically  using  finite-difference  approximations 
for  the  voltage-derivatives . These  approximations  are  also 
used  for  the  transcapacitances  and  transconductances  needed 
for  circuit  nodal  analysis. 

In  the  nodal  analysis,  SPICE  calls  the  UDCSs  iteratively, 
which  in  turn  access  the  LDMOST  model  routine  flowcharted 


(relative  to  Ildd\*<,i))  paths  exist,  both  yielding  the  LDMOST 
solution  with  Vds  = Vch  + ^ LDD  • In  the  reverse  mode  (T'y,  < ())  , 
the  BJT  analysis  is  activated  with  1 = — T *./, , yielding  the 
BJT  solution  with  VDS  = YJpp.  + V9Kb. 

Test  structures  of  LDD  LDMOST s were  fabricated  using 
lateral-charge-control  technology  at  Siliconix,  Inc.  [Wil87] . 


(5.34) 


in  Fig.  5.9.  The  model  routine,  given  V(;Sf  Vch,  and  YIJEnl 
first  determines  the  mode  of  operation  by  checking  the  sign 
of  \ </, . In  the  normal  mode  (],./,  > 0)  , low-  and  high-current 


119 


Fig.  5.9 


The  LDD  LDMOST  model  flowchart . 
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Figure  5.10  shows  measured  and  simulated  Ins  - ]'l)s 
characteristics  for  the  LDD  LDMOST . For  the  simulations, 
model  parameters  were  crudely  extracted  to  fit  measured  data 
without  any  optimization.  Other  parameters  like  geometry  and 
doping  density  are  given  technologically.  Quasi-saturation, 
in  which  IDS  is  independent  of  VGS  and  the  output  conductance 
is  high,  is  predicted  in  simulations  for  l'c;s  = 8,  10,  12 
V,  m agreement  with  the  measurements.  The  abrupt  current 
transition  from  the  triode  region  to  the  saturation  region  for 
the  lower- Vos  simulations  is  due  to  the  piecewise-continuous 
electron  velocity  model  in  (5.2) 

In  Fig.  5.11,  we  show  measured  and  model-simulated  reverse- 
recovery  characteristics  of  the  LDD  LDMOST.  When  t < 100  us, 
the  device,  in  reverse-mode  operation  (VDS  < ()  V)  where  the 
inverse-mode  BJT  is  activated,  stores  a high  concentration  of 
holes  and  electrons  in  the  p--substrate  and  n'-drift  regions . 
At  t = 100  ns,  a voltage  pulse  (-5  V to  +4  V with  100-ns 
rise  time)  is  applied  to  the  drain-source  terminals  through 
a resistor  (152  ft)  to  switch  the  device  from  the  reverse- 
ON-state  to  the  normal-OFF-state . Then  the  stored  holes  and 
electrons  in  the  substrate/drift  regions  begin  to  diminish 
predominantly  by  charging-current  JQJJt,  and  by  recombination 
current  IBl.  The  time  delay  from  initiation  of  switching  until 


the  device  recovers  to  a specified  OFF-state  is  defined 
the  "reverse-recovery  time".  The  reverse-recovery  time 
controlled  by  a)  the  rise  time  of  the  switching  pulse; 


as 

is 

b) 


IDS  [A] 
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Fig.  5.10  Measured  (small  rectangles)  and  SPICE-simulated 
(curves)  current-voltage  characteristics  of  the 
LDD  LDMOST  {Ld  = 42  //m)  . 


IDS  [mA] 
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5.11  Reverse-recovery  characteristics  for  the  LDD 
LDMOST  (L,i  = 42  //m)  : (a)  measured  (1  division 
on  the  x-axis  = 100  ns  and  1 division  on  the  y- 
axis  = 13.2  mA)  ; (b)  SPICE-simulated . The  drain- 

source  voltage,  applied  through  Rj  (=  152  Q) 
changes  from  -5  V to  +4  V with  a 100-ns  rise  timi 
starting  at  t = 100  ns. 


123 


TH ' and  c)  RpB-  Fast  switching  rise-time,  short  77/,  and 
small  RpB  reduce  the  reverse-recovery  time.  In  Fig.  5.11  the 
forward-mode  BJT  is  not  triggered  because  the  induced  ohmic 
voltage  drop  across  RpB  is  not  sufficient  to  forward-bias 
the  source-body  junction.  Note  that  the  measurement  shows 
a long  recovery  tail  due  to  stray  capacitances  on  the  probe 
station.  However  the  storage-time  simulation  corresponds 
well  with  the  measurement. 

Figure  5.12  shows  our  test  circuit  for  transient 
simulations.  This  circuit,  which  is  a basic  building  block  for 
a variety  of  DMO ST  applications,  e.g.,  two  quadrant  chopper 
circuit,  tests  both  normal—  and  reverse-mode  operations, 
including  transient  reverse  recovery.  Figure  5.13  shows 
results  of  transient  simulations.  Initially  (t  = 0)  , VCt ■S>VT; 
so  LDMOST1  is  OFF  and  LDMOST2  is  ON.  As  VGS  goes  to  zero  ( 
t = 0.15  /is)  r LDMOST2  is  turned  off  and  IDSI  decreases. 
Consequently,  the  current  through  Ll,  ILr  decreases  and 
induces  a voltage  drop  across  Ll . Then  the  voltage  at  node  P, 
YP,  is  greater  than  the  supply  voltage  VBr  and  this  voltage 
difference  turns  on  the  npn  BJT  (so  called  the  integral  diode) 
of  LDM0ST1.  when  VGS  increases  again  ( / = 0.25  //s)  , the 
npn  BJT  of  LDMOST1  should  be  turned  off,  but  undergoes  the 
reverse  recovery,  in  a mean  time  LDMOST2  is  turned  on,  but 
Iosi  = lDSi(ON)+It{R\,  where  IRRx  is  the  reverse-recovery  current 
of  LDMOST1 , because  IL  ( = I 0S2  ~ IDs\ ) changes  very  little  due 
to  Ll.  This  peak  current  occurs  at  an  instant  when  VDS2  is 
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Fig.  5.12  Test  circuit  for  the  LDD  LDMOST  transients.  Ra 
50  il,  Ll  = 5 mH,  Rl  = 250  12,  and  VD  = 20  V.  ' 
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Fig.  5.13  SPICE  simulation  results  for  the  circuit  in  Fig. 
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also  high.  This  results  in  high  power  consumption  which  can 
induce  second  breakdown  by  the  impact-ionization  discussed 
in  Section  5.3.  //?/?]  is  higher  when  the  turn-on  time  of  \r;s  is 
shorter.  In  this  case  the  second  breakdown  is  more  probable. 
The  simulations  shown  in  Figs.  5.11(b)  and  5.13,  which  reveal 
effects  of  the  inherent  BJT  in  LDMOST  circuits,  demonstrate 
a capability  of  mixed-mode  device/circuit  simulation.  This 
feature  of  our  physical  model,  which  is  not  afforded  by 
other  equivalent- (sub) circuit  models,  can  be  quite  useful  in 
computer-aided  optimal  device/circuit  design  of  HVICs . 

5 . 5 Summary 

New  physical  insights  with  regard  to  two  structures  of 
LDMOSTs  (LDD  and  RESURF  LDMOSTs)  have  been  obtained  by  using 
the  two-dimensional  device  simulator  PISCES.  These  insights 
were  used  to  develop  a physical,  semi -numerical  composite 
model  for  the  LDD  LDMOST.  The  modeling  methodology  is  also 
applicable  to  the  RESURF  LDMOST.  The  charge-based  model  is 
implemented  in  SPICE  via  UDCSs  and  supported  by  measurements 
of  test  devices.  Transient  simulations  for  both  normal- 
and  reverse-mode  operations  were  demonstrated  to  show  the 
superiority  of  our  physical  model  over  empirical  equivalent- 
(sub) circuit  models  for  HVIC  CAD. 

The  developed  physical  SPICE  model  is  useful  for  computer- 
aided  optimal  design  of  LDMOSTs  as  well  as  of  HVICs.  Device 
and  structural  model  parameters,  e.g.,  th  , B,,Br  and  L,lf  are 
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adjusted  to  derive  optimal  device/circuit  performance.  For 
example,  to  satisfy  the  specified  reverse-recovery  time,  deep 
p+_ diffusion  into  the  p-body  is  employed  to  reduce  Rvh,  or  the 
electron  irradiation  [Bal87]  can  be  used  to  reduce  77/ . Then 
subsequent  device/circuit  simulations  can  be  done  directly 
to  predict  effects  of  the  adjustments  on  performance. 


CHAPTER  6 

SUMMARY  AND  SUGGESTIONS  FOR  FUTURE  WORK 


In  this  dissertation,  new  models  for  MOS-controlled  HV 
devices  (IGBTs  and  DMOSTs)  have  been  developed  and  shown 
to  be  useful  for  optimal  device/circuit  CAD  of  HVICs . The 
unique  features  of  HV  devices  such  as  latch-up  in  the  IGBT 
and  quasi-saturation  in  the  DMOST  were  studied  extensively 
using  the  two-dimensional  numerical  device  simulator  PISCES. 
The  physical  insights  gained  from  the  numerical  simulations 
and  experimental  measurements  of  test  devices  were  used  to 
develop  our  charge-based  HV-^levice  models.  The  models  were 
implemented  in  SPICE2  via  FORTRAN  subroutines  (UDCSs) , by 
which  the  model  equations  are  solved  semi-numerically  within 
the  nodal  analysis  framework  of  SPICE.  Our  new  models  were 
verified  by  measurements  of  specially  designed  test  devices, 
using  simulations  with  model  parameters  evaluated  from  device 
structural  information  and/or  crudely  extracted  without  any 
optimization  to  fit  the  measured  data.  Furthermore  circuit 
simulations  of  HVICs  were  done  to  demonstrate  the  superiority 

of  our  model  over  equivalent-circuit  models  with  regard  to 
technology  CAD . 

Based  on  the  research  discussed  herein,  we  suggest  the 
following  topics  for  future  research  consideration. 
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First,  we  recommend  a study  of  optimum  cell  design 
of  HV  devices  based  on  our  modeling  methodology.  The 
HV-device  characteristics,  e.g.,  ON-resistance  of  the 
DMOST  and  latch-up  of  the  IGBT,  are  strongly  dependent 
on  the  (three-dimensional)  cell  design:  triangular, 
square,  rectangular,  hexagonal,  circular,  or  parallel-line 
structures.  Three-dimensional  device  simulations  may  be 
required  to  study  the  cell  design.  Our  models  effectively 
account  for  the  device  layout  of  HV  devices;  the  device  width 
is  taken  as  the  sum  of  sub-cell  widths. 

Second,  we  recommend  PISCES  modifications  to  account  for 
general  terminal  boundary  conditions  needed  in  HV-device 
simulations.  The  HV  devices  with  circuit  or  stray  inductances 
loading  the  terminals  cannot  be  simulated  by  PISCES. 

Third,  we  recommend  the  extension  of  our  DMOS  channel  model 
for  all  regions  of  operation  including  the  reverse  channel- 
drive  mode.  During  the  transients  of  HVIC  simulations,  DMOSTs 
can  operate  in  this  reverse  mode . 

Fourth,  we  recommend  the  development  of  a latch-back  model 
for  the  DMOSTs.  As  the  DMOSTs  in  HVICs  are  subjected  to  unusual 
high  current /voltage  operating  regimes  due  to  inductive  loads, 
latch-back  can  be  an  important  phenomenon. 

Fifth,  we  recommend  the  implementation  of  our  models 
directly  into  SPICE  source  code.  The  UDCS  implementation 
which  we  have  used  is  CPU-intensive  and  is  intended  for 
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preliminary  SPICE  model  verification  and  not  for  practical 
CAD  application  of  the  models.  The  direct  SPICE  implementation 
requires  the  model  routine  to  evaluate  all  the  element  values 
from  the  device  terminal  voltages,  instead  of  the  internal 
node-voltage  differences  used  in  the  UDCS  implementation.  The 
required  modification  in  the  model-routine  formalism  however 
is  straightforward. 

Sixth,  we  recommend  overall  model  refinements,  for  example 
to  account  for  parasitic  effects  in  actual  multi-cell 
devices  and  HVICs . In  particular  the  gate-drain  charging 
current  <IQgd/( ^ the  VDMOST  model  (Chapter  4)  is  possibly 
inadequately  characterized.  Under  certain  conditions,  COD  = 

^Qgd/OYcd  can  show  an  anomalous  stronger  1 £>.$•— ' dependence  than 
that  described  by  our  Q(; D ( ( 'Ds ) . 


APPENDIX  A 

BASE  CHARGE  PARTITIONING  IN  WIDE-BASE  BJT' S 


In  this  appendix  we  describe  the  base  charge  partitioning 
of  the  wide-base  p+— n~-p  BJT  structure  in  the  IGBTs  discussed 


non-quasi-static  behavior  in  bipolar  transistors  was  defined 
[Foss86]  for  the  simple  case  of  low  injection  in  a 
uniformly  doped  quasi-neutral  base  region  with  negligible 
recombination.  These  conditions  imply  no  electric  field  and 
hence  no  drift  current.  For  the  charge  partitioning  in  the 
IGBT,  we  must  account  for  high  injection,  which  implies  an 
electric  field  and  a drift  component  of  current,  and  for 
recombination  in  the  wide  base  of  the  constituent  lateral 
pnp  transistor. 

Following  the  derivation  [Foss86]  but  accounting  for 
these  IGBT  properties,  we  integrate  the  one-dimensional  hole 
continuity  equation  over  the  quasi-neutral  base  width  TT  to 


in  Chapter  3 . A base  charge  partitioning  to  best  represent 


derive 


II 


(A.l) 


o 


131 


132 


where  rp  is  the  high-injection  carrier  lifetime  and  Que  is 
a part  of  the  stored  hole  charge  defined  as 

H' 

Qbe  = J <1^  (l  - -jyr)  p(x.t)dx  ; (A. 2) 

o 

A is  the  device  area.  With  regard  to  errors  introduced  by  the 
quasi-static  approximation,  the  most  critical  term  in  (A.l)  is 
the  transport-current  integral  which  can  convey  significant 
non-quasi-static  behavior.  Thus  the  idea  [Foss86]  was  to 
reorganize  the  continuity  equation  to  make  the  transport- 
current  term  quasi-static,  thereby  minimizing  the  error 
introduced  when  the  entire  equation  is  characterized  by  the 
quasi-static  approximation. 

For  the  high-injection  case  here,  combining  the  hole  and 
electron  current  expressions,  with  p ~ n,  yields 


If(t) 

1+6 


qAD 


dp 
' dr 


(A. 3) 


where  h is  the  electron-hole  mobility  ratio,  D i is  the 
ambipolar  diffusivity,  and  IE  represents  at  any  point  the 
sum  of  the  hole  and  electron  currents.  If  the  effective  base 
contact  for  this  one-dimensional  transistor  is  assumed  to  be 
located  near  the  base-collector  junction,  as  is  the  case  for 
the  constituent  pnp  in  the  IGBT  structure  (see  Fig.  3.1),  then 
If:  is  spatially  constant,  as  indicated  in  (A. 3),  and  it  is 
in  fact  the  emitter  current  of  the  transistor.  Consequently, 
the  transport -current  integral  in  (A.l),  with  (A. 3),  is 
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quasi-static;  i.e.,  it  is  independent  of  the  transient  hole 
distribution  in  the  quasi-neutral  base  and  depends  only  on 
the  boundary  hole  densities  which  are  virtually  quasi-static. 
Since  the  non-quasi-static  nature  of  the  recombination  term, 
Qbe/th,  in  (A.l)  is  of  lesser  significance,  we  can  assume 
that  (A.l)  can  be  well  approximated  as  being  quasi-static. 

A similar,  but  reverse  double  integration  of  the  hole 
continuity  equation  yields 


\v 

— Yp  Ip  (.v.  f)  dx  + 

o 


Qbc 

TH 


+ 


''Qbc 

dt 


(A. 4) 


Qbc  = j qA  (^jpix.t)dir  (A. 5) 
o 

is  the  other  part  of  the  stored  hole  charge.  Note  that 
Qbe  + Qbc  = Q B • By  the  same  argument  given  above,  (A.  4) 
can  be  assumed,  to  first  order,  to  be  quasi-static.  Hence, 
(A.l)  and  (A. 4) , which  characterize  the  hole  transport  in 
the  quasi-neutral  base  region,  define  a partitioning  of  the 
total  base  charge  QH  between  the  emitter  and  collector  of 
the  transistor.  This  partition,  which  defines  how  the  emitter 
and  collector  share  in  the  discharging  of  the  base,  gives  a 
first-order  accounting  for  non-quasi-static  behavior  due  to 
the  finite  carrier  transit  time  in  the  base  [Foss86] . 
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To  incorporate  the  charge  partitioning  into  the  transient 
network  model  for  the  transistor  and  the  IGBT,  we  recognize 
that  the  steady-state  component  of  base  current, 


Ib 


Qbe  Qbc 

TH  TH 


+ In  e 


(A. 6) 


where  I\p  is  the  back  injection  of  electrons  into  the  emitter, 
must  flow  between  the  emitter  and  base  terminals.  Thus  (A.l) 
and  (A.  4)  imply  that 

11 

Ic  = [ Ip  ( x.  i ) dr  — (A. . 7 ) 

' ' ./  TH 

I) 

= If  - Ib 


and  that  the  partitioned  charge-based  model  for  the  pnp 
transistor  should  be  represented  in  the  IGBT  model  as  shown 
in  Fig.  3.9,  where  Ip  and  Ic  are  defined  [Fos88]  via  the  quasi- 
static approximation  and  the  ambipolar  transport  analysis. 
This  analysis  [Fos88],  which  gives 


p(.v.t)  = p((),t) 


silih  [(TT  - .r)/L,\Y 
sivh  ( \V/L  i ) 


(A. 8) 


where  La  is  the  ambipolar  diffusion  length,  also  defines  the 
cIuasi~static  partitioned  charges  Qyp  and  Q hc  . Inserting  (A.  8) 
into  (A. 2)  and  (A. 5),  with  some  mathematical  manipulation, 
yields 


Qbe(^)  = qAp(OJ)  I i 
^ (2/3 )Qg 


(A. 9) 
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and 


QBC(t)  = qAp(0.t)LA 


'La 

, /irv 

1 

'sch  ( — ) 

tr 

U,t/J 

— ( 1/3)  Qu  . 


(A. 10) 


The  time  dependences  of  and  Qbc  are  defined  by  the 

time  dependences  of  the  emitter-base  and  base-collector 
voltages  which  control  />(()./)  and  II'  [Fos88].  The 

numerical  approximations  in  (A. 9)  and  (A. 10)  follow  from  the 
corresponding  expression  for  Q B and  the  fact  that  typically 
L \ is  greater  than  IT'/ 4 in  IGBTs . 


APPENDIX  B 

MODEL  IMPLEMENTATION  IN  SPICE  VIA  UDCS'S 

In  this  appendix  we  overview  the  SPICE  implementation  of 
the  composite  LDMOST  model  developed  in  Chapter  5.  All  of  the 
elements  in  the  network  representation  of  the  LDMOST  model 
shown  in  Fig.  5.8  are  implemented  in  SPICE  via  user-defined 
controlled  sources  (UDCSs) , which  are  FORTRAN  subroutines  in 
SLICE  [Har84],  an  enhanced  version  of  SPICE2.  UDCSs  (UUID, 
UUVEP,  UUIBR,  UUICT,  UUQP,  UUQJC  in  the  following  source-code 
listing)  access  the  differentiation  routine  DEFF  to  calculate 
the  element  values,  transconductances,  and  transcapacitances 
needed  for  circuit  nodal  analysis.  The  subroutine  DEFF 
calculates  the  derivatives  via  difference  equations  based 
on  multiple  calls  of  the  model  routine  DMT  (flowcharted  in 
Fig.  5.9).  Finally,  the  subroutine  DMT  numerically  solves  the 
implicit  system  of  equations  in  the  model  using  a straight 
iteration  method.  The  general  theory  behind  the  implementation 
of  UDCS  network  models  in  SLICE  is  discussed  in  the  SLICE 
Manual  Rev.  4.08  [Har84]  and  Veeraraghavan  et  al . [Vee86] . 

The  FORTRAN  source  code,  which  details  each  UDCS  and  model 
routine  for  the  LDD  LDMOST  model  shown  in  Fig.  5.8,  is  listed 
on  the  following  pages. 
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c ************************************************************ 

c*  SPICE2  User-defined  controlled-source  subroutine  for  the  * 
c*  channel  current  in  the  LDD  LDMOST  by  Y.-S.  Kim  3/1/90.  * 

0************************************************************ 
c*  * 

c*  Return  'y'  as  value  of  source  if  'iflag'=0,  or  as  value  * 
c*  of  'iflag'-th  component  of  gradient  of  source  if  * 

c*  'iflag'>0.  * 

c*  * 

c*  Perform  one-time  initialization  of  error  checking  if  * 

c*  'iflag'=-l.  * 

c*  * 

c*  Optionally  return  'j error ' as  integer  in  the  range  * 

» c*  1..9999  to  cause  SPICE2  to  print  an  error  message  * 

c*  containing  the  value  of  'jerror',  or  in  the  range  (-1)..  * 

c*  (-9999)  to  cause  SPICE2  to  print  an  error  message  and  * 

c*  halt  execution,  the  values  outside  these  ranges  are  * 

c*  reserved  for  use  by  SPICE2 . * 

c*  * 

c*  No  other  argument  or  variables  should  be  changed,  except  * 
c*  variables  added  by  the  user.  * 

c*  * 

c*  Description  of  input  arguments:  * 

c*  * 

c*  LP  - pointer  to  settable  parameters  * 

c*  NP  - number  of  settable  parameters  (use  is  optional)  * 

c*  The  settable  parameters  (if  any,  i.  e.  if  np>0)  are  * 

c*  referenced  as  x(lp+l) , x(lp+2) , ...,  x(lp+np) . * 

c*  * 

c*  LC  - pointer  to  controlling  arguments  (voltages  or  * 

c*  currents)  * 

c*  NC  - number  of  controlling  arguments  (use  is  optional)  * 
c*  The  controlling  arguments  are  referenced  as  * 

c*  x ( lc+1 ) , x ( lc+2 ) , ...,  x(lc+nc).  * 

c*  * 

c*  The  elements  of  both  of  the  above  lists  must  be  in  one-  * 
c*  to-one  correspondence  with  the  values  specified  on  every  * 
c*  SPICE2  element  card  referencing  this  subroutine.  * 

c*  * 

c*  The  manual  has  errors:  'special  common  blank'  line  must  * 
c*  be  deleted;  'implicit  double  precision'  must  be  added;  * 
c**************** ***********************  ******  *************** 

subroutine  UUID (y , if lag , lp , np ,1c, nc , j error) 
implicit  double  precision  (a-h,o-z) 
common  /blank/  x(64) 

common  /parm/  z,xll,xl2,unch,xna,xnae,vsatch,und,xnd, 
&xld , vsatd , gqs , ta , xiccO , xiecO , rth , h , dum 
c ************************************************************ 
c*****lp=7 s , np=18 , lc=96 , nc=3 
z=  x(lp+l) 

xll=  x ( lp+2 ) 

xl2=  x ( lp+3 ) 

unch=  x ( lp+4 ) 

xna=  x ( lp+5) 

xnae=  x(lp+6) 

vsatch=  x ( lp+7 ) 

und=  x ( lp+8 ) 
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xnd= 

x ( lp+9 ) 

xld= 

x ( lp+10) 

vsatd= 

x ( lp+11 ) 

gqs= 

x ( lp+12 ) 

ta= 

x ( lp+13 ) 

xicc0= 

x ( lp+14 ) 

xiec0= 

x(lp+15) 

rth= 

x ( lp+16 ) 

h= 

x ( lp+17 ) 

dum= 

x ( lp+18 ) 

vg= 

x ( lc+1) 

vch= 

x(lc+2) 

vben= 

x(lc+3) 

c* *********************************************************** 
c*****if lag=-l : error  checking. 
c*****iflag=  0:  evaluate  y=  (value  of  source). 
c*****iflag=  1:  evaluate  y=d(value  of  source) /dx (lc+1) . 
c*****ifiag=  2:  evaluate  y=d (value  of  source) /dx (lc+2 ) . 
c*****iflag=  3:  evaluate  y=d(value  of  source) /dx (lc+3 ) . 

goto  (25,50,100,200,300) ,iflag  +2 
c* *********************************************************** 
c*  perform  one-time  initialization  or  error  checking  * 

c* *********************************************************** 
25  continue 

if  (nc  .ne.  3)  jerror  = -99010 
goto  99901 

c* *********************************************************** 
c*  y=ID  (channel  current)  * 

c* *********************************************************** 
50  continue 

call  DEFF (1, 0,y,vg, vch,vben) 
goto  99901 

c* *********************************************************** 
c*  y=dID/dVG  * 

c* *********************************************************** 
100  continue 

call  DEFF (1,1, y , vg, vch, vben) 
goto  99901 

c* *********************************************************** 
c*  y=dID/dVch  * 

c* *********************************************************** 
200  continue 

call  DEFF ( 1 , 2 , y , vg , vch, vben) 
goto  99901 

C*  ***************************************************  * * + 

c*  y=dID/dVBEn  * 

c*  *********************************************************** 
300  continue 
y=0 

goto  99901 
c 

99901  continue 
c 

return 

end 
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c*  **********************************************************  * 

c*  SPICE2  User-defined  controlled-source  subroutine  for  the  * 
c*  drift  region  voltage  drop  in  the  LDD  LDMOST . * 

c*  * 

c*  Refer  to  the  comment  in  UUID.f.  * 

c************************************************************ 

subroutine  UUVEP(y , if lag, lp, np, lc, nc, j error) 
implicit  double  precision  (a-h,o-z) 
common  /blank/  x(64) 

common  /parm/  z,xll,xl2,unch,xna,xnae,vsatch,und,xnd, 
&xld , vsatd , gqs , ta , xiccO , xiecO , rth , h , dum 
c ************************************************************ 
c*****lp=78,  np=18 , lc=96 , nc=3 


z= 

x(lp+l) 

xll= 

x ( lp+2 ) 

xl2= 

x ( lp+3 ) 

unch= 

x ( lp+4 ) 

xna= 

x(lp+5) 

xnae= 

x(lp+6) 

vsatch= 

X ( lp+7 ) 

und= 

X ( lp+8 ) 

xnd= 

X ( lp+9 ) 

xld= 

X(lp+10) 

vsatd= 

x ( lp+ 11 ) 

gqs= 

x ( lp+12 ) 

ta= 

x ( lp+13 ) 

xicc0= 

x ( lp+14 ) 

xiec0= 

x ( 1 p+ 1 5 ) 

rth= 

x ( lp+16 ) 

h= 

x ( lp+17 ) 

dum= 

x ( lp+18 ) 

vg= 

x ( lc+1 ) 

vch= 

x(lc+2) 

vben= 

x ( lc+3 ) 

c***** ****** ******************* ****************************** 

c*****iflag=-l : error  checking. 

c*****if iag=  0:  evaluate  y=  (value  of  source) . 

c*****ifiag=  l:  evaluate  y=d(value  of  source)/dx(lc+l) . 

c*****ifiag=  2:  evaluate  y=d(value  of  source) /dx (lc+2 ) . 

c*****ifiag=  3:  evaluate  y=d(value  of  source) /dx ( lc+3 ) . 

goto  (25,50,100,200,300) , if lag  +2 
c* *********************************************************** 

c*  perform  one-time  initialization  or  error  checking  * 

c************************************************************ 
25  continue 

if  (nc  .ne.  3)  jerror  = -99010 
goto  99901 

c************************************************************ 
c*  y=Vepi  (drift  region  voltage  drop)  * 

C******************************************  *************  ****■)' 
50  continue 

call  DEFF (2,0 , y , vg, vch, vben) 
goto  99901 

C*  **************************************  * -k  * + * * * * * ££*£***£***  * 

c*  y=dVepi/dVG  * 

c ************************************************************ 
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100  continue 

call  DEFF(2 , l,y,vg, vch, vben) 
goto  99901 

C*  *********  -k  ***************  -k  * * It  ******  -k  ************  -k  ********  -k  * 

c*  y=dVepi/dVch  * 

c** *************************************  * ****** * * * * * * * * * * * * * * 

200  continue 

call  DEFF(2,2,y,vg,vch,vben) 
goto  99901 

c************************************************************ 
c*  y=dVepi/dVBEn  * 

0*********************************** ***************** ******** 

300  continue 

call  DEFF (2 , 3 ,y, vg, vch, vben) 
goto  99901 
c 

99901  continue 
c 

return 

end 


0* *********************************************************** 


c*  SPICE2  User-defined  controlled-source  subroutine  for  the  * 
c*  reverse  base  current  in  the  npn  BJT  of  the  LDD  LDMOST.  * 
c*  * 


c*  Refer  to  the  comment  in  UUID.f.  * 

0* *********************************************************** 

subroutine  UUIBR (y , if lag, lp , np , lc, nc , j error) 
implicit  double  precision  (a-h,o-z) 
common  /blank/  x(64) 

common  /parm/  z,xll,xl2,unch,xna,xnae,vsatch,und,xnd, 
&xld , vsatd , gqs , ta , xiccO , xiecO , rth , h , dum 
c*************************************************** ******* ** 


C*****lp=78 , 
z= 

xll= 

xl2  = 

unch= 

xna= 

xnae= 

vsatch= 

und= 

xnd= 

xld= 

vsatd= 

gqs= 

ta= 

xicc0= 

xiec0= 

rth= 


np=18,  lc=96 , 
x(lp+l) 

X ( lp+2 ) 

X ( lp+3 ) 

X ( lp+4 ) 

X ( lp+5 ) 

X ( lp+6 ) 

X ( lp+7 ) 
x ( lp+8 ) 

X ( lp+9 ) 

x(lp+10) 
x(lp+ll) 
x ( lp+12 ) 
x ( lp+13 ) 
x ( lp+14 ) 
x ( lp+15 ) 
x(lp+16) 


nc=3 
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h= 

dum= 
vg= 
vch= 
vben= 
c************* 
c*****iflag=-l 
c*****iflag=  0 
c*****iflag=  1 
c*****ifiag=  2 
c*****iflag=  3 
goto  (25 

c*  ***********  * 

c*  perform  one 

q*  ***********  * 


x ( lp+17 ) 
x ( lp+18 ) 
x ( lc+1) 
x(lc+2) 
x(lc+3) 

*********************************************** 
: error  checking. 

: evaluate  y=  (value  of  source) . 

: evaluate  y=d(value  of  source) /dx ( lc+1) . 

: evaluate  y=d(value  of  source) /dx ( lc+2 ) . 

: evaluate  y=d (value  of  source) /dx (lc+3 ) . 
,50,100,200,300) ,iflag  +2 
*********************************************** 

-time  initialization  or  error  checking  * 

*********************************************** 


25  continue 

if  (nc  .ne.  3)  jerror  = -99010 
goto  99901 

c* *********************************************************** 


c*  y=lBR  (reverse  base  current  of  npn  TR)  * 

c* *********************************************************** 
50  continue 

call  DEFF (3 , 0 , y , vg, vch, vben) 
goto  99901 

c************************************************************ 
c*  y=dIBR/dVG  * 

c************************************************************ 
100  continue 

call  DEFF(3 , l,y,vg, vch, vben) 
goto  99901 

c* *********************************************************** 
c*  y=dIBR/dVch  * 

c* *********************************************************** 
200  continue 

call  DEFF (3, 2, y,vg, vch, vben) 
goto  99901 

C*  **********************************************************  * 

c*  y=dIBR/dVBEn  * 

c* **********************************************************  * 
300  continue 

call  DEFF (3,3 ,y,vg, vch, vben) 
goto  99901 
c 

99901  continue 


c 


return 

end 
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C*  ★ ★ * * ★ ★ ★ ★ ★ ic  ★ * ★ ★ ★ ★ ★ ★ ★ ★ ★ * ★ ★ ★ -k  ★ ★ * ★ ★ ★ * * ★ ★ ★ ★ * * * * ★ * ★ * * * * ★ ★ ★ * ★ * * * ★ * 

c*  SPICE2  User-defined  controlled-source  subroutine  for  the  * 
c*  collector  current  in  the  npn  BJT  of  the  LDD  LDMOST.  * 

c*  * 

c*  Refer  to  the  comment  in  UUID.f.  * 

C*  ★ ★ ★ ★ * ★ * ★ ★ ★ * ★ ★ ★ ★ ★ ★ * * ★ * ★ * ★ * * * * * -k  ★ ★ ★ ★ ★ ★ ★ * ★ ★ ★ ★ ★ ★ * ★ * * * * ★ ★ ★ * * ★ ★ •*  ★ 

subroutine  UUICT(y, iflag, lp,np, lc,nc, jerror) 
implicit  double  precision  (a-h,o-z) 
common  /blank/  x(64) 

common  /parm/  z,xll,xl2,unch,xna,xnae,vsatch,und,xnd, 
&xld , vsatd , gqs , ta , xiccO , xiecO , rth , h , dum 
c* *********************************************************** 


c*****lp=7  8 , np=18 , lc=96,  nc=3 


z= 

x(lp+l) 

xll= 

x(lp+2) 

xl2= 

X ( lp+3 ) 

unch= 

x ( lp+4 ) 

xna= 

x(lp+5) 

xnae= 

x ( lp+6) 

vsatch= 

x ( lp+7 ) 

und= 

x ( lp+8 ) 

xnd= 

x ( lp+9) 

xld= 

x (lp+10) 

vsatd= 

x(lp+ll) 

gqs= 

x ( lp+12 ) 

ta= 

x ( lp+13 ) 

xicc0= 

x ( lp+14 ) 

xiec0= 

x ( lp+15) 

rth= 

x ( lp+16 ) 

h= 

x ( lp+17 ) 

dum= 

x ( lp+18 ) 

vg= 

x(lc+l) 

vch= 

x(lc+2) 

vben= 

x(lc+3) 

c* *********************************************************** 


c*****ifiag=-l:  error  checking. 
c*****ifiag=  0:  evaluate  y=  (value  of  source) . 
c*****ifiag=  1:  evaluate  y=d(value  of  source) /dx (lc+1) . 
c*****ifiag=  2:  evaluate  y=d (value  of  source) /dx (lc+2 ) . 
c*****ifiag=  3:  evaluate  y=d (value  of  source) /dx ( lc+3 ) . 
goto  (25,50,100,200,300) , iflag  +2 

c* *********************************************************** 


c*  perform  one-time  initialization  or  error  checking  * 

c ************************************************************ 
25  continue 

if  (nc  .ne.  3)  jerror  = -99010 
goto  99901 

c ************************************************************ 
c*  y=ICT  (collector  transport  current  of  npn  TR)  * 

c*  **********************************************************  * 

50  continue 

call  DEFF(4,0,y,vg,vch,vben) 
goto  99901 

c ************************************************************ 

c*  y=dICT/dVG  * 

C***************************************^-),-^^***********  *****■/( 
100  continue 
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call  DEFF(4, l,y,vg,vch,vben) 
goto  99901 

C*  **********************************************************  * 

c*  y=dICT/dVch  * 

c************************************************************ 
200  continue 

call  DEFF(4,2,y,vg,vch,vben) 
goto  99901 

c* *********************************************************** 
c*  y=dICT/dVBEn  * 

c*  **********************************************************  * 
300  continue 

call  DEFF(4,3,y,vg,vch,vben) 
goto  99901 
c 

99901  continue 
c 

return 

end 


c* *********************************************************** 
c*  SPICE2  User-defined  controlled-source  subroutine  for  the 
c*  drift  region  charging  current  in  reverse  mode  of  the 
c*  LDD  LDMOST. 
c* 

c*  Refer  to  the  comment  in  UUID.f. 
c* *********************************************************** 
subroutine  UUQP (y , iflag, lp, np, lc, nc, j error) 
implicit  double  precision  (a-h,o-z) 
common  /blank/  x(64) 

common  /parm/  z,xll,xl2,unch,xna,xnae,vsatch,und,xnd, 
&xld , vsatd , gqs , ta , xiccO , xiecO , rth , h , dum 
common  /status/  omega, time, delta, delold(7) ,ag(7) ,vt, 
&xni , egfet , mode , modedc , icalc , init f , method , iord , maxord , 
Snoncon , iterno , itemno , nosolv 
common  /knstnt/  twopi , xlog2 , xloglO , root2 , rad, boltz , 
&charge , ctok , gmin , reltol , abstol , vntol , trtol , chgtol , epso , 
Sepssil , epsox 

c*****lp=262,  np=241,  lc=286,  nc=3 
if (time. eq. 0. ) then 
z=  x(lp+l) 

xll=  x ( lp+2 ) 

Xl2=  X ( lp+3 ) 

unch=  x ( lp+4 ) 

xna=  x ( lp+5 ) 

xnae=  x(lp+6) 

vsatch=  x ( lp+7 ) 

und=  x(lp+8) 

xnd=  x ( lp+9) 

xld=  x(lp+l0) 

vsatd=  x ( lp+ll ) 

gqs=  x ( lp+12 ) 

ta=  x ( lp+13 ) 
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xiccO=  x ( lp+14 ) 
xiecO=  x ( lp+15 ) 
rth=  x ( lp+16 ) 
h=  x ( lp+17 ) 

dum=  x(lp+18) 
ynoer=  x(lp+24) 
end  if 

c* *********************************************************** 
c*****ifiag=-l:  error  checking. 
c*****ifiag=  0:  evaluate  y=  (value  of  source). 
c*****ifiag=  1:  evaluate  y=d(value  of  source) /dx ( lc+1) . 
c*****iflag=  2:  evaluate  y=d(value  of  source) /dx( lc+2 ) . 
c*****ifiag=  3:  evaluate  y=d(value  of  source) /dx ( lc+3 ) . 

goto  (25,50,100,200,300) , iflag  +2 
25  continue 

return 

c************************************************************ 
c*  Initialize  the  voltages  and  make  current  =0.0  in  DC  case  * 
c* *********************************************************** 
50  continue 

vg=  x(lc+l) 

vch=  x ( lc+2 ) 

vben=  x(lc+3) 
if (time. eq. 0 . ) then 

call  DEFF (6, 0,yy, vg, vch, vben) 

qp=yy 

x(lp+19)=qp 

y=o. 

x ( lp+20) =y 
x(lp+21)=delta 
x ( lp+22 ) =qp 
x ( lp+23 ) =y 
return 
else 

c***** ******************************************************* 
c*  Update  old  charge  and  current,  iteration  at  a new  time 
c*  point  if  local  truncation  error  criteria  are  met 
c* ***************** ****************************************** 

if ( (initf.eq.6) .and. (delta. ge.x(lp+21) ) ) then 
x(lp+19) =x( lp+22) 
x ( lp+20) =x ( lp+23 ) 
end  if 

call  DEFF(6, 0,yy,vg, vch, vben) 
qp=yy 

x ( lp+22 ) =qp 
x ( lp+21) =delta 
if (delta. eq. 0. ) then 
y=0. 
else 

c*********Use  Backward  Euler  approximation  for  numerical 
c********* integration 

c* ******* *x( lp+19) =x ( tn) ,x (lp+22 )=x(tn+l) 
c*********x (lp+20) =y(tn) ,x (lp+23 )=y(tn+l) 
c*********x (lp+21) =delta (tn) , delta=delta (tn+1) 
c*********y (tn+1) =[x (tn+1) -x(tn) ]/delta 
y=(qp-x( lp+19) ) /delta 

c*********use  Trapezodial  approximation  for  numerical 
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c* ******** integration 

c*********y(tn+l)=2*[x(tn+l) -x(tn) ]/delta-y (tn) 
c y=(qp-x(lp+19) ) *2 . 0/del ta-x ( lp+20) 

end  if 
x ( lp+23 ) =y 
end  if 

c****************************************************** ****** 
c*  IF  ynoer=0  then  no  error-checking. 

c*  Error  checking  is  done  presuming  that  SLICE/SPICE2  does 
c*  check  for  current  convergence.  If  the  current  converges 
c*  and  not  the  charge,  then  this  program  increments  the 
c*  parameter  noncon. 

c********************************************* ******  ********* 
if (ynoer.ne. 0 . ) then 

toli=reltol*dmaxl (dabs (y) , dabs (x (lp+20) ) ) +abstol 
if ( (dabs (y-x(lp+20) ) ) . It . toli) then 

tol=reltol*dmaxl (dabs (qp) ,dabs (x(lp+19) ) ) +chgtol 
if ( (dabs (qp-x (lp+19) ) ) .ge. tol)  noncon=noncon+l 
end  if 
end  if 

c no  error  checking 

c ynoer=0 

return 

c* *********************************************************** 
c*  Derivative  w.r.t.  VG.  y=dQP/dVG 

c* *********************************************************** 
100  continue 

if (delta. eq. 0. ) then 
y=0. 
else 

call  DEFF (6,1, yy , vg , vch , vben ) 
y=2 . *yy/delta 
end  if 
return 

c* *********************************************************** 
c*  Derivative  w.r.t.  VCH.  y=dQP/dVCH 

c* *********************************************************** 
200  continue 

if (delta. eq.0. ) then 
y=0. 
else 

call  DEFF (6, 2 ,yy, vg, vch, vben) 
y=2 . *yy/delta 
end  if 
return 

c ************************************************************ 
c*  Derivative  w.r.t.  VBEn . y=dQP/dVBEn 

c***************************************************** ******* 

300  continue 

if (delta. eq. 0. ) then 
y=0. 
else 

call  DEFF (6, 3 ,yy, vg, vch, vben) 
y=2 . *yy/delta 
end  if 
return 
end 
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c*  SPICE2  User-defined  controlled-source  subroutine  for  the  * 
c*  LDD-drift  space-charge  region  displacement  current  of  the* 
C*  LDD  LDMOST . * 
c*  * 


c*  Refer  to  the  comment  in  UUID.f.  * 

c***************************************************** ******* 
subroutine  UUQJC (y , iflag, lp, np, lc, nc, j error) 
implicit  double  precision  (a-h,o-z) 
common  /blank/  x(64) 

common  /parm/  z,xll,xl2,unch,xna,xnae,vsatch,und,xnd, 
&xld , vsatd , ggs , ta , xiccO , xiecO , rth , h , dum 
common  /status/  omega , time , delta , delold ( 7 ) , ag ( 7 ) , vt , 
&xni , egf et , mode , modedc , icalc , init f , method , iord , maxord , 
Snoncon , iterno , itemno , nosolv 
common  /knstnt/  twopi , xlog2 , xloglO , root2 , rad, bolt z , 
Scharge, ctok,gmin, reltol , abstol , vntol , trtol , chgtol , epso, 
Sepssil , epsox 

c*****lp=2  62 , np=241 , lc=286,  nc=3 
if (time. eg. 0. ) then 


z= 

x(lp+l) 

xl  1= 

x ( lp+2 ) 

Xl2= 

x ( lp+3 ) 

unch= 

x ( lp+4 ) 

xna= 

x(lp+5) 

xnae= 

x ( lp+6) 

vsatch= 

X ( lp+7 ) 

und= 

X ( lp+8 ) 

xnd= 

X ( lp+9 ) 

xld= 

X(lp+10) 

vsatd= 

x(lp+ll) 

ggs= 

x ( lp+12 ) 

ta= 

x ( lp+13 ) 

xicc0= 

x ( lp+14 ) 

xiec0= 

x(lp+15) 

rth= 

x (lp+16) 

h= 

x ( lp+17 ) 

dum= 

x ( lp+18 ) 

ynoer= 

x ( lp+24 ) 

end  if 

c***********  *************************  **************i,^.)e.)C.f(it^^i(it 
c*****ifiag=-l:  error  checking. 
c*****ifiag=  0:  evaluate  y=  (value  of  source) . 
c*****ifiag=  l:  evaluate  y=d(value  of  source) /dx (lc+1) . 
c*****ifiag=  2:  evaluate  y=d(value  of  source)/dx(lc+2) . 
c*****ifiag=  3:  evaluate  y=d(value  of  source)/dx(lc+3  . 

goto  (25,50,100,200,300) , iflag  +2 
25  continue 

return 

c* *********************************************************** 

c*  Initialize  the  voltages  and  make  current  =0.0  in  DC  case  * 
c ************************************************************ 
50  continue 


vg=  x ( lc+1 ) 

vch=  x ( lc+2 ) 

vben=  x ( lc+3 ) 

if (time. eg. 0 . ) then 
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call  DEFF (5 , 0 , yy , vg, vch, vben) 
q j c=yy 

x( lp+19 )=qjc 
y=0. 

x(lp+20)=y 
x(lp+21) =delta 
x ( lp+22 ) =qj  c 
x ( lp+23 ) =y 
return 
else 

C*  ****************  k *****************************************  * 

c*  Update  old  charge  and  current,  iteration  at  a new  time 
c*  point  if  local  truncation  error  criteria  are  met. 
c*  **********************************************************  * 
if ( ( initf . eq. 6) . and. (delta . ge . x (lp+21) ) ) then 
x ( lp+19 ) =x ( lp+22 ) 
x(lp+20) =x( lp+23) 
end  if 

call  DEFF(5, 0,yy,vg, vch, vben) 
q j c=yy 

x( lp+22 )=qjc 
x (lp+21) =delta 
if (delta. eq. 0. ) then 
y=0. 
else 

c*********use  Backward  Euler  approximation  for  numerical 
c********* integration 

c*********x( lp+19) =x(tn) ,x (lp+22 )=x(tn+l) 
c*********x(lp+20)=y (tn) , x ( lp+23 ) =y (tn+1) 
c********+x( lp+21) =delta (tn) , delta=delta (tn+1) 
c*********y (tn+1) =[x( tn+1) -x (tn) ]/delta 

y=(qjc-x( lp+19) ) /delta 

c*********use  Trapezodial  approximation  for  numerical 
c*********integration 

c*********y (tn+1) =2* [x (tn+1) -x (tn) ]/delta-y (tn) 
c y= ( (qj c-x ( lp+19 ) ) *2 . 0/delta) -x ( lp+2  0 ) 

end  if 
x (lp+23 )=y 
end  if 

C*  **************  k *******************************************  * 

c*  IF  ynoer=0  then  no  error-checking. 

c*  Error  checking  is  done  presuming  that  SLICE/SPICE2  does 
c*  check  for  current  convergence.  If  the  current  converges 
c*  and  not  the  charge,  then  this  program  increments  the 
c*  parameter  noncon. 

C*  ******  k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k 
if (ynoer.ne. 0. ) then 

toli=reltol*dmaxl (dabs (y) ,dabs (x(lp+20) ) ) +abstol 
if ( (dabs (y-x ( lp+20) ) ) . It. toli) then 

tol=reltol*dmaxl(dabs(qjc) , dabs (x ( lp+19 ) ) ) +chgtol 
if ( (dabs (qj c-x (lp+19) ) ) .ge.tol)  noncon=noncon+l 
end  if 
end  if 

no  error  checking 
ynoer=0 
return 


c 

c 
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c*  **********************************************************  * 

c*  Derivative  w.r.t.  VG.  y=dQJC/dVG 

q*  **********************************************************  * 

100  continue 

if (delta. eq. 0. ) then 

y=o. 

else 

call  DEFF (5,1, yy , vg , vch , vben ) 
y=2 . *yy/delta 
end  if 
return 

d************************************************************ 

c*  Derivative  w.r.t.  VCH.  y=dQJC/dVCH 

c* **********************************************************  * 

200  continue 

if (delta. eq. 0 . ) then 
y=0. 
else 

call  DEFF(5,2,yy,vg,vch,vben) 
y=2 . *yy/delta 
end  if 
return 

c* *********************************************************** 

c*  Derivative  w.r.t.  VBEn . y=dQJC/dVBEn 

c ************************************************************ 

300  continue 

if (delta. eq.O. ) then 
y=0. 
else 

call  DEFF (5, 3 ,yy, vg, vch, vben) 
y=2 . *yy/delta 
end  if 
return 
end 


c ************************************************************ 


c*  SPICE2  User— defined  controlled— source  subroutine  * 
c*  for  differentiating  the  source.  * 
c*  * 
c*  Format  to  call  this  subroutine  in  UU%%  file:  * 


c*  1)  to  calculate  %%;  call  DEFF(#, 0,y,vg, vch, vben)  * 

c*  2)  to  calculate  d%%/dvg;  call  DEFF (#, l,y,vg, vch, vben)  * 

c*  3)  to  calculate  d%%/dvch;  call  DEFF(# , 2 ,y, vg, vch, vben)  * 

c*  4)  to  calculate  d%%/dvben;  call  DEFF ( # , 3 , y , vg, vch, vben)  * 

c*  where  #=1  for  %%=ID,  #=2  for  %%=VEP,  #=3  for  %%=IBR,  * 

c*  #=4  for  %%=ICT , #=5  for  %%=dQJC/dt  #=6  for  %%=dQP/dt.  * 
c ************************************************************ 
subroutine  DEFF (mf lgl , mf lg2 , y , vg, vch , vben) 
implicit  double  precision  (a-h,o-z) 
common  /parm/  z,xll,xl2,unch,xna,xnae,vsatch,und,xnd, 
&xld , vsatd , gqs , ta , xiccO , xiecO , rth , h , dum 
common  /status/  omega, time, delta, delold(7) ,ag(7) ,vt, 
&xni , egfet , mode , modedc , icalc , init f , method , iord , maxord , 
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Snoncon , iterno , itemno , nosolv 
common  /poggy/  vicon(lO) , itprev, viconl (10) ,vicon2 (10) , 
&vicon3 (10) ,vicon4 (10) 

c* *********************************************************** 
c h=ld-3 

vbcn=-vch 

c* *********************************************************** 

c*  This  part  speeds  up  the  program. 

c*  ' iterno'  is  the  iteration  number  in  the  SPICE. 

c*  'itprev'  is  the  previous  iteration  number. 

c* 

c*  At  any  iteration  number  (iterno>itprev) , once  one  of  the 
c*  UDCSs  calls  DEFF  which  in  turn  calls  DMTs , then  it  sets 
c*  itprev=iterno. 
c* 

c*  When  other  UDCSs  in  the  same  iteration  number  call  DEFF 
c*  (iterno=itprev) , it  skips  calling  of  DMTs. 
c* 

c*  In  case  of  more  than  one  devices,  which  has  several 
c*  different  sets  of  controlling  voltages,  UDCSs  with 
c*  different  controlling  voltages  need  calling  of  DMTs  even 
c*  though  iterno=itprev.  For  this  case  we  need  to  compare 
c*  controlling  voltages. 

c************************************************************ 
if ( iterno . eg . 0 ) then 
if (mflg2 .eg. 0) then 

call  DMT ( xid , vep , xibr , xict , g j c , qp , qgd , xi i , vg , vch , vben ) 

vicon(l) =xid 

vicon(2) =vep 

vicon(3)=xibr 

vicon(4)=xict 

vicon (5) =qjc 

vicon (6) =qp 

vicon(7)=qgd 

vicon (8) =xii 

end  if 

if (mflg2 .eq. 1) then 
vg=vg+0 . 01+h 

call  DMT (xid , vep , xibr , xict , qj c , qp , qgd , xii , vg , vch , vben ) 

viconl (1) =xid 

viconl (2 )=vep 

viconl (3) =xibr 

viconl (4) =xict 

viconl (5) =qjc 

viconl (6) =qp 

viconl (7) =qgd 

viconl (8) =xii 

vg=vg-0 . 01-h 

end  if 

if (mf lg2 . eq. 2 ) then 
vch=vch+h 

call  DMT (xid, vep, xibr , xict , qjc, qp , qgd, xii , vg, vch , vben) 

vicon2 (1) =xid 

vicon2 (2) =vep 

vicon2 (3)=xibr 

vicon2 (4) =xict 
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vicon2 (5) =qjc 
vicon2 (6) =qp 
vicon2 (7) =qgd 
vicon2 (8) =xii 
vch=vch-h 
end  if 

if (mf lg2 . eq. 3 ) then 
vben=vben+h 

call  DMT ( xid , vep , xibr , xict , q j c , qp , qgd , xi i , vg , vch , vben ) 

vicon3 (1) =xid 

vicon3 (2) =vep 

vicon3 (3)=xibr 

vicon3 (4)=xict 

vicon3 (5) =qjc 

vicon3 (6) =qp 

vicon3 (7) =qgd 

vicon3 (8) =xii 

vben=vben-h 

end  if 

c itprev=iterno 

c end  if 

c if (iterno.gt. itprev) then 

del=ld-6 

else  if ( (dabs (vg-vgp) .gt.del) .or. 

(dabs (vch-vchp) .gt.del) .or. 

(dabs (vben-vbenp) [. gt.del)  ) then 
call  DMT ( xid , vep , xibr , xict , q j c , qp , qgd , xi i , vg , vch , vben ) 
vicon (l)=xid 
vicon(2) =vep 
vicon(3)=xibr 
vicon(4)=xict 
vicon (5) =qjc 
vicon(6)=qp 
vicon (7) =qgd 
vicon (8) =xii 
vg=vg+0 . 01+h 

call  DMT ( xid , vep , xibr , xict , qj c , qp , qgd , xi i , vg , vch , vben ) 
viconl (1) =xid 
viconl (2) =vep 
viconl (3 )=xibr 
viconl ( 4 ) =xict 
viconl (5) =qjc 
viconl (6) =qp 
viconl (7) =qgd 
viconl (8) =xii 
vg=vg-0. 01-h 
vch=vch+h 

call  DMT ( xid , vep , xibr , xict , q j c , qp , qgd , xi i , vg , vch , vben ) 
vicon2 (1) =xid 
vicon2 (2) =vep 
vicon2 ( 3 ) =xibr 
vicon2 (4) =xict 
vicon2 (5) =qjc 
vicon2 (6) =qp 
vicon2 (7) =qgd 
vicon2 (8) =xii 
vch=vch-h 
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vben=vben+h 

call  DMT ( xid , vep , xibr , xict , q j c , qp , qgd , xi i , vg , vch , vben ) 
vicon3 (1) =xid 
vicon3 (2) =vep 
vicon3 (3)=xibr 
vicon3 (4)=xict 
vicon3 (5) =qjc 
vicon3 (6) =qp 
vicon3 (7) =qgd 
vicon3 (8) =xii 
vben=vben-h 
c itprev=iterno 

end  if 

c*****end  of  speeding  block********************************** 
c 50  continue 

c************************************************************ 
c*  Source  value;  y=ID, VEP, IBR, ICT, QJC, QP, QGD, II  * 

c* *********************************************************** 
if (mflg2 . eq. 0) then 

goto (110, 120, 130, 140, 150, 160, 170, 180) ,mflgl 
110  y=vicon(l) 
goto  400 
120  y=vicon(2) 
goto  400 
130  y=vicon(3) 
goto  400 
140  y=vicon(4) 
goto  400 
150  y=vicon(5) 
goto  400 
160  y=vicon(6) 
goto  400 
170  y=vicon(7) 
goto  400 
180  y=vicon(8) 
goto  400 
400  continue 
goto  9999 
end  if 


c :************************************************************ 

c*  Gradient  of  source ; 

c*  y=d ( ID , VEP , IBR , ICT , QJC , QP , QGD , II ) /d ( VG , Vch , VBEn , VBCn , ) 
c ************************************************************ 

if (mf lg2 . ne . 0) then 
xidl=  vicon(l) 
vepl=  vicon(2) 
xibrl=vicon(3) 
xictl=vicon ( 4 ) 
qjcl=  vicon(5) 
qpl=  vicon(6) 
qgdl=  vicon(7) 
xiil=  vicon(8) 

goto (510, 520, 530, 540, 550, 560, 570, 580) ,mflgl 
if (mf lg2 . eq. 1) then 

y=(viconl (1) -xidl) / (0 . 01+h) 
else  if (mflg2.eq.2) then 
y=(vicon2 (1) -xidl)/h 


510 


520 

530 

540 

550 

560 

570 

580 
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else  if (mflg2 . eq. 3 ) then 
y=(vicon3 (1) -xidl)/h 
end  if 
goto  800 

if (mflg2 . eq. 1) then 

y=(viconl (2) -vepl) / (0. 01+h) 
else  if (mflg2 . eq. 2) then 
y=(vicon2 (2) -vepl)/h 
else  if (mf lg2 . eq. 3 ) then 
y= ( vicon3 ( 2 ) -vepl ) /h 
end  if 
goto  800 

if (mflg2 .eq. 1) then 

y=(viconl (3) -xibrl) / (0 . 01+h) 
else  if (mflg2 . eq. 2 ) then 
y= (vicon2 ( 3 ) -xibrl) /h 
else  if (mf lg2 . eq. 3 ) then 
y=(vicon3 (3) -xibrl) /h 
end  if 
goto  800 

if (mf lg2 . eq. 1) then 

y= (viconl ( 4 ) -xictl) / ( 0 . 01+h) 
else  if (mflg2 . eq. 2) then 
y=(vicon2 (4) -xictl) /h 
else  if (mflg2 .eq. 3) then 
y=(vicon3 (4) -xictl)/h 
end  if 
goto  800 

if (mflg2 .eq. 1) then 

y= (viconl (5) -qj  cl) / ( 0 . 01+h) 
else  if (mflg2 . eq. 2) then 
y=(vicon2 (5) -qjcl) /h 
else  if (mflg2.eq.3)then 
y=(vicon3 (5) -qjcl) /h 
end  if 
goto  800 

if (mf lg2 . eq. 1) then 

y= (viconl (6) -qpl)/ (0. 01+h) 
else  if (mflg2 . eq. 2 ) then 
y=(vicon2 (6) -qpl)/h 
else  if (mflg2.eq.3) then 
y=(vicon3 (6) -qpl)/h 
end  if 
goto  800 

if (mf lg2 . eq. 1) then 
y= (viconl (7 )-qgdl)/(0. 01+h) 

else  if (mflg2.eq.2)then 
y=(vicon2 (7) -qgdl)/h 
else  if (mflg2 .eq. 3) then 
y= (vicon3 (7) -qgdl)/h 
end  if 
goto  800 

if (mf lg2 . eq. 1) then 

y= (viconl (8) -xiil) / ( o . 01+h) 
else  if (mflg2 . eq. 2) then 
y=(vicon2 (8) -xiil)/h 
else  if (mflg2 . eq. 3) then 
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y= (vicon3 ( 8 ) -xiil) /h 
end  if 
end  if 
c 

800  continue 
c 

9999  continue 
vgp=vg 
vchp=vch 
vbenp=vben 
return 
end 


c************************************************************ 
c*  This  program  is  the  LDD  LDMOST  model  subroutine  which  is  * 
c*  accessed  by  the  LDD  LDMOST  SPICE  model  UDCSs  * 

c***************************** *********************  ********** 

c*  The  format  to  call  this  subroutine;  * 

c*  call  DMT(xid,vdr,xibr,xict,qj , qp, qgd, vg, vch, vben) . * 

C*  **************************************************  ie  *******  * 
subroutine  DMT ( xid , vdr , xibr , xict , qj , qp , qgd , xii , vg , vch , 
&vben) 

implicit  double  precision  (a-h,o-z) 
common  /parm/  w,xll,xl2,unch,xna,xnas,vsatch,und,xnd, 
&xld , vsatd , dvt , ta , xiccO , xiecO , rth , h , dum 
c* **********************************************************  * 

c*  Device  parameters  with  physical  constants.  * 

c*  **********************************************************  * 


esi= 

q= 

vt0= 
xni= 
maxitl= 


1.0d-12 
1 . 6d-19 
0.0259 
1 . 45dl0 
50 


c * * * * *channel ***************** 
vfb=  -0.92 
pb=  0.39 

cox=  dum 

c cox=  3 . 452d-8 

theta=  2 . 83d-2 

eta=  9 . 0 

c*****depletion  MOST********** 
pba=  0.23 

vfba=  -0.76 

c*****drift  region************ 
bu=  2 . 9 

yj=  2 . ld-4 

pbd=  0.29 

vfbd=  -0.26 
xlo=  2 . Od-4 

c* *********************************************************** 

c*  Account  for  temperature  dependence  of  vsat  and  u. 
c*  Muller  & Kamins  p37:  vl->vsat,  vl/Ec->u. 
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c************************************************************ 

tam=300 

tj=300 

do  1000  k=l,maxitl 

tu=l . 01d6*dexp (-2 . 42*dlog (tj ) ) 

tv=143.*dexp(-0.87*dlog(tj) ) 

vtt=vtO*t j/tam 

uncht=unch*tu 

undt=und*tu 

vsatcht=vsatch*tv 

vsatdt=vsatd*tv 

da=2 . *vtt*undt/ (bu+1 . ) 

xla=dsqrt (da*ta) 

c* *********************************************************** 
c*  Channel  current  for  given  VG  and  Vch 

c*  vch  is  the  total  voltage  drop  across  the  whole  channel, 
c*  vl  is  the  voltage  drop  across  the  double-df fused  channel, 
c*  *********************************************************** 
qdlO=-dsqrt ( 4 . *q*xna*esi*pb) 
qd20=-dsqrt (4 . *q*xnas*esi*pba) 
vt=vfb+2 . *pb-qdlO/cox+dvt 
c*****when  vgcvt  id=vdr=small_value (subthreshold) 
if(vg.le.vt)  then 
xid=0 . 
goto  50 
end  if 

cdlO=-q*xna*esi/qdlO 
cd2=-q*xnas*esi/qd20 
xk=cox* (vg-vt) +2 . *eta*cdlO*pb 
unef f=uncht/ ( 1 . +theta* (vg-vt) ) 
ec=2 . *vsatcht/unef f 
qn0=-cox* (vg-vt) 
c*****calculation  of  VI, sat 
xidsat=-w*vsatcht*qnO 
a= (cox+cdlO) 

b=-xk+qn0+xidsat*2 ./w/unef f/ec 
c=xidsat*2 . *xll/w/unef f 
bbac=b*b-4 . *a*c 
if (bbac. It. 0)  then 

print  *, 'During  cal.  of  vl,sat,  D<0' 
print  No  posssible  solution  for  vl,sat' 
print  *, ' Try  different  vsat(ch)' 
stop 
end  if 

vlsat=(-b-dsqrt (bbac) )/2./a 
c*****calculation  of  Vch, sat 
a=(cox+cd2) 

b=-2 • *xk+2 . * (cdl0-cd2) *vlsat+xidsat*2 ./w/uneff/ec 
c=2 . *xk*vlsat+ (-cox-2 . *cdl0+cd2 ) *vlsat*vlsat 
& +xidsat*2 . * (xl2-vlsat/ec) /w/unef f 
bbac=b*b-4 . *a*c 
if (bbac. It . 0)  then 

print  *, 'During  cal.  of  Vch, sat,  D<0' 
stop 
end  if 

vchsat=(-b-dsqrt (bbac) )/2 ./a 
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c*****when  vch<0,  NEED  to  DEVELOP  a MODEL**** 
if(vch.lt.O)  vch=-vch 
if (vch. lt.vchsat)  then 
bO=l ./ec* (-cdl0+cd2 ) 

bl=xll* (-cox-2 . *cdl0+cd2) -xl2* (cox+cdlO) 

& +l./ec* ( (cdlO-2 . *cd2-cox) *vch+xk+qnO) 

b2=xll* (2 . * (cdl0-cd2) *vch+2 . *xk) +xl2* (xk-qnO) 

& +vch/ec* ( -xk-qnO+ ( cox+cd2 ) *vch) 

b3=xll* (-2 . *xk*vch+ (cox+cd2 ) *vch*vch) 

al=bl/bO 

a2=b2/b0 

a3=b3/b0 

vl=cubic (al , a2 , a3 ) 
xid0=w*uneff/2 ./xll/ (l.+vl/ec/xll) 
xid=xidO*vl* (-2 . *qn0+2 . *eta*pb*cdlO- (cox+cdlO) *vl) 
else 

xid=xidsat 
end  if 

if (vch. It. 0)  then 
xid=-xid 
vch=-vch 
end  if 

c*****end  of  calculation  of  channel  curent. 

50  continue 

c* *********************************************************** 
c*  NORMAL  MODE  * 

c*  Calculation  of  LDD-region  voltage  drop  * 

c ************************************************************ 
if(vch.ge.O)  then 
if (dabs(xid) .le.l.d-9)  then 
vdr=0 
goto  300 
end  if 

ecd=2 . *vsatdt/undt 
vbi=vtt*dlog(xnd*xnas/xni/xni) 
xne=xnd* (xnd+xnas) /xnas 
de=l./yj*dsqrt(2 . *esi/q/xne) 
c*****determining  vdl,xidl 

c*****Above  xidl,  LDD  region  is  divided  into  two  sections, 
c*****a  and  B. 

C*****carriers  in  section  A has  field-dependent  velocity. 
C*****carriers  in  section  B has  saturated  velocity. 
al=-3./de 

a2=-3 . * (xld*ecd-vbi) 

a3=-4 . *vbi*dsqrt (vbi) +3 ./de*vbi+3 ./de*xld*ecd 

x=cubic(al,a2,a3) 

vdl=x*x-vbi 

if  (vdl.le.O)  then 

print  *,'vdl<0  in  cal.  of  idl' 
stop 
end  if 

xidl=w*yj  * ( 1 . -de*dsqrt (vdl+vbi) ) *xnd*q*vsatdt 
c*****determining  xid2 

c*****Above  xid2 , LDD  region  has  only  section  B. 
xk2=0 . 5 
al=-3 . /de 
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a2=-3 . * (xk2*yj*ecd-vbi) 

a3=-4  . *vbi*dsqrt  (vbi)  +3  ./de*vbi+3  ./de*xk2*yj  *ecd 
x=cubic (al , a2 , a3 ) 
vc2=x*x-vbi 
if  (vc2 . le. 0)  then 

print  *,'vc2<0  in  cal.  of  id2' 
stop 
end  if 

xid2=w*yj  * ( 1 . -de*dsqrt (vc2+vbi) ) *xnd*q*vsatdt 
c*****i-v  for  region  I 

if  (xid. le. xidl)  then 

al= (xid/2 ./w/yj/xnd/q/vsatdt-1 . ) *3./2./de 
a3= (xid*xld/w/yj/q/undt/xnd* ( 1 . -vbi/xld/ecd) 

& -2 ,/3 . *de*vbi*dsqrt (vbi) +vbi) *1 . 5/de 

x=cubic(al, 0,a3) 
vdr=x*x-vbi 

c*****i-v  for  region  II  and  III 

c*****vel  is  the  voltage  at  the  boundary  bet.  of  section 
c*****a  and  B. 

c*****xlel  is  the  length  of  section  A. 
c*****yel  is  the  depth  of  section  B. 
c*****el  is  the  electric  field  at  the  boundary, 
else 

if  (xid. Ie.xid2)  then 

vel= (yj -xid/w/q/xnd/vsatdt) * (yj -xid/w/q/xnd/vsatdt ) * 
& q*xne/2 ./esi-vbi 
if  (vel.lt.O)  vel=0 

xlel=w*yj *q*undt*xnd/xid* ( vel-2 . /3 . *de* ( 

& (vel+vbi) *dsqrt (vel+vbi) -vbi*dsqrt (vbi)  ) ) -vel/ecd 
else 
vel=vc2 
xlel=xk2*yj 
end  if 

c yel=l. 0* (yj-dsqrt (2 . 0*esi/q/xne* (vdl+vbi) ) ) 

el=ecd 

eyl=dsqrt (2 . 0*q/esi*vel/ ( 1 . 0/xnd+l . 0/xnas) ) 
vdr=2 . 0*vel 
do  200  i=l,50 

eyd=dsqrt ( 2 . 0*q/esi*vdr/ ( 1 . 0/xnd+l . 0/xnas ) ) 
ej=(eyl+eyd)/2 . 0 

ve2=(xld-xlel) * (xld-xlel)/2 . O/esi/yj/w/vsatdt* 

& (xid-w*q*xnd*yj  *vsatdt+esi*ej  *w*vsatdt ) + 

& (xld-xlel) *el 

vdr=vel+ve2 

if (i • eq. 50)  print  *,'NO  CONV.  in  Vep' 
if (dabs (vdr-vdro) .It. 0.1)  goto  210 
vdro=vdr 
200  continue 
210  continue 
end  if 

c*****end  of  calcualtion  of  V(drift  reqion) 

300  continue 

qj=q*w* (xld+yj ) *xnas*dsqrt (dabs ( 

& _q ’ 0*esi*xnd/q/xnas/ (xnas+xnd) * (vch+vdr+2 . *pbd)  )) 

xibr=0 

c*****end  of  normal  mode 
end  if 
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c************************************************************ 

c*  REVERSE  MODE 

c*  IBR( reverse  base  current)  is  predominantly  the 

c*  recombination  current  in  the  substrate  base  region. 

c* 

c*  ICT  is  the  base  transport  current, 
c* 

c*  QP(drift  charge)  is  given  by  integrating  n(x)  over  the 
c*  substrate  base  region, 
c* 

c*  Vdr (drift  region  voltage  drop)  is  given  by  integrating 
c*  E(x)  (from  Jn  and  Jp)  over  the  substrate  base  region, 
c* 

c*  Width  of  the  substrate  base  is  xb=xl2+xld. 
c*  Area  of  the  substrate  base  is  sb=w*yj . 

c* *********************************************************** 
xr=0 . 2 

vj pp=-xr*vch 

xict=xiccO*dexp (vben/vtt ) -xiecO*dexp (2 . 0*vjpp/vtt) 
if(vch.lt.O)  then 
xb=xl2+xld 
sb=w*yj 
wx=xb/xla 
gj=0 

xnO=xnas*dexp (vjpp/vtt) 
xnwl=xla*dsinh(wx) * (-xict) /q/sb/da/ 

& (bu*dcosh(wx) +1. ) 

xnw2= (bu+dcosh (wx) )/(bu*dcosh(wx)+l. ) *xnO 
xnw=xnwl+xnw2 

qp=q*sb*xla/dsinh (wx) * (dcosh (wx) -1 . ) * (xnO+xnw) 
xibr=qp/ta 

c*******calculating  Vepi 

if (xnO . It . xnas)  xnO=xnas 
if (xnw. It . xnas)  xnw=xnas 

xi=(l*+bu) *q*sb*da/xla/dsinh(wx) * (xnw*dcosh (wx) -xnO) 
pp=(xnw-xnO*dexp(-wx) )/2 ./dsinh(wx) 
qq= (xnO*dexp (wx) -xnw) /2 . /dsinh (wx) 

if  (pp*qq. eq. 0)  print  *, 'warning! ! ! pp*qq=o  in  DMT . f ' 
qqpp=qq/pp 
if  (qqpp.lt.O)  then 
dqp=dsqrt ( -qq/pp ) 
xp=xla/2./dsqrt(-pp*qq) *( 

& dlog (dabs ( (dexp (wx) -dqp) / (dexp (wx) +dqp) ) ) 

& -dlog (dabs ( (1. -dqp)/ (l.+dqp) ) ) ) 

else 

dpq=dsqrt(pp/qq) 

xp=xla/ dsqrt (pp*qq) * ( datan (dpq*dexp (wx) ) - 
& datan (dpq)  ) 

end  if 

vdrO=xi/q/sb/ ( 1 . +1 ./bu) /undt*xp 
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vdrl=vtt* (bu-1. ) / (bu+1 . ) *dlog (dabs (xnO/xnw) ) 
vdr2=vtt*dlog(dabs (xnw*xnas/xni/xni) ) 
vdr= ( 1-xr) /xr*vjpp-vdr0-vdrl-vdr2 
end  if 

c*****end  of  reverse  mode 
qgd=0 
xii=0 

if (k. eq.maxitl)  print  *,'No  Temp,  convergence' 
tjp=tj 
vd=vch+vdr 
xit=xid-xibr 
tj=dabs (tam+rth*vd*xit) 
if (dabs (tj-tjp) .It. 1.0)  goto  1100 
1000  continue 
1100  continue 
return 
end 

c*****function  to  solve  for  cubic  eq. 
function  cubic (al , a2 , a3 ) 
implicit  double  precision  (a-h,o-z) 
qq= (3.*a2-al*al)/9. 

rr= (9 . *al*a2-27 . *a3-2 . *al*al*al) /54 . 
if  ( (qq*qq*qq+rr*rr) .gt. 0)  then 

print  *, 'During  function  call  "cubic",  D>0' 
stop 
end  if 

if  (qq.gt.0)  then 

print  *, 'During  function  call  "cubic",  qq>0' 
stop 
end  if 

cth=rr/dsqrt ( ~qq*qq*qq) 
th=dacos (cth) 

cubic=2 . *dsqrt (-qq) *dcos (th/3 .+4 ./3 . *3 . 14159) -al/3 . 

return 

end 


REFERENCES 


[App79]  Appels,  J.  A.  and  Vaes,  H.  M.  J.,  "High  voltage 
thin  layer  devices  (RESURF  devices),"  in  IEDM  Tech. 
Dig.  , pp.  238-241,  1979. 

[Bal87]  Baliga,  B.  J.,  Modern  Power  Devices.  New  York: 

Wiley,  1987. 

[Bali84  ] Baliga,  B.  J.,  Adler,  M.  S.,  Love,  R.  P., 

Gray,  P.  V.  , and  Zommer,  N.  D.,  "The  insulated 

gate  transistor:  A new  three-terminal  MOS- 

controlled  bipolar  power  device,"  IEEE  Trans. 
Electron  Devices,  vol . ED-31,  pp.  821-828,  June 

t wr. 

[Cho88]  Chow,  T.  P.,  Baliga,  B.  J.,  and  Pattanayak,  D. 

N.,  "Counterdoping  of  MOS  channel  (CDC) -A  new 
technique  of  improving  suppression  of  latching  in 
insulated  gate  bipolar  transistors,"  IEEE  Electron 
Device  Lett . , vol.  9,  pp.  29-31,  Jan.  1988. 

[Chow87]  Chow,  T.  P.,  Pattanayak,  D.  N.,  Baliga,  B.  J., 

and  Adler,  M.  S.,  "Latching  in  lateral  insulated 
gate  bipolar  transistors,"  in  IEDM  Tech.  Dig.,  pp. 
774-777,  1987.  

[Cla86]  Claessen,  H.  R.  and  Van  Der  Zee,  P.,  "An  accurate 
DC  model  for  high-voltage  lateral  DMOS  transistors 
suited  for  CAD,"  IEEE  Trans.  Electron  Devices,  vol. 
ED-33,  pp.  1964-T970,  Dec.  

[Col81 ] Colak,  S.,  "Effects  of  drift  region  parameters  on 
the  static  properties  of  power  LDMOST, " IEEE  Trans. 
Electron  Devices,  vol.  ED-28,  pp.  1455-1466,  Dec. 


[Cou88  ] Coughran,  W.  M. , Jr.,  Pinto,  M.  R.  , and 

Smith,  R.  K.,  "Computation  of  steady-state  CMOS 
latch-up  characteristics,"  IEEE  Trans.  Computer- 
Aided  Design,  vol.  7,  pp . 307-323,  Feb.  1988. 

[Dar86]  Darwish,  M.  N.,  "Study  of  the  quasi-saturation 
effect  in  VDMOS  transistors,"  IEEE  Trans.  Electron 
Devices,  vol.  ED-33,  pp.  1710-1716,  Nov.  1986. 


159 


160 


[Darw84] 

[Darw87] 

[Fos86] 

[Fos88 ] 

[Foss86] 

[Goo83] 

[Har84 ] 
[Hef 86] 

[Huc82 ] 
[Huc84 ] 
[Kuo86] 

[ Lam7  0 ] 


Darwish,  M.  N.  and  Board,  K.,  "Lateral  resurfed 
COMFET,"  Electron.  Lett.,  vol . 20,  pp.  519-520, 

June  1984. 

Darwish,  M.  N.  and  Shibib,  M.  A.,  "DC  and  transient 
analysis  of  lateral  insulated  gate  conductivity 
modulated  transistors,"  in  Proc.  ECS  Symp.  on  High 
Voltage  and  Smart  Power  Devices,  pp.  295-309,  1987. 

Fossum,  J.  G.  and  McDonald,  R.  J.,  "Charge-control 
analysis  of  the  COMFET  turn-off  transient,"  IEEE 
Trans.  Electron  Devices,  vol.  ED-33,  pp . 1377-1382 , 
Sept . 1986 . 

Fossum,  J.  G.,  McDonald,  R.  J.,  and  Shibib,  M.  A., 
"Network  representations  of  LIGBT  structures  for  CAD 
of  power  integrated  circuits,"  IEEE  Trans.  Electron 
Devices,  vol.  35,  pp.  507-515,  Apr.  1988. 

Fossum,  J.  G.  and  Veeraraghavan,  S.,  "Partitioned- 
charge-based  modeling  of  bipolar  transistors  for 
non-quasi-static  simulation, " IEEE  Electron  Device 
Lett . , vol.  EDL-7 , pp.  652-654";  DecT  1986 . 

Goodman,  A.  M.  , Russell,  J.  P.,  Goodman,  L. 
A.,  Nuese,  C.  J.,  and  Neilson,  J.  M.,  "Improved 
COMFETs  with  fast  switching  speed  and  high-current 
capability,"  I EDM  Tech.  Dig.,  pp . 79-82,  1983. 

Harris  Semiconductor  Corp.,  SLICE  Manual  Rev.  4.08, 
Melbourne,  FL,  Jan.  1984. 

Hefner,  A.  R.,  Jr.  and  Blackburn,  D.  L., 

"Performance  trade-off  for  the  insulated  gate 
bipolar  transistor:  buffer  layer  versus  base 

lifetime  reduction,"  presented  at  IEEE  Power 
Electronics  Specialists  Conference,  1986. 

Hu,  C.  and  Chi,  M.-H.,  "Second  breakdown  of  vertical 
power  MOSFET's,"  IEEE  Trans.  Electron  Devices,  vol 
ED-29,  pp.  1287-1293,  Aug.  TWT. 

Hu,  C .,  Chi,  M.-H.,  and  Patel,  V.  M.,  "Optimal  design 
of  power  MOSFET's,"  IEEE  Trans.  Electron  Devices, 
vol.  ED-31,  pp.  1693-1700,  Dec.  mi; 

Kuo,  D.-S.  and  Hu,  C.,  "Optimization  of  epitaxial 
layers  for  power  bipolar-MOS  transistor, " IEEE 
Electron  Device  Lett . , vol.  EDL-7,  pp.  510-5T27 


Lampert,  M.  A.  and  Mark,  P.,  Current  Infection  in 
Solids,  New  York:  Academic  Press";  1970 . 


161 


[Luc88 ] Lu,  C.-Y.,  Tsai,  N.-S.,  Dunn,  C.  N.,  Riffe,  P. 

C.,  Shibib,  M.  A.,  Furnanage,  R.  A.,  and  Goodwin, 
C.  A.,  "An  analog/digital  BCDMOS  technology  with 
dielectric  isolation-devices  and  processes,"  IEEE 
Trans.  Electron  Devices,  vol . 35,  pp . 230-239,  Feb. 
1988 . “ 

[McD85]  McDonald,  R.  J.,  Fossum,  J.  G.,  and  Shibib,  M.  A., 
"A  physical  model  for  the  conductance  of  gated  p-i-n 
switches,"  IEEE  Trans.  Electron  Devices,  vol.  ED-32, 
pp.  1314-1320,  July  1985. 

[Men86]  Mena,  J.  G.  and  Salama,  C.  A.  T.,  "High-voltage 
multiple-resistivity  drift-region  LDMOS , " Solid-St 
Electron. , vol.  29,  pp.  647-656,  1986.  

[Mul77 ] Muller,  R.  S.  and  Kamins,  T.  I.,  Device  Electronics 
for  Integrated  Circuits . New  York:  Wiley,  1977 . 

[Nak85]  Nakagawa,  A.  and  Ohashi,  H.,  "600-  and  1200-V 

bipolar-mode  MOSFET' s with  high-current  capability, " 
IEEE  Electron  Device  Lett . , vol.  EDL-6,  pp . 378-380, 
July  1985. 

[Pat86 ] Pattanayak,  D.  N.,  Robinson,  A.  L.,  Chow,  T.  P., 
Adler,  M.  S.,  Baliga,  B.  J.,  and  Wildi,  E.  J.,  "n- 

channel  lateral  insulated  gate  transistors:  Part  I- 
Steady-state  characteristics,"  IEEE  Trans.  Electron 
Devices,  vol.  ED-33,  pp . 1956-1963,  Dec.  1986. 

[Pin84]  Pinto,  M.  R.,  Rafferty,  C.  S.,  and  Dutton,  R. 

W.,  "PISCES-II  user's  manual,"  Stanford  Electronics 
Labs.,  Stanford  University,  Stanford,  CA,  1984. 

[Poc76 ] Pocha,  M.  D.  and  Dutton,  R.  w.,  "A  computer-aided 
design  model  for  high-voltage  double  diffused  MOS 
(DMOS ) transistors,"  IEEE  J.  Solid-State  Circuits, 
vol.  SC-11,  pp.  718-726,  Oct.  V5TZ~. 

[Rob86]  Robinson,  A.  L.,  Pattanayak,  D.  N.  , Adler,  M. 

S.,  Baliga,  B.  J.,  and  Wildi,  E.  J. , "Lateral 
insulated  gate  transistors  with  improved  latching 
characteristics,  " IEEE  Electron  Device  Lett. , vol 
EDL-7,  pp . 61-63,  Feb.  1986 . — 

[Rus83 ] Russell,  J.  P.,  Goodman,  A.  M.,  Goodman,  L.  A.,  and 
Neilson,  J.  M.,  "The  COMFET-A  new  high  conductance 
MOS-gated  device,"  IEEE  Electron  Device  Lett.,  vol. 
EDL-4 , pp.  63-65,  Mar . 1983.  ‘ 

Seitchik,  J.  A.,  Chatter  jee,  A.,  and  Yang,  P., 
An  analytic  model  of  holding  voltage  for  latch- 
up  _ in  epitaxial  CMOS,"  IEEE  Electron  Device  T.Atl- 
vol.  EDL-8,  pp.  157-159,  Apr.  1^87. ~ 


[Sei87] 


[Sod84 ] 

[Sze81] 
[Vee86 ] 

[W1187] 


162 

Sodini,  C.  G.,  Ko,  P.,  and  Moll,  J.  L.,  "The 
effect  of  high  fields  on  MOS  device  and  circuit 
performance,"  IEEE  Trans.  Electron  Devices.  vol.ED- 
31,  PP.  1386-T393,  Oct. IJWT. 


Sze,  S.  M.,  Physics  of  Semiconductor  Devices,  2nd 
ed.  New  York!  Wiley,  1961 . 


Veeraraghavan,  S.,  Fossum,  J.  G.,  and  Eisenstadt, 
W.  R.,  "SPICE  simulation  of  SOI  MOSFET  integrated 
circuits,"  IEEE  Trans.  Computer-Aided  Desian.  vol 
cad-5,  pp.  "653-638,  oct  1986. 3~ 


Williams,  R.  K.,  Timmes,  F.  X.,  Busse,  R.  W., 
and  Siu,  I.  K.,  "Two-dimensional  device  simulation 
of  high-voltage  lateral  DMOST's,"  in  Proc.  ECS 

318P327n  H1987  V°ltaqe  and  Smart  Power  Devices,  pp. 


BIOGRAPHICAL  SKETCH 


Yeong-Seuk  Kim  was  born  in  Kyung-nam,  Korea,  on  July  3, 
1957.  He  received  the  B.S.  and  M.S.  degrees  in  electronics 
from  Seoul  National  University,  Korea,  in  1980  and  1982, 
respectively.  Since  October  1985,  he  has  been  working  toward 
the  Ph.D.  degree  in  electrical  engineering  at  the  University 
of  Florida,  Gainesville.  His  doctoral  research  involves 
modeling  of  semiconductor  power  devices. 

From  1982  to  1985,  he  was  with  the  Central  Research 
Laboratories  of  Gold  Star,  Seoul,  Korea,  where  he  designed 
integrated  circuits  for  consumer  electronics . 


163 


I certify  that  I have  read  this  study  and  that  in  my 
opinion  it  conforms  to  acceptable  standards  of  scholarly 
presentation  and  is  fully  adequate,  in  scope  and  quality,  as 
a dissertation  for  the a degree  of  Dpctor  of  Philosophy. 


Jtfrry  C/.  Fossum,  Chairman 
rofessor  of  Electrical  Engineering 


I certify  that  I have  read  this  study  and  that  in  my 
opinion  it  conforms  to  acceptable  standards  of  scholarly 
presentation  and  is  fully  adequate,  in  scope  and  quality,  as 

degree  of  Doctor  of  Philosophy. 


a dissertation  for  the 


7L 

Sheng  S . Li'" 

Professor  of  Electrical  Engineering 


I certify  that  I have  read  this  study  and  that  in  my 
opinion  it  confopns  to  acceptable  standards  of  scholarly 
presentation  and  is  fully  adequate,  in  scope  and  quality,  as 
a dissertation  for  the  degree  of  Doctor  of  Philosophy. 

IzL.  

Dorothea  E . Burk 

Professor  of  Electrical  Engineering 

1 certify  that  I have  read  this  study  and  that  in  ay 
opinion  it  confopns  to  acceptable  standards  of  scholarly 
presentation  and  is  fully  adequate,  in  scope  and  quality,  as 
a dissertation  for  the  degree  of  Doctor  of  Philosophy. 

Li/luv^  

Khai  D . T . Ngo  ^ 

Assistant  Professor  of  Electrical 
Engineering 


I certify  that  I have  read  this  study  and  that  in  my 
opinion  it  conforms  to  acceptable  standards  of  scholarly 
presentation  and  is  fully  adequate,  in  scope  and  quality,  as 
a dissertation  for  the  degree  of  Do^r  of  Philosophy. 

5 


Kevin  S . Jones 

Assistant  PptJfessor  of  Materials  Science 
and  Engineering 


This  dissertation  was  submitted  to  the  Graduate  Faculty 
of  the  College  of  Engineering  and  to  the  Graduate  School  and 
was  accepted  as  partial  fulfillment  of  the  requirements  for 
the  degree  of  Doctor  of  Philosophy. 

August  1990 

Winfred  M.  Phillips 
Dean,  College  of  Engineering 


Madelyn  M.  Lockhart 
Dean,  Graduate  School 


